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VORWORT 


Dieses Buch ist für alle geschrieben, die einen Commodore 64 besitzen und gern 
etwas mehr über einige mathematische Verfahren wissen möchten. Wahrscheinlich 
wissen Sie, was für ein Programm Sie schreiben wollen, sind aber vielleicht nicht 
ganz sicher, was Sie an Mathematik dazu brauchen. Können Sie etwas mit COS, 
ABS oder SGN anfangen? 

Hier sind alle mathematischen Funktionen des Commodore 64 beschrieben. Ihre 
Anwendung ist in kurzen Programmen erläutert, die Sie wörtlich übernehmen und in 
Ihre eigenen Programme einbauen können. 

Aber das Buch ist nicht nur eine Einführung in diese elementaren mathematischen 
Funktionen. Es enthält auch Einführungen und Programme zu so unterschiedlichen 
Gebieten wie Codes und Kryptographie, Zufallszahlen, Trigonometrie, Primzahlen 
und Auswertung statistischer Daten. Diese Informationen können Sie sowohl für 
ernsthaftes Programmieren wie für Spiele verwenden. 

Mein Dank gilt Ann, Kora und Inga für ihre Geduld beim Schreiben dieses Buches. 

Czes Kosniowski 
Newcastle upon Tyne, September 1983 
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BEMERKUNGEN ZU DEN PROGRAMMEN 

Der Commodore 64 verwendet ‘Steuerzeichen’, um Funktionen wie Cursorbewe¬ 
gung und Zeichenfarbe zu kontrollieren. Diese Steuerzeichen werden gewöhnlich 
revers dargestellt. Ein revers dargestelltes Herz in Anführungszeichen löscht zum 
Beispiel den Bildschirm und bringt den Cursor in die Home-Position. Um Schwierig¬ 
keiten bei den Programmlistings zu vermeiden, sind keine Steuerzeichen benutzt 
worden. Statt dessen sind die zugehörigen CHR$-Codes verwendet worden. Die 
gebräuchlicheren sind in der folgenden Liste zusammengestellt. Weitere sind im 
Anfang F des Commodore 64 Handbuchs enthalten. 

CHR$(5) Weiß 

CHR$(17) Cursor nach unten 

CHR$(28) Rot 

CHR$(30) Grün 

CHR$(31) Blau 

CHR$(145) Cursor nach oben 

CHR$(147) Bildschirm löschen und Cursor Home 

CHR$(154) Hellblau 

CHR$(157) Cursor nach links 

CHR$(158) Gelb 

CHR$(159) Türkis 

Beim Eingeben der Programme können Sie die CHR$-Codes durch Steuerzeichen 
ersetzen. 
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KAPITEL 1 


EINFACHE FUNKTIONEN 


ÜBERSICHTLICHE DARSTELLUNG VON ZAHLEN 

Ganze Zahlen sind Zahlen, die keine Stellen nach dem Dezimalpunkt haben. Wenn 
Sie die folgenden Zeilen in einem Programm verwenden, erscheinen ganze Zahlen 
übersichtlich angeordnet auf Ihrem Bildschirm. 

L = LEN(STR$(X)) 

PRINT TAB(25-L) X 

Die Funktion STR$(X) wandelt die Zahl X in einen String um, LEN berechnet seine 
Länge, und TAB führt den Cursor an die richtige Stelle des Bildschirms. 

9 

123 
- 10 
89 

Bei nicht-ganzen Zahlen mißlingt diese Darstellung. 

89 

1.2 

-13.89 

.126 

Die Zahlen sind rechtsbündig; es wäre aber schöner, wenn die Dezimalpunkte 
untereinander ständen. Das erreicht man mit den Funktionen INT(X) und ABS(X). 
Die Funktion INT(X) bestimmt den ganzzahligen Anteil von X, also die größte ganze 
Zahl, die kleiner oder gleich X ist. Zum Beispiel 


INT(1.21) = 1 
INT(2) = 2 
| NT(2.1) = 2 
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INT(—2) = -2 
INT(—2.1) = -3 
INT(9.1) = 9 
INT(—9.2) = -10 


Die Funktion ABS(X) bestimmt den absoluten Wert von X, d. h. die Zahl ohne das 
Vorzeichen + oder Zum Beispiel 

ABS(9.1) = 9.1 
ABS(-9.1) = 9.1 

Das folgende Programm stellt beliebige Zahlen übersichtlich auf Ihrem Bildschirm 
dar: 

Y = INT(ABS(X)) : L = LEN(STR$(Y)) 

IF Y = 0 AND X <> 0 THEN L = L-1 
PRINT TAB (25—L) X 


Eine typische Darstellung sieht so aus: 
3 

.23 

-89.14 

6712399.1 

2.23871 

- 1.22 

-.13 


In der ersten Zeile des Programms kümmert sich ABS um die negativen Zahlen und 
INT um die nicht-ganzen Zahlen. Beachten Sie, daß die Funktion INT allein nicht 
dasselbe leistet. (Nehmen Sie etwa die Zahl -9.1) Die zweite Zeile behandelt 
Zahlen, die größer als -1 und kleiner als 1 sind. 

Das obige Programm zeigt eine einfache Anwendung der Funktionen INT und ABS. 
Es funktioniert nicht für Zahlen nahe 0 (absoluter Wert kleiner oder gleich 0.01) und 
sehr große Zahlen (absoluter Wert größer oder gleich 1000000000). Ansonsten 
mißglückt die Darstellung, wenn die wissenschaftliche Notation E vorkommt. 

3 

-89.14 

1E-04 

9.9E+16 
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Wenn Sie wollen, können Sie unserem Programm zwei Zeilen hinzufügen, um 
Zahlen zu behandeln, die in der wissenschaftlichen Notation E angezeigt sind. 


RUNDEN VON ZAHLEN 

Die Funktion INT läßt sich zum ‘Runden’ von Zahlen gebrauchen. Wenn Sie 
beispielsweise DM 565,58 auf einem Bankkonto haben und 9% Zinsen im Jahr 
bekommen, beträgt Ihr Guthaben nach einem Jahr 

565.58 + 565.58*9/100 

Mit Ihrem Commodore 64 können Sie ausrechnen, daß das 616,4822 ergibt. Die 
Bank würde diesen Wert natürlich auf DM 616,48 abrunden. Entsprechend würde 
etwa ein Betrag von 76,6752 auf DM 76,68 au/gerundet. Dieses Runden bewerk¬ 
stelligt Ihr Commodore 64 mit der Zeile 

X = INT(X*100 + 0.5)/100 

Zur Umwandlung in Pfennig wird X erst mit 100 multipliziert. Das Hinzuaddieren von 
0,5 bewirkt dann ein Aufrunden, wenn die Pfennig-Bruchteile mindestens Einhalb 
ausmachen. Die Funktion INT ignoriert alle Bruchteile, und die Division durch 100 
wandelt die Zahl wieder in DM um. 

Allgemein ergibt die Programmzeile 


B = INT(A*10|D+0,5)/10|D 
den Wert von A, gerundet auf D Dezimalstellen. 


BANKKONTEN 

Das oben ausgelistete Programm zur übersichtlichen Darstellung von Zahlen auf 
em Bildschirm läßt sich zum Beispiel für ein Saldierprogramm verwenden. Auf 
diese Weise erhält man etwa folgende Darstellung: 
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ART AUSGABEN 

ÜBERTRAG 

869162 46.22 

869164 169 

869165 15.01 

SCHECK 

Wie Sie sehen, fehlt die Null, wenn der Betrag keine Pfennige oder ein Vielfaches 
von zehn Pfennig enthält. Um das zu ändern, müssen wir unsere Zahlen in STRings 
umwandeln und die nötige Anzahl von Nullen hinzufügen. Nachfolgende Pro¬ 
grammzeilen zeigen, wie man das macht. 

X$ = STR$(X) : L = LEN(STR$(X*100)) : M = LEN(X$) 

IF M = L THEN X$ = X$ + “0“ 

IF M = L-2 THEN X$ = X$ + “.00“ 

PRINT TAB(25-L) X$ 


EINNAHMEN SALDO 

596.61 

550.39 

381.39 
366.38 

75.7 442.08 


Beachten Sie, daß wir die Funktionen INT und ABS nicht brauchen. Der eben 
dargestellte Kontoauszug würde nun wie folgt aussehen: 


ART 

ÜBERTRAG 

AUSGABEN 

EINNAHMEN 

869162 

46.22 


869164 

169.00 


869165 

15.01 


SCHECK 


75.70 


SALDO 

596.61 

550.39 

381.39 
366.38 
442.08 


ÜBERZOGENE KONTEN 

Bankkonten werden gelegentlich überzogen (oder geraten in die roten Zahlen). Das 
geschieht, wenn der Saldo negativ (kleiner als Null) wird. So bedeutet ein Saldo von 
-DM 64,00, daß Sie DM 64,00 überzogen haben. Die folgenden Programmzeilen 
zeigen, wie man mit Hilfe der Funktion ABS Kontoauszüge ausdrucken und 
überzogene Beträge kennzeichnen kann (mit DB für Debit). 

X$ = STR$(ABS(X)) : L = LEN(STR$(X*100)) : M = LEN(X$) 

IF M = L THEN X$ = X$ + “0“ 

IF M = L-2 THEN X$ = X$ + “.00“ 

PRINT TAB(25-L) X$; 
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IF X < 0 THEN PRINT „DB“; 
PRINT 


Zum Beispiel: 

ART AUSGABEN EINNAHMEN 

ÜBERTRAG 

869166 52.80 

869167 422.00 - 


SALDO 
442.08 
389.29 
32.72 DB 


FARBIGE KONTOAUSZÜGE 

Die Funktion SGN(X) ist die Vorzeichen-Funktion, die das Vorzeichen (positiv, 
negativ oder Null) der Zahl X ergibt. Das Ergebnis ist +1 für eine positive Zahl, -1 
für eine negative Zahl und 0 für die Zahl 0. Zum Beispiel: 

SGN(9.21) = 1 
SGN(-9.1) = -1 
SGN(0) = 0 

Die Funktion SGN wird insbesondere verwendet, wenn das Programm verschie¬ 
dene Unterprogramme abarbeiten soll, je nachdem, ob das Vorzeichen einer Zahl 
positiv, negativ oder Null ist. So würde das Programm aufgrund der Programmzeile 

ON SGN(X) + 2 GOSUB 1000, 1100, 1200 

das Unterprogramm 1000 ausführen, falls X negativ ist, das Unterprogramm 1100, 
falls X gleich Null ist, und das Unterprogramm 1200, falls X positiv ist. 

Eine interessante Anwendung der Funktion SGN besteht darin, die Zeichenfarbe 
auf einfache Weise zu ändern. Auf dem Commodore 64 bedeutet CHR$(28) rot und 
CHR$(30) grün. Also ist CHR$(29+SGN(X)) rot oder grün, je nachdem, ob X 
negativ oder positiv ist. Das nächste Programm ergänzt unser Saldierprogramm um 
diesen Farbwechsel. Das erste POKE wählt eine Hintergrundfarbe, vor der rot und 
grün deutlich sichtbar sind. 
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POKE 53281,7 

X$ STR$(ABS(X)) : L = LEN(STR$(X*100)) : M = LEN(X$) 

IF M = L THEN X$ = X$ + “0“ 

IF M = L-2 THEN X$ = X$ + “.00“ 

PRINT CHR$(29+SGN(X)) TAB(25-L) X$; 

IF X < OTHEN PRINT “DB“; 

PRINT CHR$(30) 

Dieses kurze Programm druckt Zahlen mit zwei Dezimalstellen (zum Beispiel DM 
und Pfennig) senkrecht untereinander aus. Dabei wird der ABSolute Wert einer 
negativen Zahl rot und mit DB dahinter ausgedruckt. 
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KAPITEL 2 


TRIGONOMETRIE 


MASSSTABTREUE ZEICHNUNGEN 

Nur selten kann man die Höhe großer Gebäude, Berge, Bäume usw. direkt messen. 
Man kann aber die Höhe eines Gebäudes oder Baumes bestimmen, indem man aus 
einiger Entfernung den Winkel zwischen der Horizontalen und dem höchsten Punkt 
des Objekts mißt (mit einem Klinometer, das im Grunde nur ein veredelter Winkel¬ 
messer ist) und dann den Abstand bestimmt, den man zum Objekt hat Siehe 
Abbildung 1. 



Zur Bestimmung der Höhe einer kilometerweit entfernten und in Wolken gehüllten 
Bergspitze könnten Sie dieses Verfahren nicht anwenden. Die Wolken würden 
stören, und der horizontale Abstand könnte nicht gemessen werden. In solchen 
allen würde ein Instrument wie das Tellurometer weiterhelfen. Es lokalisiert die 
Bergspitze mit Hilfe von Radar. Außerdem mißt es Winkel und Abstand zwischen 
hnen und der Spitze. Mit einer maßstabtreuen Zeichnung können Sie dann die 
Hohe des Berges bestimmen. Siehe Abbildung 2. 

Ein weiteres Beispiel ist die Bestimmung der Breite eines großen Teiches oder 
Gees; siehe Abbildung 3. 
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Abbildung 2. 



Abbildung 3. 


Mit den gemessenen Werten könnte man eine maßstabtreue Zeichnung anteiligen 
und die fragliche Länge schätzen. 

Hier ist ein ähnliches Beispiel: Ein Navigator befindet sich in einer bestimmten 
Position A. Er ist 150 km westlich von der Stadt B und 188 km von der Stadt C 
entfernt. Von seiner Position aus ist der Winkel zwischen beiden Städten 23 Grad. 
Wie weit sind die beiden Städte voneinander entfernt? Eine maßstabtreue Zeich¬ 
nung könnte auch diese Aufgabe lösen. 
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Maßstabtreue Zeichnungen beantworten zwar die oben erwähnten Fragen, aber sie 
liefern nur eine große Näherung. Sie sind weder genau noch immer praktikabel. 
Einen anderen Zugang bietet die Trigonometrie, und dabei hilft Ihnen Ihr Commo- 
dore 64. 


DIE TRIGONOMETRISCHEN FUNKTIONEN 

Die drei wichtigsten trigonometrischen Funktionen sind SIN (die Sinus-Funktion), 
COS (die Cosinus-Funktion) und TAN (die Tangens-Funktion). Alle stellen Verhält¬ 
nisse verschiedener Seiten eines rechtwinkligen Dreiecks dar. Das unten gezeigte 
Dreieck ist ein solches rechtwinkliges Dreieck. Der Winkel an der linken Ecke ist mit 
dem Symbol X bezeichnet. Die drei Seiten des Dreiecks heißen Ankathete von X, 
Gegenkathete von X und Hypootenuse (längste Seite). 

Gegenkathete 
Ankathete 

Gegenkathete 
Hypotenuse 

Ankathete 
Hypotenuse 


TAN(X) = 
SIN(X) = 
COS(X) = 


Es ist sinnvoll, sich die folgenden Werte zu merken: 


SIN(0°) = 0 
SIN(30°) = 0.5 
SIN(45°) = 1/SQR(2) 
SIN(60°) = SQR(3)/2 
SIN(90°) = 1 


COS(0°) = 1 
COS(30°) = SQR(3)/2 
COS(45°) = 1/SQR(2) 
COS(60°) = 1/2 
COS(98°) = 0 


TAN(0°) = 0 
TAN (30°) = 1/SQR(3) 
TAN (45°) = 1 
TAN(60°) = SQR(3) 


Kennt man den Winkel X und eine der drei Seiten eines rechtwinkligen Dreiecks, so 
kann man die anderen beiden Seiten bestimmen. Sind beispielsweise der Winkel X 
und die Länge der Ankathete bekannt, so ergeben sich die beiden anderen Seiten 
aus folgenden Formeln: 


Gegenkathete = TAN(X) * Ankathete 
Hypotenuse = Ankathete / COS(X) 
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GEGENKATHETE 



Abbildung 4. 

Man kann die trigonometrischen Funktionen auch anhand eines Kreises mit dem 
Radius 1 beschreiben. Dazu zeichnet man den gegebenen Winkel wie in Abbil¬ 
dung 5 ein. Die Werte der verschiedenen trigonometrischen Funktionen lassen sich 
aus der Abbildung ablesen. 

Mathematisch werden Längen horizontal von links nach rechts und vertikal von 
unten nach oben gemessen. Das erklärt, warum z. B. COS(X) in Abbildung 6 einen 
negativen Wert hat. 

Sie erhalten SIN, COS und TAN eines Winkels X durch Eingeben von 

PRINT SIN(X) etc., 

wobei für X der gewünschte Wert einzusetzen ist. Das einzige Problem dabei ist, 
daß der Commodore 64 wie die meisten Mikrocomputer die Winkel im Bogenmaß, 
nicht in Grad erwartet. Glücklicherweise kann man Grad leicht in Radiant (engl.: 
radian) umrechnen und umgekehrt. 

Zunächst: Was ist ein Radiant? Zeichnen Sie einen Kreis vom Radius 1. Tragen Sie 
auf dem Umfang eine Strecke von der Länge des Radius ab. Der von diesem Bogen 
aufgespannte Winkel ist 1 Radiant oder ungefähr 57°. Siehe Abbildung 7. 
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1 RADIANT 


Abbildung 7. 


Die Zahl n (oder pi) ist bemerkenswert und berühmt. Sie ist definiert als das 
Verhältnis des Umfangs eines Kreises zu seinem Durchmesser. Der (ungefähre) 
Wert von tt ist in Ihrem Commodore 64 gespeichert. Tippen Sie einfach 

PRINT n 

ein, und Sie erhalten den gespeicherten Wert. In einem Kreis vom Radius 1 ist der 
Durchmesser gleich 2. Also ist der Umfang des Kreises gleich 2*tt, und deshalb hat 
ein ganzer Kreis 2*n- Radiant. Da ein ganzer Kreis 360 Grad hat, ergibt sich 

360 Grad = 2*zr Radiant, 

180 Grad = n Radiant. 

Mit folgenden einfachen Formeln läßt sich Grad in Radiant und Radiant in Grad 
umrechnen: 

X Grad = XV/180 Radiant, 

Y Radiant = Y*180/rr Grad. 

Mit dem nächsten Programm können Sie die Seiten eines rechtwinkligen Dreiecks 
bestimmen. Sie müssen dazu einen Winkel und die Länge einer Seite eingeben. 
Das Programm berechnet dann die beiden anderen Seiten. 
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10 REM 1RECHTWINKL I GE DREIECKE 
E0 POKE53E80,6sPR INT CHR$( 147>SPRINT" 

»♦♦RECHTWINKLIGE DRE IECKE***"CHR*< 17> 
30 PR INT "DIESES PROGRAMM ERMOEGLICHT ES, 
SEITEN" 


40 PR INT"EINES RECHTWINKLIGEN DREIECKS 2 
U" 


50 PR INT"FINDEN, ANGENOMMEN SIE KENNEN E 
INE SEITE UND EINEN WINKEL"CHR*< 1 7 > 


PRINT" 
PRINT" 
PRINT" 
PRINT“ 
PR INT" 
75 PRINT" 
80 PRINT" 
PRINT" 
PRINT" 
PRINT" 
96 PRINT" 
100 


54 

55 
G0 
65 
70 


35 

90 

95 


* 


* " 
** " 
♦ *" 
♦ 11 
* " 


♦GEGENKATHETE" 

^ II 


HYPOTENUSE* 

* 

* *" 

♦WINKEL *" 

********** " 

ANKATHETE" CHR*<17> 

REM EINGABE DETAILS 

INPUT "WINKEL IN GRAD! ";X 

IF X <=0 OR X > =90THEN PRINT "FEHLER-K 


1 10 
120 

EIN DREIECK" .'GOTO 1 10 
125 
130 


PRINTCHR*< 147 > 

PRINT CHR* < 17 >"WELCHE 


NNT? " 


SEITE IST BEKA 


135 PRINT" 1 ( GEGENKATHETE? 

140 PRINT" 2 < ANKATHETE) 

150 PRINT" 3 < HYPOTENUSE)“CHR$< 

17 > 

155 INPUT " 1,2 ODER 3 ";T 
160 IF T< 1 OR T>3 OR T<>INT<T) THEN 155 
170 PRINTCHR$< 17)"GEBEN SIE DIE LAENGE D 
ER SEITE EIN." 

180 INPUT "LAENGE: ";L 

190 IF L <=0 THEN130 

£00 REM UMWANDLUNG IN BOGENMASS" 

£10 X =X * */ 180 

££0 REM AUFTEILEN 

£30 ON T GOSUB 300,350,400 
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NOCH EIN 


240 PRINTCHR:T< 17> "DAS WAR'S - 
VERSUCH? <J/N)s" 

250 GET G$ s IF G$ = " " THEN250 
260 IF Gf = "J"THEN RUN 
270 IF G$="N"THEN PRINT CHR*<147> " AUF 

WIEDERSEHEN"s END 
300 REM GEGENKATHETE BEKANNT 
310 PRINT CHR*C 17 > "ANKATHETE :" L/TANC 
X) 

320 PR I NT ” HYPOTENUSE : " L/'S I NC X > 

330 RETURN 

350 REM ANKATHETE BEKANNT 

360 PRINTCHR4T 17>"GEGENKATHETE TANCX) 
*L 

370 PR INT"HYPOTENUSE !"L/COSCX> 

330 RETURN 

400 REM HYPOTENUSE BEKANNT 

410 PRINT CHR$C 17; 1 "GEGENKATHETE SINC 

X> *L 

420 PRINT "ANKATHETE L*COSCX> 

430 RETURN 


INVERSE FUNKTIONEN 

Angenommen, wir kennen die Seiten eines rechtwinkligen Dreiecks - können wir 
dann die verschiedenen Winkel bestimmen? Die Antwort ist ja, doch dazu benöti¬ 
gen wir die inversen trigonometrischen Funktionen. Ist ein Winkel X gegeben, so 
ergibt TAN(X) eine Zahl, den Tangens des Winkels X. Wenn umgekehrt eine Zahl N 
gegeben ist, gibt es einen Winkel, dessen Tangens gleich dieser Zahl ist. Einen 
solchen Winkel könnte man den inversen Tangens von N nennen. Gewöhnlich wird 
er mit ATN(N), das heißt Arcustangens von N, bezeichnet. 

Sehen Sie sich das Dreieck in Abbildung 8 an. 

Wenn wir die Werte von A und B kennen, können wir den Wert des Winkels X 
bestimmen. Wir wissen, daß TAN(X) = B/A ist, also X = ATN(B/A). Sie können die 
betreffenden Werte in diesen Ausdruck einsetzen und die Antwort von Ihrem 
Commodore 64 ausdrucken lassen. Die Antwort würde natürlich in Radiant angege¬ 
ben. Für eine Antwort in Grad müssen Sie das Ergebnis mit 180/w multiplizieren. 
Die trigonometrischen Funktionen SIN und COS haben ebenfalls inverse Funktio¬ 
nen; sie werden mit ASN (Arcussinus) und ACS (Arcuscosinus) bezeichnet. 
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Abbildung 8. 



ASN(N) ist der Winkel, dessen Sinus gleich N ist; entsprechend ist ACS(N) der 
Winkel, dessen Cosinus gleich N ist. Leider sind in Commodore-64-Basic wie bei 
vielen anderen Mikrocomputern diese Funktionen nicht vordefiniert. Sie können 
aber leicht aus der Funktion ATN abgeleitet werden. 

Um zu verstehen, wie man ASN aus ATN ableitet, betrachten Sie das rechtwinklige 
Dreieck mit Hypotenuse 1 in Abbildung 9. 
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Angenommen, wir kennen den Wert von B und wollen den Winkel X finden. Wir 
wissen, daß SIN(X) = B ist, also X = ASN(B), aber wie erwähnt ist ASN nicht im 
Commodore 64 vorhanden. Wenn uns aber der Wert von A bekannt wäre, könnten 
wir ATN verwenden, weil auch X = ATN(B/A) gilt. Um A zu bestimmen, benutzen 
wir den Satz des Pythagoras. 

Der Satz des Pythagoras ist Ihnen sicher bekannt. Wörtlich besagt er, daß das 
Quadrat der Hypotenuse eines rechtwinkligen Dreiecks gleich der Summe der 
Quadrate der beiden anderen Seiten ist. In Symbolen ausgedrückt: 


C 2 = A 2 + B 2 


Dabei bezeichnet C die Länge der Hypotenuse. Da unsere Hypotenuse die Länge 1 
hat, erhalten wir 

1 = A*A + B*B 

oder 

A*A = 1 - B*B 

also 

A = SQR(1 - B*B) 

Wegen X = ATN(B/A) ergibt sich 
X = ATN(B/SQR(1 - B*B)) 

Andererseits gilt auch ASN(B) = X und deshalb 
ASN(B) = ATN(B/SQR(1 - B*B)) 

Auf ähnliche Weise können wir folgende Formel für ACS(A) herleiten: 

ACS(A) = tt/2 - ATN(A/SQR(1 - A*A)) 

Beachten Sie, daß ACS(A) = tt/ 2 - ASN(A) ist. Anhang H des Commodore-64- 
Handbuchs enthält andere Beispiele mathematischer Funktionen, die nützlich sein 
können, aber in Commodore-64-Basic nicht vordefiniert sind. 
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NICHT-RECHTWINKLIGE DREIECKE 


Die ersten beiden Beispielaufgaben aus dem Abschnitt über maßstabtreue Zeich¬ 
nungen können mit dem Programm ‘Rechtwinklige Dreiecke’ gelöst werden. Im 
dritten Beispiel allerdings ist das Dreieck (im allgemeinen) nicht rechtwinklig. 

Ein Dreieck hat drei Winkel und drei Seiten. Wenn wir drei beliebige Werte davon 
kennen (mit Ausnahme von drei Winkeln), können wir die anderen Werte bestim¬ 
men. Nehmen wir an, wir kennen zwei Seiten und einen Winkel. Die dritte Seite und 
die beiden anderen Winkel können wir dann mit Hilfe einer Formel berechnen. 
Bezeichnen wir die drei Winkel unseres Dreiecks mit X, Y und Z, und die drei Seiten 
mit SX, SY und SZ, wobei die Seite SX dem Winkel X gegenüberliegt, etc. Siehe 
Abbildung 10. 



Abbildung 10. 


Die verschiedenen Seiten und Winkel stehen über folgende Gleichungen miteinan¬ 
der in Beziehung: 

Der Cosinus-Satz: 

SZ*SZ = SX*SX + SY'SY - 2*SX*SY*COS(Z) 

SY*SY = SX*SX + SZ*SZ - 2*SX*SZ*COS(Y) 

SX*SX = SY*SY + SZ*SZ - 2*SY*SZ*COS(X) 

Der Sinus-Satz: SIN(X)/SX = SIN(Y)/SY = SIN(Z)/SZ 

Beachten Sie, daß COS(Z) = 0 ist, wenn Z ein rechter Winkel ist (also 90 Grad hat). 
Die erste Gleichung hat dann folgende Form: 

SZ*SZ = SX*SX + SY*SY 
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Das ist gerade der Satz des Pythagoras. 

Das nächste Programm bestimmt die übrigen Winkel und Seiten, wenn eine der 
folgenden Angaben bekannt ist: 

Seite Seite Seite: Sie kennen alle drei Seiten und wollen die drei Winkel 
bestimmen. 

Seite Seite Winkel: Sie kennen zwei Seiten und einen Winkel, der nicht zwischen 
ihnen liegt (einen nicht eingeschlossenen Winkel), und wollen die dritte Seite und 
die übrigen Winkel bestimmen. 

Seite Winkel Seite: Sie kennen zwei Seiten und den Winkel dazwischen (den 
eingeschlossenen Winkel) und wollen die dritte Seite und die übrigen Winkel 
bestimmen. 

Seite Winkel Winkel: Sie kennen zwei Winkel und eine nicht dazwischen liegende 
Seite (eine nicht eingeschlossene Seite) und wollen die übrigen Seiten und den 
dritten Winkel bestimmen. 

Winkel Seite Winkel: Sie kennen zwei Winkel und die dazwischen liegende Seite 
(die eingeschlossene Seite) und wollen die übrigen Seiten und den dritten Winkel 
bestimmen. 

Beachten Sie, daß im zweiten Fall (Seite Seite Winkel) gewöhnlich zwei verschie¬ 
dene Dreiecke möglich sind, je nachdem, ob der der dritten Seite gegenüberlie¬ 
gende Winkel größer oder kleiner als 90 Grad ist. Abbildung 11 verdeutlicht das. 
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*** DREIECKE 


10 REM HDRE IECKE 
15 POKE53ES0,6 
£0 PRINTCHR*< 147> " 

*** "CHR$< 17 > 

30 PR INT"DIESES PROGRAMM FINDET DIE FEHL 
ENDEN " 

40 PR INT"SEI TEN UND WINKEL EINES DREIECK 
ES" 

45 PRINT 

50 PR INT"WELCHE INFORMATIONEN HABEN SIE? 

11 

55 PRINT 

60 PRINT"1) SSS s ALLE 3 SEITEN" 

70 PR INT"£) SSW : £ SEITEN UND NICHT EIN 
GESCHL. WINKEL" 

80 PR INT"3 > SWS : £ SEITEN UND EINGESCHL 
.WINKEL” 

30 PR INT"4) SWW : £ WINKEL UND NICHT EIN 
GESCHL. SEITE 

100 PR INT"5) WSW : £ WINKEL UND EINGESCH 
L.SEITE" 

110 REM WAHL 
115 PRINT 

1£0 INPUT "BITTE NUMMER EINGEBEN "; N 
1£5 IFN=1ORN=£QRN=40RN=30RN=5GOTO140 
130 PRINTCHR$< 17)"BITTE 1,£,3,4 ODER 5 V 
ERSUCHEN":GOTO1£0 

140 REM ARCUS - SINUS FUNKTION DEFENIERE 
N< AUF ZWEI STELLEN HINTERM KOMMA) 

150 DEF FNAS< K ) ® !NT< 18000*ATN< X/SGR< i-x* 

H) ) Pif + . 5 >.100 
1S0 REM AUFTEILEN 
165 PR. INTCHR*< 147 > 

170 PRINT50N N GOSUB 310,510,710,910,101 
0 

180 PR I NT:PR INTCHR$ < 17) "EIN NEUER START 
< J.''N) " 

130 GET G$:IF G$ = ""THEN190 
200 IF G$ = "J"THENRUN 

210 PRINTCHR$< 147):PR INT"AUF WIEDERSEHEN 
" : END 
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3 SEITEN BEKANNT 


300 REM 3 SEITEN 
310 PRINT" *** 

+ " 

320 M=1: GOSUB 1110!SH =S 
330 M=2:GOSUB 1110:SY=S 
340 M=3:GOSUB 1110:SZ=S 
350 A=< SY*SY+SZ *SZ-SX*SX>/< 2*SY*SZ) 

360 IF ABS< A > > = 1 THENPRINT"KE IN DREIECK" 
365 PR INT!PR INT 

370 PR INT"WINKEL GEGENUEBER SEITE 1 = "9 
0-FNAS': A > ! PR INT 

380 A = <SX*SX+SZ *SZ-SY*SY>/<2*SX*SZ> 

330 PR INT"WINKEL GEGENUEBER SEITE 2 = "9 
0-FNASC A >iPR INT 

400 A =< SK*SX+SY*SY-SZ *SZ >/< 2*SX*SY> 

4 10 PR INT"WINKEL GEGENUEBER SEITE 3 = "9 
0-FNASCA >!PR INT 
420 RETURN 

500 REM 2 SEITEN UND NICHT EINGSCHL. WIN 
KEL 

510 PRINT"** 2 SEITEN UND NICHT EINGSCHL 
.WINKEL **" 

520 PR INT"GEBEN SIE DIE SEITE EIN, VON D 
ER DER GEGENUEBERLIEGENDE WINKEL 

BE - " 

525 PRINT"KftNNT IST" 

530 M=l!GOSUB 1 110!SX=SiGOSUB 1210:AX=A 
540 M=2:GOSUB 1110:SY=S 

550 A=SIN<AX)*SY/SX: IF ABS<A)>1ORA=0THEN 
PR INT:PR I NT"KE IN DREIECK":RETURN 

559 PRINT 

560 PRINT"IST DER WINKEL GEGENUEBER SEIT 

E 2 GROE- SSER <>) ODER KLEINER << > "! 

565 PR INT"ALS 90 GRAD 7" 

570 INPUT"BITTE > ODER < EINGEBEN"; A$ 
580 IF A*< >"< " AND A*< >">"THEN570 
530 AY=FNAS< A YS IFA* = " > " AND AY<90 THENAY= 
90+AY 

600 PR INT!PR INT:PR INT "WINKEL GEGENU 

EBER SEITE 2 = "AY 
6 10 AZ = it -AX -AY* ■$/ 180 
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620 PRINTsPRINT"LAENGE VON SEITE 3 
= " SH * S I NC A2 > /S I NC AH > 

630 PR INTSPRINT"WINKEL GEGENUEBER SEITE 
3 = " INTC 18000*A2/rf+.5J/100 SPRINT 
640 RETURN 

700 REM 2 SEITEN UND EINGESCHL. WINKEL 
710 PRINT"*** 2 SEITEN UND EINGESCHL. WI 
NKEL ***" 

720 M=lseOSUB 1110s SH = S 
730 M=2 s GOSUB 1110:SY=S 
740 M=3sGOSUB 1210sA2=A 

750 S2= SQRC SH*SH+SY*SY-2*SH*SY*C0SCA2 )> 
760 IF S2 =0 THENPRINTSPRINT"KE IN DREIECK 
"sRETURN 

770 PRINTsPRINTS PRINT"LAENGE SEITE 2 
= " S2 

780 A = C SY*SY + S2 *S2 -SH*SH ~>/< 2*SY*S2 > 

730 PRINT-WINKEL GEGENUEBER SEITE 1 = "3 
0-FNASC A > 

800 A=< SH*SH +S2 *S2 -SYtSY) Y( 2*S*S2 ) 

810 PRINTCHRSC17)"WINKEL GEGENUEBER SEIT 
E £ = " 30-FNASC A > 

320 RETURN 

900 REM 2 WINKEL UND NICHT EINGESCHL. SE 
ITE 

910 PRINT"* £ WINKEL UND NICHT EINGESCHL 
. SEITE *" 

915 PRINT 

9£0 PR INT"GEBEN SIE DEN WINKEL, DER DER 
BEKANNTEN SEITE GEGENUEBER LIEGT 
925 PR INT"2UERST EIN." 

930 M= ls GOSUB 12 10 s AH=A s GOSUB U10sSH=S 
940 M=£sGOSUB 1210sAY=A 

950 A = <f-AH-AYs IFA< =0 THEN PR INT "KE IN DRE 
IECK"s RETURN 

960 PR INTSPRINTs PRINT"LAENGE SEITE £ 

= "SH*S IN( AY?/S INC AH ) 

970 PRINTs PR INT"WINKEL GEGENUEBER SEITE 
3 = "INT< 13000*A/rf+.5)/100 
380 PR INT:PRINT"LAENGE SEITE 3 
= " SH *S I N< A > /S I NC AH ) 
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990 RETURN 

1000 REM £ WINKEL UND EINGESCHL. SEITE 
1010 PRINT"**** £ WINKEL UND EINGEoCHL. 

SEITE ****" 

10£0 M=1 * GOSUB 1£ 10 5 AM =A 
1030 M=£•GOSUB 1£10:AY=A 
1040 M=3: GOSUB U10:SZ=S 

1050 fl^-AX-AY: IF A<0 THEN PRINT "KEIN D 
RE IECK":RETURN 

1060 PR INT:PR INT: PR INT"WINKEL GEGENUEBER 
SEITE 3 = "INTC 13000*A/rf+.5>/100CHR$C 17 

) 

1070 PR INT"LAENGE SEITE 1 
S2 *S I NC AX 5 /S I NC A ) CHR$ C 1 < ) 

1080 PR INT"LAENGE SEITE £ 

SZ *S INC AY)/S INC A > 

1090 RETURN 

1100 REM EINE SEITE BERECHNEN 
1110 PR INT* S=0:PR INT"LAENGE DER SEITE" M 
"EINGEBEN : ";:INPUT S 

1 i£0 IF s<=0 THEN PR INT"KE IN DREIECK":GO 

TO 1 1 10 
1130 RETURN 

1£00 rem EINEN WINKEL BERECHNEN 

12 10 prINT:A = 0:PR INT "WINKEL GEGENUEBER 

32IX£ 11 M " : " ' : INPUT A 

1230 IF A<=0.001 OR A >180 THEN PR INT"KE I 
N DREIECK”:GOTO 1£10 
1230 A=A*<i/190 : RETURN 


BRECHUNG 

Betrachtet man Gegenstände durch Glas oder Plastik hindurch, sehen sie oft 
verzerrt aus. Gewässer scheinen seichter als sie wirklich sind. Das liegt an der 
Brechung des Lichts. Tritt ein Lichtstrahl aus einem Medium (Luft) in ein anderes 
(Glas, Wasser. . .) ein, so wird er gebeugt oder gebrochen. Der Winkel, unter dem 
der Strahl auf das Glas eintrifft, heißt Einfallswinkel; der Winkel nach der Brechung 
heißt Brechungswinkel — siehe Abbildung 12. 
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Für ein gegebenes Material besteht eine feste Beziehung zwischen Einfallswinkel 
und Brechungswinkel. Sie wird durch das Snelliussche Gesetz ausgedrückt, das 
besagt, daß das Verhältnis der Sinuswerte für jedes Material (im Vergleich zur Luft) 
konstant ist. Dieses Verhältnis heißt Brechungsindex (Brechzahl). 

Brechungsindex = SIN (Einfallswinkel) 

SIN (Brechungswinkel) 

Der Brechungsindex für Glas ist ungefähr 1.5, für Wasser beträgt er 1.333, für 
Diamant ist er 2.417. 



Abbildung 12. 


Mit dem folgenden Programm können Sie den Brechungswinkel berechnen, wenn 
Sie den Einfallswinkel und den Brechungsindex kennen: 

10 REM 3BRECHUNG 

£0 POKE53S80,6sPRINTCHR$< 147)" *** 

** BRECHUNG *****" 

30 PRINTs PR INT"DIESES PROGRAMM ERMOEGLIC 
HT ES, DEN AUS-FALLSWIhKEL EINES LICHT"; 

40 PR INT“STRAHLS ZU BE- RECHNEN, WENN 
ER 
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AUF EIN ANDERES MEDIUM TRIFFT" 

100 REM EINGABE DETAILS 

110 PRINT:PRINT"BITTE EINFALLSWINKEL EIN 
GEBEN < IN GRAD > ’* 

1£0 INPUT “WINKEL ! "’X 

130 IF X<=0ORX>=90THENPRINT" !! IRRTUM 

! ! ":GOTO 120 

140 PR INT:PRINT"BITTE BRECHUNGSINDEX EIN 
GEBEN" 

145 INPUT "BRECHUNGSINDEX • "JR 

150 IF R<=0THENPRINT"!! IRRTUM M" :GOT 

0145 

160 REM UMRECHNUNG IN BOGENMASS 

170 X =X * <i/130 

180 REM BERECHNUNG 

190 Y=SIN<X)/R:Y=Y/SQR< 1-Y*Y) 

200 PRINTsPRINTJPRINT-AUSFALLSUINKELi "A 
TN 1 ' Y>*180/ir "GRAD" 

£10 PR INTSPRINT"PROZENTSATZ DES EINFALLS 
WINKELS : " INT< ATNC Y) * 100/X+ .5 ) 

£40 PR INT:PRINTsPRINT"NEUER START < J/N>? 

250 GET A$'• IF A$ = " " THEN£50 
£60 IF A$ = " J "THENRUN 

270 IF A* = "N"THENPRINTCHR$< 147 >"AUF WIED 
ERSEHEN”SEND 
£80 GOTO250 

TOTALREFLEXION 

Ein Stück Glas oder eine Wasseroberfläche verhalten sich gelegentlich wie ein 
gewöhnlicher Spiegel: Sie reflektieren alles. Das geschieht, wenn der Einfallswinkel 
zu groß ist und der Lichtstrahl voll zurückgeworfen wird. Der kleinste Winkel, bei 
dem dies auftritt, heißt Grenzwinkei der Totalreflexion. Er ergibt sich aus dieser 
einfachen Formel: 

SIN (Grenzwinkel) = 


Also läßt sich der Grenzwinkel mit Hilfe der oben beschriebenen Funktion ASN aus 
dem Brechungsindex bestimmen. 
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KAPITEL 3 


TRIGONOMETRIE DER ERDE 

DIE ERDE 

Die Erde ist fast eine Kugel. Nur am Nord- und Südpol ist sie leicht abgeflacht. Der 
Äquator hat einen Radius von 6378 Kilometern. Der Pol-Radius ist 6357 Kilometer. 
Da der Unterschied nur 0,3% beträgt, ist die Erde kaum von einer Kugel zu 
unterscheiden. Der mittlere Radius beträgt 6371 Kilometer. 

Wir sind daran gewöhnt, uns unter der kürzesten Verbindung zwischen zwei 
Punkten eine 'gerade Linie' vorzustellen. Das gilt aber nur für eine Ebene, bzw. für 
jede ‘flache’ Oberfläche. Auf einer Kugel wie der Erde ist die kürzeste Verbindung 
zwischen zwei Punkten ein Teil eines sogenannten Großkreises. Ein Großkreis ist 
ein Kreis, dessen Mittelpunkt im Mittelpunkt der Erde liegt. 

Großkreise durch den Nord- und Südpol heißen Längenkreise. Jedem Längenkreis 
ist ein Winkel zugeordnet. Der durch Greenwich in England verlaufende Längen¬ 
kreis wird mit 0° (0 Grad) bezeichnet. Die anderen ergeben sich jeweils aus dem 
Winkel zwischen Längenkreis, Erdmittelpunkt und Längenkreis durch Greenwich, 
gemessen am Äquator - siehe Abbildung 13. 

Gewöhnlich zählt man die geographische Länge von 0° bis 180°, sowohl nach Osten 
als auch nach Westen. 

An den Längengraden erkennt man, wie weit westlich oder östlich von Greenwich 
ein Punkt auf der Erde liegt. Wie weit nördlich oder südlich vom Äquator ein Punkt 
liegt, sieht man an den Breitenkreisen. Der Äquator hat den Breitengrad 0°. Alle 
Kreise auf der Erde, die parallel zum Äquator verlaufen, heißen Breitenkreise. 
Jedem Breitenkreis ist der Winkel zugeordnet, den er am Erdmittelpunkt mit dem 
Äquator bildet - siehe Abbildung 14. 

Die geographische Breite zählt von 0° bis 90° nach Norden und Süden. Der 
Breitengrad des Nordpols ist 90° Nord, der des Südpols 90° Süd. 

Jeder Punkt auf der Erde ist durch seinen Breiten- und Längengrad festgelegt. 
Beispielsweise liegt Newcastle upon Tyne (England) ungefähr 55° nördlicher Breite 
und 1,5° westlicher Länge. Genauere Werte sind 54° 58’ nördliche Breite und 1 ° 36’ 
westliche Länge, wobei ’ als Minuten gelesen wird und eine Minute 1/60 Grad 
beträgt. 

Breiten- und Längengrade bilden einen Satz von Koordinaten auf der Erdober¬ 
fläche. 
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Es ist nicht leicht, die (kürzeste) Entfernung (entlang eines Großkreises) zwischen 
zwei Punkten auf der Erde zu berechenn. Wie weit ist z, B. Newcastle upon Tyne 
von Paris entfernt? (Paris liegt ungefähr bei 49° nördlicher Breite, 2° östlicher 
Länge.) Mit einem geeigneten Programm stellt das auf dem Commodore 64 kein 
Problem dar. 

Das folgende Programm berechnet die Entfernung zwischen zwei Punkten auf der 
Erdoberfläche. Die Mathematik hinter diesem Programm beruht auf den Cosinus- 
und Sinus-Gesetzen, die im vorigen Kapitel behandelt wurden. 

10 REM 4ER0TRIGONOMETRIE 

£0 POKE53S80,6:PRINTCHRf< 147> " *** 

ERDTRIGONOMETRIE ***" 

£5 PRINT" <GROSSKREISNAVIGATION) " 

30 PR INT:PR INT"DIESES PROGRAMM BERECHNET 
DIE KUERZESTE ENTFERNUNG ZWISCHEN"; 
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40 PR INT"ZWEI PUNKTEN RUF DER ERDE" 

50 REM EINGABE DATEN 
60 FOR 1 = 1 TOS 

70 PR INT:PRINTCHR*< 144>"POS ITION"; I:PR IN 
T 

80 PRINTCHRfC 154 >:INPUT"LAENGE " ; A< I > 
90 IFA< IX 0 OR A< I>>30 THEN PRINT: PRINT 
"!! ZWISCHEN 0 UND 90 GRAD !!":GOTO80 
100 INPUT "N ODER S ";A$ 

110 IF A$ = "N"ORA$ = "S"THEN1£0 

115 PR INT:PR INT" ! ! NORD ODER SUED !! ":PR 

INT:GOTO 100 

120 IFA$ = "S"THENA( I> = -A< I > 

130 PRINT:INPUT"BREITE "JB(I> 

140 !FB( 1X0 OR B< I > > 1 80THENPR INT: PR INT" 

!! ZWISCHEN 0 UND 180 GRAD !!“:GOTO130 

150 INPUT "O ODER W ”A4 

160 IF A$ = "O" ORA$ = "W"THEN170 

165 PR INT:PR INT" ! ! OST ODER WEST ! ! ":PR I 

NT:GOTO 150 

170 IF A$ = "O " THENBC I>=-B< I > 

130 NEXT 

190 PRINT CHR$< 147>:PR INT:PR INT"ENTFERNU 
MG IN KILOMETERN ODER MEILEN ?" 

200 INPUT "K ODER M "JA$ 

210 IF A* = "K”0RA4 = "M"THEN220 
215 GOTO200 

220 R =6371:B* = "KILOMETER": IFA$ = "M"THENB$ 
= "MEILEN":R =3960 
230 REM BERECHNUNG 

240 A 1 =A< 1 >*if/180:A2=A< 2 > * it/I 80 : B 1 =B< 1 > * 
c/180 : B2=B< 2 > + tf/130 

250 B =ABS(B< 1 >-B< 2 >>: IFB>180THENB=180-B 
260 A=ABS<A'; 1 >-A< 2> >*«f/360:B=B*,t/360 
270 X = COS< A 1 ) *S IN< B > *COS< A2 > *S I N< B > +S IN< 
R > *S I N( A > 

280 D = 2tR*ATN< SQR< X/< 1 -X > > > 

290 PRINTCHRSC 17>CHR$< 17>CHR$< 144>“ENTFE 
RNUNG =“; INT<D * 100 +.5)/100; B$ 

300 PR INT:PR INT:PR INTCHR$( 154 >"NEUER STA 



THEN310 


RT <J/N> ?" 

3 10 GET s IFA*=” " 

320 IFA* = " J "THENRUN 

330 IFA* = "N"THENPR INTCHR*< 147>CHR*< 154 V 
AUF LJIEOERSEHEN"CHR*< 154 V END 
340 GOTO310 
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KAPITEL 4 


POTENZEN 


Das Quadrat einer Zahl X ist das Produkt der Zahl mit sich selbst, also X*X. Es wird 
mit X 2 oder Xt2 bezeichnet. Man nennt es auch die zweite Potenz von X. 

Die Potenzen von X sind entsprechende mehrfache Produkte der Zahl X mit sich 
selbst. Die erste, zweite, dritte, vierte und fünfte Potenz von 2 sind die Zahlen 

2, 4, 8, 16 und 32. 

Man bezeichnet diese Potenzen entsprechend mit 
2\ 2 2 , 2 3 , 2 4 und 2 5 . 

Mit dem folgenden einfachen Programm können Sie die Potenzen von 2 (oder jeder 
anderen Zahl) berechnen: 

10 REM 5P0TENZEN 

£0 POKE53£80,6 : PRINTCHR$( 147)CHR$<5>" 

*** POTENZEN *** M 

30 PR INT:PR INT"DIESES PROGRAMM ERRECHNET 
DIE POTENZEN EINER EINGEGEBENEN ZAHL" 

40 PRINT:INPUT "ZAHL EINGEBEN " ; X 
45 IFX=0THEN40 

50 PRINT:INPUT "WIEVIEL POTENZEN ";N 
60 IF N< 1ORN< >INT<N)THENPRINT:PR INT"NUR 
NATUERLICHE ZAHLEN BITTE":GOTO50 
70 REM BERECHNUNG 

30 PR INTCHR$< 147):PR INT sPR INT" I "X"tl" 

88 PRINT 
30 Y= 1 

100 FOR 1 = 1 TON:Y=Y#X:PR INT I ; Y:NEMT 
110 PR INT:PR INT:PR INT"NEUER START <J/N> 

? ” 

130 GET A$: IFA$ = ""THEN1H0 
130 IF A $ = "J"THENRUN 

140 IFA$ = "N"THENPRINTCHRS( 154)CHR$< 147)" 
ftlJ F WIEDERSEHEN" : END 
150 GOTO 1£0 
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Falls die Zahl N zu groß ist, wird OVERFLOW ERROR gemeldet. 

Mit einigen zusätzlichen Zeilen können Sie ein mögliches Überlaufen testen und 
vermeiden: 

10 REM 5P0TENZEN 1 

20 POKE53280,6 ! PRINTCHR$< 147>CHR$<5> " 

+** POTENZEN *** M 

30 PR INT s PR INT"DIESES PROGRAMM ERRECHNET 
DIE POTENZEN EINER EINGEGEBENEN ZAHL" 

40 PRINT!INPUT "ZAHL EINGEBEN "JX 
45 IFX =0THEN40 

50 PRINT:INPUT »WIEVIEL POTENZEN ";N 
60 IF N< 10RN< >INT<N)THENPRINTiPR INT"NUR 
NATUERLICHE ZAHLEN B ITTE " '• GOTO50 
6 1 IF<N+1> *LOG< ABS<X > X 12S*L0G< 2 JTHEN70 

62 PR INT!PRINT"RECHNERKAPAZITAET WIRD UE 
BERSCHRITTEN !" 

63 F'R INT: PR INT "NEU START < J/N> ?“ 

64 GETA$! IFA$ = ""THEN64 

65 IFA* = "N"THEN70 

66 IFAf = "J"THENRUN 

67 G0T064 

70 REM BERECHNUNG 

80 PRINTCHR$< 147)!PR INT!PR INT“ I "X " 11 " 

86 PRINT 
80 Y = 1 

100 FOR I=1T0N:Y=Y*X s PRINTI;YsNEXT 
110 PRINT!PRINT!PRINT"NEUER START <J/N> 

? " 

120 GET A*!IFA$=""THEN120 
130 IFAf = "J"THENRUN 

140 IFA$ = "N"THENPRINTCHR$< 154>CHR$< 147>" 

AUF WIEDERSEHEN":END 
150 GOTO 120 

Der Commodore 64 kann die N-te Potenz einer Zahl X auch mit PRINT X|N 
berechnen. 

10 REM 6POTENZEN 2 

20 POKE53280,6 s PRINTCHR$< 147>CHR$(5)" 

*** POTENZEN *** “ 
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30 PR INT:PR INT"DIESES PROGRAMM ERRECHNET 
DIE POTENZEN EINER EINGEGEBENEN ZAHL" 

40 PRINT:INPUT "ZAHL EINGEBEN ";X 
45 IFX=0THEN40 

50 PR INT: INPUT "WIEVIEL POTENZEN ";N 

60 IF N< 1ORN< >INT<N)THENPRINT:PR INT"NUR 
NATUERLICHE ZAHLEN BITTE":GOTO50 

61 IF< N+l > +LOGC ABS( X> X 1£6*L0G< 2)THEN70 
6£ PR INT:PRINT"RECHNERKAPAZITAET WIRD UE 
BERSCHRITTEN !" 

63 PR INT:PRINT"NEUSTART <J/N> ?" 

64 GETA$: IFA$ = ""THEN64 

65 IFAf="N"THEN70 

66 IFAf = "J"THENRUN 

67 G0T064 

70 REM BERECHNUNG 

80 PR INT:PR INT" I "X"TI" 

86 PRINT 

100 FÜR 1 = 1 TON:PR INTI;X11:NEXT 

110 PR INT:PR INT:PR INT"NEUER START <J/N> 

7 " 

1 £0 GET A$: IFA$ = ""THEN120 
130 IFA$ = ”J"THENRUN 

140 IFA$ = 11 N"THENPR INTCHR$< 154>CHR*< 147 > " 

AUF WIEDERSEHEN":END 
150 GOTO 120 

Wie Sie vielleicht wissen oder leicht nachprüfen können, braucht N keine ganze 
Zahl zu sein. 

Aber was bedeutet 2|1.7? Es sollte natürlich eine Zahl zwischen 2f1 und 2j2 sein, 
also zwischen 2 und 4. Tatsächlich hat 2j'1.7 ungefähr den Wert 3.24900959. 
Näherungsweise kann man den Wert von 2|1.7 (ohne Computer) mit Hilfe von 
Millimeterpapier bestimmen. Tragen Sie die erste bis fünfte Potenz von 2 ein - 
siehe Abbildung 15. Zeichnen Sie eine glatte Kurve durch diese Punkte - siehe 
Abbildung 16. Der ungefähre Wert von 2f1.7 läßt sich aus dieser Darstellung 
ablesen. 


Ähnlich könnte man die Potenzen anderer Zahlen graphisch darstellen. Das ist 
natürlich unnötig, da Ihr Commodore 64 Ihnen die Antwort sofort gibt. 

Negative Potenzen von Zahlen sind ebenfalls sinnvoll; sie sind einfach durch 
folgende Regel definiert: 
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X N = 1/X N . 


Also ist 2 2 = 1/2 2 , das heißt 1/4 oder 0.25. Die nullte Potenz einer Zahl ist 
konventionsgemäß gleich 1. 

Potenzen von Zahlen gehorchen der folgenden einfachen Multiplikationsregel: 

__ ^M + N 

So ist zum Beispiel 
3 4 * 3 2 = 3 6 
und 

10 2 * 10 3 = 10 '. 


QUADRATWURZELN 

Für eine Zahl X hat die Zahl X" 2 einen besonderen Namen: sie heißt Quadratwurzel 
von X. Die Quadratwurzel einer Zahl ist diejenige Zahl, deren Quadrat die ursprüng¬ 
liche Zahl ergibt. So ist die Quadratwurzel von 9 gleich 3 und die Quadratwurzel von 
2 ungefähr 1.41421356, wie Sie leicht nachprüfen können, indem Sie diese Zahl mit 
sich selbst multiplizieren. 
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40 


Die Schreibweise X 1 ' 2 für die Quadratwurzel von X paßt zur Multiplikationsregel für 
Potenzen: 

X 1 ' 2 * X’' 2 = X 1 

Das Ergebnis ist natürlich X. 

Die Quadratwurzel einer Zahl X kann man auch mit dem Commodore 64 bestim¬ 
men, indem man PRINT SQR(X) eingibt. 


IMAGINÄRE ZAHLEN 

Das Quadrat einer Zahl ist immer positiv, und deshalb sollte man nicht erwarten, die 
Quadratwurzel einer negativen Zahl finden zu können. Tatsächlich wird Ihnen der 
Commodore 64 auf die Frage PRINT SQR(-I) mit ILLEGAL QUANTITY ERROR 
antworten. 

Vielleicht haben Sie aber schon von imaginären und komplexen Zahlen gehört, die 
den Quadratwurzeln negativer Zahlen entsprechen. Mathematiker lassen sich nie¬ 
mals von scheinbar Unmöglichem wie der Quadratwurzel von -1 abschrecken. 
Man schafft einfach ein neues Symbol für diese Zahl. So schreibt man I für die 
Quadratwurzel von -1. Daran ist nichts Seltsames; I existiert ebenso wie negative 
Zahlen. 

Man kann I zu sich selbst und anderen Zahlen addieren und mit anderen Zahlen 
multiplizieren. Auf diese Weise entstehen Zahlen wie 

2 + 3*1, 1.41412*1, 9 — I, 10 — 8*1. 

Zahlen, die I enthalten, heißen komplexe Zahlen, im Unterschied zu gewöhnlichen 
oder reellen Zahlen. Jede komplexe Zahl läßt sich in der Form 

X + Y*l 

schreiben, wobei X und Y reelle Zahlen sind. X heißt Realteil und Y Imaginärteil der 
komplexen Zahl. 

Sobald wir einmal das Symbol I für SQR(—1) gewählt haben, können wir die 
Quadratwurzeln anderer negativer Zahlen bestimmen. 

SQR(X) = SQR(ABS(X)*SGN(X)) 

= SQR(ABS(X)) * SQR(SGN(X)) 

Zum Beispiel gilt 
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SQR<-9) = SQR(9) * SQR(-I) 

= 3*1. 

Natürlich arbeitet der Commodore 64 nicht in dieser Weise, denn er kennt keine 
komplexen Zahlen. 

Das folgende Programm zeigt, wie man komplexe Zahlen auf dem Commodore 64 
darstellen kann. So lassen sich einfache Rechnungen durchführen. Beachten Sie, 
daß komplexe Zahlen nach dieser Regel multipliziert werden: 

(A + B*l) + (C + D*l) = A*C + A*D*I + BTC + BTD*I 
= A*C + A*D*I + B*C*I + B*D*I*I 
= A*C + A*D*I + B*C*I + B*D*-1 
= A*C - B*D + (A*D + B*C)*I 

10 REM 7KOMPLEXE ZAHLEN 

S0 POKE53280,6sPRINTCHR$(147)" *** 

KOMPLEXE ZAHLEN ***" 

30 PR INT: PR INT''D IESES PROGRAMM MULT1PLIZ 
IERT UND DIVI - DIERT KOMPLEXE ZAHLEN. 

40 PR INTSPRINT"KOMPLEXE ZAHLEN WIE GEFRA 
GT EINGEBEN" 

60 PR INT:PR INT"ERSTE ZAHL" 

70 PRINT:INPUT-REALTEIL ";A 

80 INPUT"IMAGINAERTEIL "JB 
30 PRINTS PRINT"ZLEITE ZAHL“ 

100 PRINT:INPUT"REALTEIL ";C 

110 INPUT"IMAGINAERTEIL ";D 

180 PR INTSPRINT"WAS SOLL GEMACHT WERDEN 

7 " 

130 PR INT:PR INT"1> SOLLEN DIE ZAHLEN MUL 
TIPLIZIRT WERDEN 7" 

140 PR INT”2) SOLL DIE ERSTE VON DER ZWEI 
TEN ZAHL DIVIDIERT WERDEN 7" 

150 PR INT: INPUT"BITTE 1 ODER 2 EINGEBEN" 

.SN 

160 IFN=1ORN=2THEN170 
165 GOTO 150 

170 ON N GOSUB 300,350 

180 PR INT:PR INT s PR INT"EIN NEUER START <J 

/N > 7 " 
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190 GETA$:IFAf=""THEN190 
200 IF A$ = " J "THENRUN 

£10 IF A$="N"THENPRINTCHR$< 147>"AUF WIED 
ERSEHEN"!END 
220 GOTO 190 

300 REM MULTIPLIKATION 

310 PRINTCHR$( 147 > s PR INT!PR INT:PR INT“PRO 
DUKT DER 2AHLEN :" 

320 PRINT!PRINT"REALTEIL :"JA*C-B*D 

330 PRINT"IMAGINAERTEIL !“;A*D+B*C 

340 RETURN 

350 REM DIVISION 

360 H =C *C +D*D ! IFX=0THENPR INTCHR* < 147VDI 
VISION DURCH NULL NICHT MOEGLICH"RETURN 
370 PRINTCHR*< 147>:PR INT:PR INT:PR INT“QUO 
TIENT DER ZAHLEN :" 

380 PR INT:PR INT“REALTEIL :";<A*C+B* 

D ) /X 

390 PRINT“IMAGINAERTEIL :";<-A*D+B*C 
400 RETURN 


QUADRATISCHE GLEICHUNGEN 

Quadratische Gleichungen ergeben sich bei der Lösung ganz unterschiedlicher 
Probleme. Die allgemeine Form einer quadratischen Gleichung ist 

A‘X 2 + B*X + C = 0 

wobei A, B und C bekannte Zahlen sind und A von Null verschieden ist. Das 
Problem besteht darin, diejenigen Zahlen X zu finden, die die Gleichung erfüllen. 
Sie heißen die Wurzeln der quadratischen Gleichung. Zum Beispiel erfüllen die 
Werte X=1 und X=2 die folgende quadratische Gleichung: 

X 2 - 3*X + 2 = 0 

wie Sie leicht überprüfen können. 

Gewöhnlich hat eine quadratische Gleichung zwei Wurzeln, und häufig sind die 
Wurzeln komplexe Zahlen. 
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Es gibt eine ganz einfache Formel für die Wurzeln einer quadratischen Gleichung, 
und zwar sind die Wurzeln gegeben durch 

(-B + SQR(B 2 -4*A*C)) 


2*A 


und 


(-B-SQR(B 2 -4*A*C)) 

2*A 

Die Art der Wurzeln ergibt sich aus dem Teil der Formel, der die Quadratwurzel 
enthält, 

SQR(B 2 - 4*A*C) 
die sogenannte Diskriminante. 

Wenn B 2 - 4*A*C > 0 ist, gibt es zwei reelle Wurzeln. 

Wenn B 2 - 4*A*C = 0 ist, sind die beiden Wurzeln reell und stimmen überein. 

Wenn B 2 - 4*A*C < 0 ist, gibt es zwei komplexe Wurzeln. 

Das nächste Programm berechnet die Wurzeln einer quadratischen Gleichung. 


10 REM 8QUADRATISCHE GLEICHUNGEN 
20 POKE53280,6sPRINTCHR$< 147 V *** QU 
BORAT ISCHE GLEICHUNGEN ***" 

30 PR INT!PRINT"DIESES PROGRAMM LOEST QUA 
DRATISCHE GLEICHUNGEN 

40 PR INT:PRINT"GLEICHUNGSSCHEMA * A*X*X 
+ B*X + C =0 " 

50 PR INT: INPUT"WERT VON A“;A 

60 ipA =0THENPRINT”A KANN NICHT NULL SEIN 

";GOTQ50 

70 PR INT! INPUT 11 WERT VON B “ ; B 
80 PR INT! INPUT "WERT VON CDC 
90 REM BEGINN DER RECHNUNG 
100 D=B*B-4*A*C 
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110 ON SGN< D > +£ GOSUB 800,300,400 

120 PR INTs PR INT:PR INT"EIN NEUER START <J 

/N) ? " 

130 GETA$:IFA*=""THEN130 
140 IFA* = "J"THENRUN 

150 IFA* = "N"THENPRINTCHR*< 147 > "AUF WIEDE 
RSEHEN":END 

£00 REM ZWEI KOMPLEXE WURZELN 
210 PR INT:PR INTsPR INT s PR INT"ES GIBT ZWEI 
KOMPLEXE WURZELN" 

220 X = -B/A/2:Y=ABS< SQRC -D)/A/2) 

230 PR I NT s PR INTX ; " + " ; Y; *• * I '■ 

£40 PR INT: PR INTX; " - " ; y; ■■ * i » 

£50 RETURN 

300 REM ZWEI GLEICHE WURZELN 
3 10 PR INTSPRINTSPRINT:PR INT"ES GIBT ZWEI 
GLEICHE WURZELN" 

320 PRINTS PRINT-B/A/2 
330 RETURN 

400 REM ZWEI REALE WURZELN 
410 PRINT:PR INT:PRINTS PRINT"ES GIBT ZWEI 
REALE WURZELN" 

4£0 D=SQR< D > : X = C -B-D >/'£ : IFB< 0THENX=( -B+D 
)/2 

430 PR INT:PR INTX/A:PR INTC/X 
440 RETURN 


Falls es zwei reelle Wurzeln gibt, berechnet das Programm mit Hilfe der Formel für 
die Wurzeln zunächst eine Wurzel, und zwar die mit dem größten Betrag. Die 
andere Wurzel ergibt sich dann daraus, daß das Produkt der beiden Wurzeln gleich 
C/A ist. 


LÖSEN ANDERER GLEICHUNGEN 

Quadratische Gleichungen sind recht einfach zu lösen. Das gilt nicht für andere 
Gleichungen wie 

9*X 5 - 3*X 4 + X 3 - X 2 + 5*X - 4 = 0 
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Wenn die Gleichung nur nichtnegative ganze Potenzen enthält wie in diesem 
Beispiel, wird sie polynomische Gleichung genannt. Die größte von Null verschie¬ 
dene Potenz darin heißt der Grad des Polynoms. Die obige Gleichung hat den Grad 
5. Eine quadratische Gleichung ist eine polynomische Gleichung vom Grad 2. 
Quadratische Gleichungen haben gewöhnlich 2 Wurzeln; polynomische Gleichun¬ 
gen vom Grad N haben gewöhnlich N Wurzeln. 

Mit Ausnahme einiger Sonderfälle gibt es keine allgemeinen Formeln zur Lösung 
polynomischer Gleichungen. Das hat für die Mathematik sehr interessante Konse¬ 
quenzen gehabt, aber das ist ein anderes Thema, auf das wir hier nicht eingehen 
werden. 

Wir können jedoch unseren Commodore 64 dazu benutzen, durch wiederholtes 
Raten die Wurzeln einer polynomischen Gleichung zu bestimmen. Im wesentlichen 
prüft der Computer viele verschiedene Zahlen daraufhin, ob sie die Gleichung 
erfüllen. 

Das nächste Programm liefert ein Verfahren zur Bestimmung einer reellen Wurzel 
einer polynomischen Gleichung. Es ist mathematisch grob und findet gelegentlich 
auch dann keine Wurzel, wenn es eine gibt. Die zugrundeliegende Methode läßt 
sich jedoch an diesem Programm erkennen: 


10 REM 9P0LYN0ME 

2Q POKE53280,6sPRINTCHR$< 147? " *** POL 

YNOMISCHE GLEICHUNGEN ***" 

30 PR INT!PRINT"DIESES PROGRAMM VERSUCHT 
LOESUNGEN ZU FINDEN FUER GLEICHUNG 
35 PR INT"GEN WIE s" 

40 PRINT!PR INT“A*XtN + B*Xt<N-l> + ... + 

C*X + D = 0" 

50 PRINT 

60 INPUT "GRAD DES POLYNOMS UN 

70 IF N<2 OR N< >INT< N> THENPRINT"NATUERL 

ICHE ZAHLEN EINGEBEN"sGOTO60 

80 D IMAC N> s PR INT - . PR INT"KOEFF IZENTEN TERM 

NACH TERM EINGEBEN" 

85 PRINT 
90 FOR 1=0 TO N 

100 IF I<N THENPRINT" KOEFFIZENT FUER Xt 
" J MID* C STR* < N-I>,2 > ; 

110 IF I=N THENPRINT" KONSTANTE < D)"f 
120 INPUTAC I> 

130 IFAC 0 > =0THENPR INT 11 KE INE NULL BITTE"! 
GOTO 100 
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140 NEXT 

150 PR INT : PR INT"GEBEN SIE DIE GRENZEN EI 
N, IN DENEN DIE WURZELN GESUCHT WERDEN"; 

160 PRINT" SOLLEN" 

170 PRINT:INPUT"UNTERER WERT";A 
180 INPUT"OBERER WE R T";B 

130 IFA >=BTHENPR INT: PR INT“ERSTER WERT MU 
SS KLEINER SEIN":GOTO 170 
200 REM SUCHE 

210 S=B-A:T=0:TEST =-1:D = 1E-9 

220 PRINTCHR$< 147):PR INT:PR INT"TEST LAUF 

";T +1; 10tT; "DIVISIONEN **":GOSUB300 

230 IF TEST THEN S=S/10:T=T+1:IF T<4 THE 
N 220 

240 IFT=4THENPRI NT:PR INT"KANN LEIDER KE I 
NE WURZEL FINDEN" 

250 PR INT:PR INT:PR INT"NEUER START < J/N) 

7 " 

260 GETG$: IFG$ = " "THEN260 
270 IFG$ = "J"THENRUN 

280 IFG$ = "N"THENPRINTCHRf( 147 >"AUF WIEDE 
RSEHEN":END 
290 GOTO260 

300 REM SCHRITT BEI SCHRITT SUCHE 

310 X=A:GOSUB400:X1=A:Y1=Y 

320 FORX = <A+S>TOBSTEPS 

330 GOSUB400:X2=X:Y2=Y 

340 IFY1*Y2< =DTHENGOSUB500 

350 Y1=Y2:X1=X 

360 NEXT 

370 RETURN 

400 REM BERECHNUNG DES POLYNOMS 
410 Y=A<0>:FOR1 = 1 TON:Y=Y*X+A( I):NEXT 
420 RETURN 

500 REM SCHLUSS BERECHNUNG 

510 PR INT:PR INT"SCHLUSS BERECHNUNG":B$ = " 

II 

520 IFHBS<Y1)<DTHENX= X1:GOTO590 
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530 IF PBS< Y2)<DTHENK =X1 :GOTO590 

540 2 = X : X =< X 1 +X2 ) /2 :GOSUB400 

550 IF i4BS< Y)<DTHENB* = "WAHRSCHEINLICH” :G 

OTO530 

570 IFY+Y2>0THENX2=X:GOTO540 
580 X 1 =XiGOTO540 

530 PR INT: PR INTB$ 11 WURZEL BE I " ; X : TEST=0 : 

X =B 

600 RETURN 

Manchmal sagt das Programm von einer Zahl, daß sie vielleicht eine Wurzel ist. Um 
zu prüfen, ob sie tatsächlich eine ist, oder wie nahe sie einer Wurzel ist, beantwor¬ 
ten Sie zunächst die Frage nach einem weiteren Programmdurchlauf mit N, und 
geben Sie dann folgendes ein; 

X=Z:GOSUB 400:PRINT Y 

Der Wert der daraufhin ausgegebenen Zahl sagt Ihnen, wie nahe Sie an einer 
Wurzel sind. 


DAS NEWTONSCHE VERFAHREN 

Das im vorigen Abschnitt angegebene Verfahren zur Bestimmung von Wurzeln 
einer Gleichung läßt sich mit dem sogenannten Newtonschen Verfahren verbes¬ 
sern. 

Angenommen, wir wollen Wurzeln der polynomischen Gleichung 
9*X 5 - 3*X 4 + X 3 - X 2 + 5*X - 4 = 0 

finden. Wir bezeichnen das Polynom mit P(X) und schreiben P’(X) für das folgende 
Polynom: 


5*9*X 4 - 4*3*X 3 + 3*X 2 - 2*X + 5. 


P’(X) ist die Ableitung von P(X), aber das braucht uns nicht zu kümmern. Wenn Y 
näherungsweise eine Wurzel der Gleichung P(X) = 0 ist, so ist 

Y - P(Y)/P'(Y) 
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gewöhnlich eine bessere Näherung, falls P’(Y) verschieden von Null ist. 

Das nächste Programm verwendet diese Technik zur Bestimmung von Wurzeln 
polynomischer Gleichungen. 

10 REM 10POLYNOME £ 

£0 POKE53380,6:PRINTCHR$( 147)" *** POL 

YNOMISCHE GLEICHUNGEN * * *" 

£5 PRINT" *** NEWTON'S METHODE *** 

11 

30 PR INTSPRINT"DIESES PROGRAMM FINDET LO 
ESUNGEN NACH NEWTON'S METHODE FUER "; 

35 PRINT"GLEICHUNGEN WIE 

40 PR I NT s PR INT" A t-X tN + B*Xt<N-l> + ... + 

C*X + D = 0“ 

50 PRINT 

60 INPUT "GRAD DES POLYNOMS ";N 
70 IF N< £ OR f-K > INT< N> THENPR INT"NATUERL 
ICHE ZAHLEN EINGEBEN":GOTO60 
80 DIMA< N):PR INT:PR INT"KOEFFIZENTEN TERM 
NACH TERM EINGEBEN" 

85 PRINT 
80 FOR 1=0 TO N 

100 IF I< N THENPRINT" KOEFFIZENT FUER Xt 
MID$( STR$< N-I),£?; 

110 IF I=N THENPRINT" KONSTANTE <D>"; 

130 INPUTA< I ) 

130 IFAC 0 > =0THENPRINT"KE INE NULL BITTE": 

GOTO 100 
140 NEXT 

150 PR INTSPRINT"GEBEN SIE EINEN ANGENOMM 
EN WERT FUER DIELAJRZELN EIN"; 

160 PR INT s INPUT"WERT" ; X 
170 REM BERECHNUNG P'<X> 

180 0 IM B<N):FORI=0TONSB< I>=CN-l)*A( I)SN 
EXT 

190 GOSUB500 s REM WURZELN FINDEN 
200 REM ENDE 

210 PR INTSPRINT SPRINT"NEUER START <J/N> 

? " 

220 GETG$:IFG$=""THEN3S0 
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230 IFG$ = "J"THENRUN 

240 IFG$ = "N"THENPR INTCHR$< 147 ) “AUF WIEDE 
RSEHEN"!END 
250 GOTO210 

400 REM LOESUNG POLYNOM 

410 Y=A< 0>:FORI=1TON:Y=Y*X+A< IXNEXT 

420 Y 1 =B( 0 ) s FOR I = 1 TO< N- 1 ) : Y1 = Y1 *X +B< I > ’• N 

EXT 

430 RETURN 

500 REM LOESUNGEN FINDEN 
505 PR INTCHR$< 147) 

510 J = 1 :D= 1E-9:GOSUB400 
520 IFRBSC YXDTHEN590 

530 IFY1=0THENPRINT"DIV IS ION DURCH NULL 

- NEUER VERSUCH":RETURN 

540 2 =X:X =X-Y/Y1iJ=J+1:GOSUB400 

550 IFHBS't YXDTHEN590 

560 IFABS< Z -XX DTHENBf = " WAHRSCHEINLICH " 
:GOTO590 

570 IFJ >1000THENPRINT"LEIDER KEINE WURZE 
L GEFUNDEN" 

580 GOTO530 

590 PR INT:PRINTB*"WURZEL BE I " JX 
600 RETURN 


DIE EXPONENTIALFUNKTION 

Funktionen wie 2 X oder 10 x heißen Potenzfunktionen, weil die Variable X als Potenz 
eingeht. Potenzfunktionen erhält man auf dem Commodore 64 durch PRINT 2|X 
etc. Eine wichtige Potenzfunktion ist im Commodore 64 eingebaut: die Exponential¬ 
funktion EXP(X). 

Die Exponentialfunktion beruht auf Potenzen der Zahl E, die ungefähr den Wert 
2.71828183 hat. Verwechseln Sie dieses E nicht mit dem E, das in der Gleitkomma¬ 
darstellung von Zahlen vorkommt. Die Zahl E ist definiert durch 

e = i + -L + — + -L + — + — + — + ... 

1! 2! 3! 4! 5! 6! 
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Dabei bedeutet. . daß die Summe endlos fortgesetzt ist, und das Symbol ! steht 
für Fakultät, definiert durch 

N! = N * (N -1) * (N-2) * ... * 2 * 1 

also das Produkt der ganzen Zahlen von 1 bis N. Z.B. ist 4! gleich 4*3*2*1, also 24. 
Die Exponentialfunktion ist definiert durch 

EXP(X) = E x 

= E jX. 

Insbesondere gilt also E = EXP(1). Geben Sie Ihrem 64 versuchsweise folgendes 
ein: 

E = EXP(1) 

PRINT Ef5, EXP(5) 

Aus den Eigenschaften von Potenzfunktionen ergeben sich die Eigenschaften der 
Exponentialfunktion: 

EXP(X)*EXP(Y) = EXP(X+Y) 

EXP(X)/EXP(Y) = EXP(X-Y) 

Mit Hilfe einer ganz einfachen Formel läßt sich EXP(X) für jedes X berechnen. Diese 
Formel sieht so aus: 

exp(X) = i + JL + .xi + .XL + -Xl + -Xi + .xi + 

1! 2! 3! 4! 5! 6! 


DIE LOGARITHMUSFUNKTION 

Welche Zahl X erfüllt EXP(X) = 3? Da EXP(1) = 2.71828183 gilt, muß X etwas 
größer als 1 sein. Tatsächlich lautet die Antwort: ungefähr 1.09861229. 

Die Zahl X, die EXP(X) = N erfüllt, heißt (natürlicher) Logarithmus von N. Sie wird 
gewöhnlich mit LN(N) bezeichnet. Wie die meisten Mikrocomputer schreibt der 
Commodore sie jedoch als LOG(N). 

Die Logarithmusfunktion hat folgende Eigenschaften: 

LOG(X*Y) = LOG(X) + LOG(Y) 

LOG(X/Y) = LOG(X) - LOG(Y) 

LOG(X N ) = N*LOG(X) 
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Wegen dieser Eigenschaften sind Logarithmen wichtig für Multiplikationen, Divisio¬ 
nen etc. ohne Computer. 

Zwischen der Exponentialfunktion und der Logarithmusfunktion bestehen die Be¬ 
ziehungen 

EXP(LOG(N)) = N 
LOG(EXP(N)) = N 

Eine einfache Anwendung der Logarithmusfunktion ist das Testen großer Zahlen. 
Z.B. erfüllt die Zahl X die Beziehung 

X< 10|N 

dann und nur dann, wenn die Beziehung gilt: 

LOG(X)<N*LOG(10). 

Eine solche Umformulierung ist nützlich, weil die Zahl 10|N selbst einen OVER¬ 
FLOW ERROR hervorrufen kann. 


WURZELN ANDERER FUNKTIONEN 

Ist eine Funktion wie zum Beispiel X*EXP(X) + 1 gegeben, so versteht man unter 
einer Wurzel der Funktion eine Zahl, für die die Funktion den Wert 0 ergibt. 

Zwei Programme zur Bestimmung der Wurzeln einer quadratischen Gleichung 
wurden weiter oben vorgesteilt. Das erste davon läßt sich so verändern, daß es 
Wurzeln anderer Funktionen bestimmt. Wenn Sie das Programm starten, werden 
Sie zur Eingabe von zwei Zeilen aufgefordert: 

100 DEF FNA(X) = Funktion von X eingeben 
GOTO 100 

Geben Sie in der ersten Zeile die Funktion ein, von der Sie die Wurzel bestimmen 
wollen. 

10 REM 11WURZELN 

20 POKE53280 ,6 :PRINTCHR$C 147>" **+ WURZ 

ELN ANDERER FUNKTIONEN ***“ 

30 PR INT s PR INT"DIESE PROGRAMM VERSUCHT W 
URZELN VON FUNKTIONEN ZU FINDEN. 
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40 PR INT"BITTE FUNKTION IN DER FOLGENDEN 
FORM EINGEBEN:" 

50 PR INT SPRINT "a|l 00 DEF FNA< X > = X*EXP<X 

> 

60 PR INT SPRINT" LPANACH" 

70 PR INT: PR INT" sf GOTO 100a" : END 
100 DEF FNA( X ) =X*EXP< X ) *LOG< ABS< X0 ) + 1 > 
150 PRINT CHR$< 147):PR INT:PR INT"GEBEN SI 
E DIE GRENZEN EIN, IN DENEN EINE"; 

160 PR I NT 11 KURZ EL GESUCHT WERDEN SOLL !” 

170 PRINT:INPUT"UNTERER WERT ";A 

180 INPUT”OBERER WERT " ■ B 

190 IF A>B THEN PR INT:PR INT"ERSTER WERT 

MUSS KLEINER SEIN" 

£00 REM SUCHE 

210 S =B-A:T=0:TEST =-i:D = lE-9 
2£0 PRINTCHR$( 147):PR INT:PRINT"TEST LAUF 
"■T+l;" ** 10tT; " DIVISIONEN * * 11 : GOSU 

B300 

£30 IF TEST THEN S=S/10sT=T+1:IF T<4 THE 
N GOTO220 

240 IF T = 4 THEN PR INT:PRINT"LEIDER KEINE 
WURZEL GEFUNDEN" 

250 PR INT:PR INT:PR INT"NEUER START <J/N) 

? " 

260 GET G$: IF G$ ="" THEN 260 
270 IF G$="J" THEN RUN 

280 IF G$ = "N" THEN PRINT CHR$< 147VAUF W 
IEDERSEHEN":EKO 

300 REM SCHRITT FUER SCHRITT SUCHE 

310 X =A:Y=FNA( X):X1=A:Y1=Y 

320 FOR X =A + S TO B STEP S 

330 Y=FNA<X >:X£ =X s Y2=Y 

340 IF Y1 *Y2< =D THEN GOSUB500 

350 Y1=Y2:X1=X 

360 NEXT 

370 RETURN 

500 REM SCHLUSSBERECHNUNG 

510 PRINTsPRINT"SCHLUSSBERECHNUNG” 

520 IF ABSCY1XD THEN X=X1:GOTO590 
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530 IF RBS< YH XD THEN X=X2sGOTO590 
540 Z =Xs X=< X1+X2 >/2s Y=FNA<X) 

550 IF ABS( YXD THEN590 

560 IF RBS< 2-X ) < D THEN B$="WRHRSCHEINLIC 
H "s GOTO590 

570 IF Yt-Y2 >0 THEN X2 =X : GOTO540 
580 X1=Xs GOTO540 

590 PR INT:PRINTB$"WURZEL BEI ";X:TEST=0: 
X =B 

600 RETURN 


Wenn Sie ein Ergebnis mit dem Zusatz PROBABLY (vielleicht) erhalten, stoppen 
Sie das Programm und geben 

PRINT FNA(Z) 

ein, um zu ermitteln, wie nahe Sie an einer Wurzel sind. 
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KAPITEL 5 
FOLGEN 


Folgen (und Reihen) sind wichtige Begriffe, die überall Vorkommen. Eine Folge ist 
einfach eine Liste von Zahlen, zum Beispiele 

1, 2, 3, 4, 5, 19, 7, 8, 12 
20, 18, 16, 14, 12, 10, 8 
1 , 0 . 1 , 0 . 01 , 0 . 001 , 0 . 0001 . 0 . 00001 . 0.000001 

Die einzelnen Zahlen oder Elemente der Folge heißen Glieder der Folge. Gewöhn¬ 
lich ist eine Folge durch eine bestimmte Regel erzeugt, wie etwa die letzten beiden 
oben. Die zweite ist aus der Formel 22 - 2*N für N = 1 bis 7 entstanden, und die 
letzte aus 10/10' für N = 1 bis 7. Ihr Commodore 64 kann leicht Folgen erzeugen, 
wie das nächste Programm zeigt. Setzen Sie Ihre eigene Formel (von N abhängig) 
in der zweiten Zeile ein. 

10 REM SEQUENCE GENERATOR 

20 DEF FNA(N) = (von N abhängige Formel einsetzen) 

30 FOR N = 1 TO 10 

40 PRINT FNA(N);:IF N < 10 THEN PRINT CHR$(157);“,”; 

50 NEXT 
60 PRINT 

Hier sind einige mit diesem Programm erzeugte Folgen. Können Sie die jeweils 
zugrundeliegende Formel erkennen? Überprüfen Sie Ihre Vermutung, indem Sie 
die Formel in das obige Programm einsetzen. (Die Antworten werden weiter unten 
in diesem Abschnitt gegeben; die nachfolgenden Abschnitte enthalten weitere 
Beispiele.) 

(a) 1, 6, 11, 16, 21, 26, 31, 36, 41, 46 

(b) 1, 2, 4, 8, 16, 32, 64, 128, 256, 512 

(c) 1, 4, 9, 16, 25, 36, 49, 64, 81, 100 

(d) 1, 2, 4, 7, 11, 16, 22, 29, 37, 46 

(e) 4, 4, 8, 12, 20, 32, 52, 84, 136, 220 

Für die ersten Folgen sind die Formeln nicht schwer zu bestimmen. Für (a) ist es 
5*N-4, für (b) ist es 2|N/2, und für (c) ist es N*N. Die vierte (d) ist nicht so leicht zu 
erraten, es ist (N*N - N + 2)12. Auf die Formel für die fünfte Folge (e) schließlich 
kann man unmöglich kommen, wenn man sie noch nie gesehen hat; sie lautet: 
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1+(-1)|N 

N*(-1)|N 

INT(SIN(N)*10) 


ARITHMETISCHE FOLGEN 

Eine arithmetische Folge oder arithmetische Progression ist eine Folge, in der jedes 
Glied die Summe des vorhergehenden Gliedes und einer Konstanten ist. Die 
allgemeine Formel für eine arithmetische Folge ist 

A + (N - 1)*D 

wobei A das erste Glied der Folge und D die gemeinsame Differenz ist. Hier sind 
noch mehr Beispiele für arithmetische Folgen: 

5,10,15, 20, 25, 30, 35, 40, 45, 50 (A = 5, D = 5) 

1,1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5 (A = 1 , D = 0.5) 

0, 2, 4, 6, 8, 10, 12,14,16, 18 (A = 0, D = 2) 

Das unten abgedruckte Programm kann Ihnen helfen, arithmetische Folgen zu 
untersuchen. Sie geben das erste Glied der Folge, die gemeinsame Differenz und 
die gewünschte Anzahl der Glieder ein. Beachten Sie, daß keine Formel erforderlich 
ist, weil der Commodore 64 die Folge iterativ berechnet. Zusätzlich berechnet das 
Programm die Summe der Folgeglieder und druckt das Ergebnis aus. 

10 REM 1SFOLGEN ARITHMETRISCHE 
20 POKE53280,6:PRINTCHR$< 147 >" *** ft 

RITHMETISCHE FOLGEN ***" 

30 PR INT!PRINT"DIESES PROGRAMM BERECHNET 
ARITHMETISCHE FOLGEN" 

40 PRINT: INPUT"ERSTE GLIED DER FOLGE "fA 
50 PR INT> INPUT"GEMEINSAME DIFFERENZ ”*D 
60 PR INT: INPUT"ANZAHL DER GLIEDER "JN 
70 IF N< 1 OR N< >INT( N> THEN 60 
80 REM DIE FOLGE 

90 PR INTCHR$ < 147 > s PR INT s PR INT"DIE FOLGE 

IN DEN GRENZEN 1 BIS"?N 

100 TERM=ft:SUM=0 

105 PRINT 

1 10 FOR 1 = 1 TO N 
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1S0 IF38-POS<0XLEN< STR*< TERM) > THENPRIN 
T 

130 PR INTCHR$ ( 26 >TERM 

140 SUM=SUM+TERM:TERM=TERM+D 

150 NEXT I 

160 PR INTs PR INT!PRINT"DIE SUMME = "SUM 
170 PR INTSPRINTSPRINT:PR INT"NEUER START 

< J/'N> " 

160 GETG$ s IF6$ = ""THEN180 
130 I FG$ = "J"THENRUN 

200 IFG$ = "N"THENPRINTCHR$C 147 >"AUF WIEDE 
RSEHEN"SEND 
210 GOTO 180 


WAS MÖCHTEN SIE LIEBER? 

Angenommen, der Verleger dieses Buches bietet Ihnen eine Steile als freier 
Mitarbeiter an und fragt Sie dann, wie Sie am liebsten bezahlt werden möchten: 
’Was möchten Sie lieber? Mit DM 3000 für sechs Monate anfangen und dann alle 
sechs Monate DM 120 mehr bekommen, oder mit DM 6120 für ein Jahr anfangen 
und dann jedes Jahr DM 240 mehr bekommen?' Die in einem Jahr erhaltenen 
Geldbeträge bilden in beiden Fällen eine arithmetische Folge. Aber ein Angebot ist 
viel besser als das andere: das erste. Erkennen Sie warum? Sehen Sie sich die 
folgenden Berechnungen an: 


Erstes Angebot 

Erstes Jahr (Die ersten 6 Monate DM 3000) 
(Die zweiten 6 Monate DM 3120) 
DM 6120 

Zweites Jahr (Die ersten 6 Monate DM 3240) 
(Die zweiten 6 Monate DM 3360) 
DM 6600 

Drittes Jahr (Die ersten 6 Monate DM 3480) 
(Die zweiten 6 Monate DM 3600) 
DM 7080 


Zweites Angebot 


DM 6120 

DM 6360 

DM 6600 


Sie sehen, daß das erste Angebot besser ist, wenn Sie die Steile länger als ein Jahr 
behalten. Beachten Sie, daß beim ersten Angebot das Gehalt für die ersten 6 
Monate jeden Jahres um DM 240 wächst, demnach die jährliche Erhöhung DM 480 
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beträgt. Das jährliche Gehalt beim ersten Angebot ist eine arithmetische Folge nach 
der Formel: 

6120 + (N —1)*480, 

beim zweiten dagegen nach der Formel: 

6120 + (N-1)*240 

Angenommen, Sie bekämen ein drittes Angebot: 'Sie fangen mit DM 1440 im 
Quartal an und bekommen alle drei Monate DM 60 mehr.’ Welches Angebot würden 
Sie jetzt vorziehen? Die Antwort ist Ihnen hoffentlich klar. Die Rechnung sähe so 
aus: 


Erstes Jahr 


Zweites Jahr 


Drittes Angebot 
(Die ersten 3 Monate DM 1440) 
(Die zweiten 3 Monate DM 1500) 
(Die dritten 3 Monate DM 1560) 
(Die vierten 3 Monate DM 1620) 
(Die ersten 3 Monate DM 1680) 
(Die zweiten 3 Monate DM 1740) 
(Die dritten 3 Monate DM 1800) 
(Die vierten 3 Monate DM 1860) 


DM 6120 


DM 7080 


Die Erhöhung für ein Quartal beträgt von Jahr zu Jahr DM 240, aber das bekommen 
Sie alle drei Monate. Also ist die jährliche Erhöhung DM 960. 


GEOMETRISCHE FOLGEN 

Eine andere verbreitete Art von Folgen ist die geometrische Folge oder geometri¬ 
sche Progression. In einer geometrischen Folge ist das Verhältnis zweier aufeinan¬ 
der folgender Glieder konstant, zum Beispiel 

2, 6, 18, 54, 162 

wo jedes Glied (außer dem ersten) das Dreifache des vorhergehenden ist. Die 
allgemeine Formel für eine geometrische Folge ist 
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A * Rf(N-l) 


wobei A das erste Glied und R das gemeinsame Verhältnis ist. 

Hier sind noch mehr Beispiele für geometrische Folgen: 

4, 2,1,0.5, 0.25, 0.125, 0.0625 (A = 4, R = 0.5) 

2,-4,8,-16,32,-64,128 (A = 2, R = -2) 

Das nächste Programm kann Ihnen helfen, geometrische Folgen zu untersuchen. 
Sie geben das erste Glied der Folge, das gemeinsame Verhältnis und die 
gewünschte Anzahl der Glieder ein. Beachten Sie, daß keine Formel erforderlich ist, 
weil der Commodore 64 die Folge iterativ berechnet. Zusätzlich addiert das Pro¬ 
gramm alle Glieder der Folge und druckt das Ergebnis aus. 

10 REM 13F0LGEN GEOMETRISCHE 
20 POKE53280,SsPRINTCHR$(147 >" *** G 

EOMETRISCHE FOLGEN ***" 

30 PRINTS PRINT"DIESES PROGRAMM BERECHNET 
GEOMETRISCHE FOLGEN" 

40 PRINTsINPUT"ERSTE GLIED DER FOLGE 
1 A 

50 PR INTs INPUT"GEMEINSAMES VERHAELTNIS " 

;r 

60 PRINTsINPUT"AN2AHL DER GLIEDER " 

;n 

70 IF N< 1 OR N< > INTC N> THEN 60 
80 REM DIE FOLGE 

30 F'R I NTCHR$< 147 ) s PR INT s PR I NT "DIE FOLGE 

IN DEN GRENZEN 1 BIS"?N 

100 TERM=A s SUM=0 

105 PRINT 

110 FOR 1=1 TO N 

120 IF38-P0S< 0 X LEN< STR$( TERM) > THENPRIN 
T 

130 PRINTCHR$< 26 >TERM 

140 SUM=SUM+TERMSTERM=TERM*R 

150 NEMT I 

160 PR INTs PR INT s PR INT"DIE SUMME = "SUM 
170 PRINTsPRINTsPR INT:PRINT"NEUER START 

< J / N > ” 

180 GETG$ : IFG$ = 11 "THEN180 
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190 I FG$ = "J"THENRUN 

200 IFGf = "N"THENPRINTCHR$< 147 > "AUF WIEDE 
RSEHEN":END 
210 GOTO 180 


ZINSEN 

Am 1. Januar zahlt eine Frau DM 100 in eine Bank ein, die jedes Jahr (am 
Jahresende) 6% Zinsen gewährt. Auf weichen Betrag sind die DM 100 nach 10 
Jahren auf der Bank angewachsen? 

Nach einem Jahr hat die 

100 + 100 * 6/100 
= 100 + 100 * 0.06 
= 100 * 1.06 

also DM 106. Nach zwei Jahren hat sie 

100 * 1.06 + (100 * 1.06) * 0.06 
= 100 * 1.06 * 1.06 

Wie Sie jetzt herausfinden können, hat sie nach 10 Jahren 
100 * 1.06X10. 

Der Gesamtbetrag nach jedem Jahr bildet eine geometrische Reihe wie folgt: 

100 * 1.06, 100 * 1.06|2, 100 * 1.06T3, 100 * 1.06|4, 

100 ’ 1.06|5, 100 * 1.06|6, 100 * 1.06|7, 100 * 1.06|8, 

100 * 1.06|9, 100 * 1,06|10 

Beginnt man allgemeiner mit einem Betrag A und erhält 1% Zinsen im Jahr, dann ist 
der ursprüngliche Betrag nach N Jahren auf folgende Summe angewachsen: 

A * (1 + l/100)fN 
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TÄGLICHE ZINSEN 

Werden DM 1000 bei einer Sparkasse angelegt, die am Ende jeden Jahres 6% 
Zinsen zahlt, so beträgt das Guthaben nach einem Jahr 

1000 * 1.06 


vorausgesetzt, es finden keine weiteren Einzahlungen oder Auszahlungen statt. 
Wenn jedoch die Bank alle 6 Monate Zinsen zahlt (und auch die schon gezahlten 
Zinsen verzinst, also Zinseszins zahlt), ist das Guthaben am Jahresende 

1000 * (1.03)|2. 


Allgemeiner ausgedrückt: Zahlt die Sparkasse 6% Zinseszins N-mal im Jahr, so 
sind die DM 1000 am Ende eines Jahres auf den folgenden Betrag angewachsen: 

1000 * (1 + 0.06/N)jN 

Nachstehende Tabelle zeigt die Beträge, die sich bei unterschiedlich häufiger 
Zinseszinszahlung ergeben: 


N 6% Zinseszins 

1 jährlich 
2 halbjährlich 
4 vierteljährlich 
6 zweimonatlich 
12 monatlich 
52 wöchentlich 
365 täglich 
8760 stündlich 


Guthaben am Ende des Jahres 
(auf Pfennig gerundet) 

DM 1060.60 
DM 1060.90 
DM 1061.36 
DM 1061.52 
DM 1061.68 
DM 1061.80 
DM 1061.83 
DM 1061.84 


Diese Tabelle hat der Commodore 64 mit dem folgenden einfachen Programm 
erstellt: 


10 REM 14ZINSESZINS 

20 POKE53280,6!PRINTCHR$< 147>" *** 

ZINSESZINS FORMEL ***'• 

25 PR INT s PR INT"WIE VERAENDERT SICH DER B 

ETRAG VON 1000 DM, WENN IM JAHR ”; 

26 PRINT"MEHRFACH ZINSEN GEZAHLT WERDE 
N":PR INT 
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30 PR INTSPRINT"BITTE ANZAHL DER ZINSPERI 
ODEN ANGEBEN ." 

40 PRINT: INPUT "2 INSPERIODEN ";Ns IFN<=0OR 
N< > I NTC N;>THEN40 

50 T=1000*< 1 +0.06/N5 tN: PR INT: PR INTINT< T* 
1005/100 

60 PR INT SPRINT"NEUER START <J/N5 ?" 

70 GETG*:IFG$= H "THEN70 
80 IFGt = "J"THENRUN 

90 IFG* = "N"THENPRINTCHR*< 147> "AUF WIEDER 
SEHEN":END 
100 GOTO70 

Um andere Zinssätze einsetzen zu können, sind einige Zusätze und Änderungen 
erforderlich. Veränderungen und Zusätze im nächsten Listing sind durch ein Stern¬ 
chen am Anfang einer Zeile gekennzeichnet. 

10 REM 15ZINSESZINS VAR IABELER ZINSSATZ 
30 POKE533S0,6sPRINTCHR$( 147 5" *** VARIA 
BLE ZINSESZINS FORMEL ***" 

30 PRINTs PR INT"WIE VERAENDERT SICH EIN B 
ETRAG VON 1000 DM, WENN IM JAHR "; 

40 PR INT"MEHRFACH ZINSEN GEZAHLT WERDE 
N":PRINT 

50 PR INT SPRINT"BITTE ZINSSATZ EINGEBEN" 

60 PRINT:INPUT"ZINSSATZ “S 1 

70 PR INT:PR INT"BITTE ANZAHL DER ZINSPERI 

ODEN ANGEBEN ." 

80 PR INT: INPUT"ZINSPERIODEN “ i N: IFN<=0OR 
N<>INT<N5THEN80 

90 T=1000*< 1 + I /100/N 5 tN:PR INT s PR INTINT< T 
* 100 5/100 

100 PR INT:PR INT:PRINT"NEUER START ( J/N5 

7 " 

110 GETG$: IFG$ = ""THEN110 
130 IFG$="J"THENRUN 

130 IFG$ = "N"THENPR INTCHR< : < 147 5 "AUF WIEDE 
RSEHEN"SEND 
140 GOTO 1 10 
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VERDOPPELN ODER AUFGEBEN 


Manche halten es für möglich, bei Glücksspielen niemals zu verlieren. 

Nehmen Sie zur Veranschaulichung dieses Glücksspiel: ‘Eine unverfälschte Münze 
wird geworfen, gleichzeitig machen Sie Ihren Einsatz. Zeigt die Münze Kopf, 
erhalten Sie das Doppelte Ihres Einsatzes zurück.’ 

Das folgende Argument soll begründen, daß Sie nie zu verlieren brauchen: Setzen 
Sie erst DM 1. Wenn Sie gewinnen, hören Sie auf. Wenn Sie verlieren, spielen Sie 
mit einem Einsatz von DM 2 weiter. Verdoppeln Sie Ihren Einsatz nach jedem 
verlorenen Spiel und spielen Sie weiter. Sobald Sie gewinnen, hören Sie auf und 
haben auch dann mehr als Ihren Einsatz zurück. 

Nehmen Sie an, Sie verlieren die ersten vier Spiele und gewinnen das fünfte. Es 
ergibt sich die Tabelle: 



Einsatz 

Verlust 

Gewinn 

Erster Wurf 

DM 1 

DM 1 


Zweiter Wurf 

DM 2 

DM 2 


Dritter Wurf 

DM 4 

DM 4 


Vierter Wurf 

DM 8 

DM 8 


Fünfter Wurf 

DM 16 


DM 16 

INSGESAMT 


DM 15 

DM 16 


NETTOGEWINN = DM 1 

Die untenstehende Folge ist eine geometrische Folge. Überzeugt Sie dieses 
Argument davon, daß man niemals zu verlieren braucht? 


FIBONACCI-FOLGEN 

Am Anfang dieses Kapitels hatten wir die Folge 

(e) 4, 4, 8, 12, 20, 32, 52, 84, 136, 220 

die durch die Formel erzeugt wird: 

4’INT(((0.5 + SQR(5)/2)|N - (0.5 - SQR(5)/2)|N) / SQR(5)). 

Es gibt jedoch ein einleuchtenderes Verfahren, diese Folge zu erzeugen. Jedes 
Glied außer den ersten beiden ist die Summe der zwei vorhergehenden Glieder. 
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4 + 4 = 8 
4 + 8 = 12 
8 + 12 = 20 

und so weiter. 

So erzeugte Folgen heißen Fibonacci-Folgen. Im Jahr 1202 beobachtete Leonardo 
von Pisa, genannt Fibonacci, eine solche Folge im Zusammenhang mit der Kanin¬ 
chenzucht. 

Hier sind zwei weitere Fibonacci-Folgen: 

2, 5, 7, 12, 19, 31,50 

3, 5, 8, 13, 21, 34, 55 

Das nächste Programm erzeugt Fibonacci-Folgen bis zum Überdruß: 

10 REM 16F I BONACC I FOLGEN 

20 POKE53280,6:PRINTCHR$< 147>" *** 

FIBONACCI FOLGEN ** + " 

30 PR INT3 PR INT"DIESES PROGRAMM BERECHNET 
FIBONACCI- FOLGEN" 

40 PR INT:PRINT"GEBEN SIE ZWEI GANZE ZAHL 
EN EIN. TRENNENS IE DIESE DURCH "; 

45 PR INT"EIN KOMMA." 

50 PRINT:INPUT"ZAHLEN";U,V 

S0 IF U< >INTCU>ORV< >INTCV>THENPRINT:PR IN 

T"GANZE ZAHLEN BITTE" 

70 PRINT:INPUT-WIEVIELE FOLGEGLIEDER”;N 
80 I FN< 1 ORN< > INT< N > THEN70 
90 REM DIE FOLGE 

100 PRINTCHR$< 147):PR INT:PR INT"FIBONACCI 
-FOLGE MIT GRUNDGLIEDERN"JU;",";V 
105 PRINT 
110 FOR 1 = 1 TON 

120 IF 3S-P0SC 0 X LEN( STR$< U))THENPRINT 
130 PRINTCHR$< 29);U 
140 W=U+V:U=V:V=W 
150 NE.7TI 

180 PR INT:PR INT:PRINT:PR INT"NEUER START 

>: J/N > 7" 

170 GETGt:IFG$=""THEN170 
130 IFG$ = ”J"THENRUN 


64 



190 IFG$ ="N"THENPRINTCHR$< 147 ) "AUF WIEDE 
RSEHEN”SEND 
200 GOTO 170 

Es folgt eine Übungsaufgabe, die Sie vielleicht lösen möchten: Schreiben Sie ein 
kurzes Programm für Ihren Commodore 64. 

‘Schreibe zwei beliebige ganze Zahlen auf. Erzeuge eine Fibonacci-Folge durch 
paarweises Addieren von Gliedern zur Bildung eines weiteren Gliedes. Bestimme 
das Verhältnis jedes Gliedes zum vorhergehenden. Wohin tendiert das Verhältnis, 
wenn die Anzahl der Glieder groß wird? Berechne den Wert von 0.5 + SQR(5)/2.’ 
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KAPITEL 6 
ZAHLENBASEN 


Gewöhnlich schreiben wir Zahlen im Dezimalsystem auf. Zum Beispiel steht 1432, 
in Worten: Eintausendvierhundertzweiunddreißig, für den umständlicheren Aus¬ 
druck 

1*10 3 + 4*10 2 + 3*10 + 2. 

Eine noch etwas umständlichere Schreibweise wäre 


1*10 3 + 4*10 2 + 3*10' + 2*10° 


denn 10' = 10 und 10° = 1. Anders gesagt ist die Zahl 1432 als Summe von 
Vielfachen von Zehnerpotenzen zu interpretieren. Die ganzen Zahlen 1,4, 3 und 2 
heißen die Stellen der Zahl; 1 ist dabei die Tausenderstelle, 4 die Hunderterstelle, 3 
die Zehnerstelle und 2 die Einerstelle. Diese Darstellung der Zahl nennt man ihre 
Dezimaldarstellung; man sagt, daß die Zahl zur Basis 10 dargestellt ist. Das Wort 
dezimal ist abgeleitet vom lateinischen decem, zehn. 

Das Dezimalsystem hat 10 zur Basis. Man kann aber auch andere Basen als 10 
nehmen. Es ist nicht nur interessant, sondern auch nützlich, Zahlen mit verschiede¬ 
nen Basen auszudrücken. Beispielsweise hat sich die Darstellung von Zahlen zur 
Basis 2 als außerordentlich wichtig für Computer und alles damit Zusammenhän¬ 
gende erwiesen. 

Jede ganze Zahl größer als 1 kann als Basis dienen, und jede Zahl kann zu jeder 
beliebigen Basis dargestellt werden. Ihrem Commodore fällt es im übrigen leicht, 
Zahlendarstellungen von einer Basis auf eine andere umzurechnen. 

Sei N irgendeine positive Zahl und B eine ganze Zahl größer als 1. Um N zur Basis B 
auszudrücken, müssen wir N wie folgt schreiben: 

N = X m *B m + X m ^*B m -' + . .. + X^B + X„ 

wobei jede der Zahlen X 0 , X!, .... X m eine ganze Zahl zwischen 0 und B-1 ist. 
(Sehen Sie sich an, was herauskommt, wenn Sie 10 für B einsetzen). Die Stellen 
X 0 , X,, etc. heißen Koeffizienten der Zahl N zur Basis B. 

Kleine Werte der Basis B führen zu Zahlendarstellungen von großer Länge. Aber sie 
haben den Vorteil, daß sie mit weniger Koeffizienten auskommen. Der Extremfall ist 
B = 2. Das zugehörige System heißt binäres Zahlensystem (vom lateinischen 
binarius, zwei). Schreibt man eine Zahl im Binärsystem, können nur 0 und 1 als 
Koeffizienten Vorkommen. Zum Beispiel 
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86 = 64 + 16 + 4 + 2 

= 1 ’ 2 6 + 0 * 2 5 + 1 * 2 4 + 0 * 2 3 + 1 * 2 2 + 1*2 + 0 


Die Zahl 86 hat also die Binärstellung 1010110. Binärzahlen werden von Computern 
verwendet, weil sie als Folgen von Nullen und Einsen ausgedrückt werden können. 
Der Grund ist, daß 0 und 1 in einem Computer einfach den Stellungen ‘Aus’ und 
‘Ein’ eines Schalters entsprechen. 

Bei Basen größer als 10 braucht man zusätzliche Symbole. Es liegt nahe, die 
Buchstaben A, B, C etc. des Alphabets zu verwenden. Eine häufig benutzte Basis 
ist 16, Eine zur Basis 16 dargestellte Zahl heißt hexadezimal. Der Vorteil dieser 
Basis ist es, daß sie wenige Koeffizienten braucht und sich Hexadezimalzahlen 
trotzdem leicht in Binärzahlen umwandeln lassen. 

Die Umwandlung einer Zahl von der Basis 10 zur Basis B ist ganz einfach. 
Angenommen, wir wollen die Zahl N von der Basis 10 zur Basis B umwandeln. Erst 
subtrahieren wir von N alle Vielfachen von B. 

M = INT(N/B) : R = N - B * M 

Den Rest notieren wir und nennen ihn R 0 . Dann wiederholen wir das Verfahren mit 
M, indem wir N = M setzen. Den neuen Rest nennen wir R,. So fahren wir fort, bis 
der Wert von M schließlich 0 erreicht. Falls R s der letzte so erhaltene Rest ist, ist die 
ursprüngliche Zahl zur Basis B gleich 

R s . . . R 2 RiRo 

Wir wollen nun ein konkretes Beispiel durchrechnen. Angenommen, wir möchten 
die Zahl 29 zu Basis 3 schreiben. Dann rechnen wir wie folgt: 

Schritt 1. 

N = 29 
M = INT(29/3) 

= 9 

R 0 = 29 - 3*9 
= 2 

Schritt 2. 

N = M 

= 9 

M = INT(9/3) 

= 3 

Ri = 9 - 3*3 
= 0 
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Schritt 3. 

N = M 
= 3 

M = INT (3/3) 
= 1 

R 2 = 3 - 3*1 
= 0 


Schritt 4 
N 

M 

r 3 


= M 
= 1 

= INT (1/3) 
= 0 

= 1-3*0 
= 1 


Das Verfahren bricht nach vier Schritten ab, wenn M Null erreicht. Der Wert von 29 
zur Basis 3 ist also 1002. 

Das unten abgedruckte Programm wandelt ganze Zahlen von einer Basis in eine 
andere um. Beispielsweise kann man damit Zahlen von der Basis 10 auf die Basis 
16 umschreiben. Für Basen größer als 10 stellen die Buchstaben A, B, C etc. die 
Zahlen 10, 11, 12 etc. dar. 

10 REM 17BASIS UMWANDLER 

20 POKE53280,6sPRINTCHR$< 147)" ** 

* BAS IS-UMWANDLER ***" 

30 PRINT:PR INT"DIESES PROGRAMM WANDELT 2 
AHLEN VON EINERBASIS ZUR ANDEREN." 

40 PR INT s INPUT"UM2 UWANDELNDE BASIS “JA 
45 AA=54+As IFA< 1 1 THENAA =47+-A 
50 IFA< 20RA >35QRINT< A X >ATHENPRINT"NEUE 
EINGABE"s GOTO40 

60 PRINTS INPUT"UM2UWANDELNDE ZAHL "JM* 

65 N*=M* 

70 REM KONTROLLE N$ 

30 IFN*= ""THENPRINT"KEINE ZAHL"sGOTO60 

90 L =LEN< N$ > ■ I =0 

100 1=1+15N=ASC<MID$<N*,I,1>) 

110 IFN<4S0RCN>57 AND N<65> OR N>AA THEN 
PR INT"KE INE ZAHL"sGOTO60 
120 IF I < L THEN 100 
124 DIM A< L) 
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1£5 F0RI=1T0L:N=ASC<MID*<N$,I,1>>:IF N<5 
8 THEN A< I>=N-48 

186 IF N>64 THEN A<I>=N-55 

187 NEXT 

150 REM UMWANDLUNG VON BASIS ft ZU BASIS 
10 

160 N=VAL< N$ > 

170 I FA< >10THENN=0:FORI = 1TOL!N=A< I)+N*A: 
NEXT 

180 PRINTCHR$< 147 > s PR INT: PR INT"DEZIMALFO 
RM DER ZAHL: "N 

190 PRINT*. INPUT*'Z IELBAS IS "JB:PRINT 

800 IFB< >INT< B)ORB< 20RB >35THENPRINT"UNMO 

EGLICH-NEUER VERSUCH M :GOTO 190 

810 REM UMWANDLUNG N NACH B 

280 N$ = "" 

230 M=INT<N/B>:R=N-B*M:N=M 
240 IFR<10THENN$=CHR*<48+R>+N$ 

250 IFR >9THENN$ =CHR$< 55+R > +N$ 

8S0 IFN<>0THEN830 

870 PR INT:PR INTM$" MIT DER BAS IS"A" ERG I 

BT":PRINTN*" MIT DER BASIS"B 

280 PR INT:PR INT:PR INT"NEUER START <J/N>" 

890 GETA$: IFA*< >"N"AFOA$< >"J"THEN890 
300 IFA$ = ''J "THENRUN 

310 PRINTCHR*( 147) "AUF WIEDERSEHEN”:POKE 
53280,14:END 


64ER ZAHLEN 

Will man eine Zahl zwischen 0 und 255 als Binärzahl darstellen, braucht man 
höchstens 8 Koeffizienten, Zum Beispiel 


255 = 1*2 7 + 1*2 6 + 1*2 5 + 1 *2 4 + 1*2 3 + 1*2 2 + 1*2 + 1 
128 = 1*2 7 + 0*2 6 + 0*2 5 + 0’2 4 + 0*2 3 + 1*2 2 + 0*2 + 0 

Diese 8 Koeffizienten, oder 8 Bits auf einem Computer, nennt man ein Byte. Der 
Commodore 64 speichert ganze Zahlen in zwei Bytes, die High Byte und Low Byte 
genannt werden. Das High Byte stellt Vielfache von 256 dar. Beispielsweise würde 
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die Zahl 999 mit High Byte 3 und Low Byte 231 gespeichert, denn 999 = 3*256 + 
231. Zahlen in Bytes zu speichern bedeutet nichts anderes als Zahlen zur Basis 256 
zu speichern. 

Durch PEEKen können Sie sehen, wie der 64 High Byte und Low Byte einer ganzen 
Zahl speichert. Erinnern Sie sich, daß Integervariablen auf dem Commodore 64 
durch das Prozentzeichen (%) hinter dem Variablennamen gekennzeichnet sind. 
Geben Sie folgendes ein, und drücken Sie die Return-Taste am Ende jeder Zeile: 

NEW : CLR 

X% = 999 (beliebige ganze Zahl eingeben) 

PRINT “HIGH BYTE=“ PEEK(2053), “LOW BYTE=” PEEK (2054) 

Die größte ganze Zahl, die der Commodore 64 speichern kann, ist 32767, also 
127*256 + 255. Zahlen mit einem High Byte von 128 oder größer sind negative 
Zahlen. Negative Zahlen werden nämlich gespeichert, indem man erst den ABSo- 
lutbetrag der Zahl betrachtet, High und Low Byte berechnet und dann den Wert des 
High Byte von 255 subtrahiert und den des Low Byte von 256 subtrahiert. Auf diese 
Weise erhält z. B. -1 ein High Byte von 255 und ein Low Byte von 255. 

Eine Zahl wird wie folgt aus High und Low Byte berechnet. Bezeichne H das High 
Byte und L das Low Byte. 

NUMBER = H*256 + L 

IF H > = 128 THEN NUMER = —((255—H)*256 + (256—L)) 

Die letzte Zeile kann auch so geschrieben werden: 

IF H > = 128 THEN NUMBER = H*256 + L - 256*256 

Versuchen Sie, Zahlen in die Adressen 2053 und 2054 zu POKEn, geben Sie 
PRINT X% ein und vergleichen Sie die Antworten. Tippen Sie z. B. folgendes ein, 
und drücken Sie die Return-Taste am Ende jeder Zeile: 

NEW : CLR : X% = 0 

POKE 2053,98 : POKE 2054,99 (Variableninhalt ersetzen) 

N = PEEK(2053)*256 + PEEK(2054) 

IF PEEK(2053) > = 128 THEN N = N - 256*256 
PRINT X%, N 
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KLEINE ZAHLEN 

Das Programm im vorletzten Abschnitt verarbeitet positive ganze Zahlen. Jede Zahl, 
ob ganz oder nicht ganz, hat jedoch eine Darstellung in jeder Basis. Beispielsweise 
wird die Dezimalzahl 0.25 binär als 0.01 dargestellt, und die Dezimalzahl 0.125 binär 
als 0.001. Um das nachzuvollziehen, wollen wir uns zunächst klarmachen, was die 
Dezimalzahl 0.25 bedeutet. Diese Zahl steht für 2 Zehntel und 5 Hundertstel, also 

0.25 = 2/10 + 5/100 
= 2*10 1 + 5 * 10 2 . 

Um sie zur Basis B da’zustellen, müssen wir sie in dieser Form schreiben: 

Y^B ’ + Y 2 *B 2 + . . . 


wobei wie üblich B - ' 1 =1/B, B 2 =1/B 2 , etc. 

Die Dezimalzahlen 0.25 und 0.125 können wie folgt geschrieben werden: 

0.25 = 1/4 

= 0 * 2 -’ + 1*2 2 
0.125 =- 1-8 

= 0*2 1 + 0*2 2 + 1*2 3 

Daraus erklärt sich die Binärdarstellung dieser Zahlen. 

Sehen Sie als weiteres Beispiel, wie die Zahl 0.6 in (negativen) Potenzen von 2 
dargestellt wird: 

0.6 = 1*2 1 + 0 * 2~ 2 + 0 * 2- 3 + 1*2 4 
+ 1*2 5 + 0*2 6 + 0 * 2 7 + 1 * 2 8 
+ 1 * 2 9 + 0*2 10 + 0*2 " + 1*2 12 
+ 1 * 2~ 13 + 0*2 14 + 0 * 2' 5 + 1*2 16 
+ . . . 

Tatsächlich braucht man unendlich viele Stellen, um 0.6 exakt binär aufzuschreiben. 
Es ergibt sich die Darstellung 

0.10011001100110011001100110011001 . . . 

Der Commodore 64 speichert nur 32 dieser Stellen, angefangen bei der ersten von 
Null verschiedenen; zusätzlich rundet er auf, wenn die 33. signifikante Stelle 
ungleich Null ist. Also speichert der 64 die Zahl 0.6 binär als 

0.10011001100110011001100110011010 
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Das nächste Programm druckt die Binärdarstellung einer Dezimalzahl zwischen 0 
und 1 so aus, wie sie vom Commodore 64 gespeichert wird. 

10 REM 18BINAER UMWANDLER 

20 POKE53280,6sPRINTCHR$( 147 )" *** UMWAN 
DLUNG DEZIMAL ZU BINAER ***“ 

30 PRINTs PR INT"DIESES PROGRAMM ZEIGT DIE 
BINAERFORM EINER ZAHL ZWISCHEN m f 
40 PR INT"NULL UND EINS 
50 N=0.'PRINTs INPUT"ZAHL E INGEBEN“ ; N 
60 IFN<=0ORN>=1THENPRINTSPRINT"!!ZWISCHE 
N NULL UND EINS!!"sGOTO50 
70 N$ = "0" 

30 FOR 1 = 1T032 

30 N=N*2SN$=N$+MID$< STR*C INTC N> ) , 2,1 > : N= 

N- INT< N> 

100 NEXT I 

110 PRINTCHR$< 147)sPR INTsPR INT"DIE BINAE 
RFGRM DER ZAHL ISTs" 

115 PRINTsPRINTN* 

120 PR INT s PR INT s PR INT s PR INT"NEUER START 
<J/N> ?'* 

130 GET GfslF G$ = "" THEN GOTO130 
140 IF G*="J” THEN RUN 

150 IF G$ ="N" THENPRINTCHR$< 147)"AUF WIE 
DERSEHEN" : EtO 
160 GOTO 130 


GLEITKOMMAZAHLEN 

Ganze Zahlen werden im Commodore 64 in zwei Bytes gespeichert. Aber reelle 
Zahlen werden in 5 Bytes gespeichert, selbst wenn ihr Wert eine ganze Zahl ist. 
Wenn Sie eine Zahl nicht mit dem Prozentzeichen zur ganzen Zahl erklären, wird sie 
als reelle Zahl in 5 Bytes gespeichert. Eine Zahl kann binär wie folgt dargestellt 
werden: 


1.X 1 X 2 X 3 ...X m *2 N 
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wobei X,, X 2 . . X m entweder 0 oder 1 sind und N eine ganze Zahl ist (positiv, 
negativ oder Null), Die ganze Zahl N heißt binärer Exponent der Zahl, der Rest heißt 
binäre Mantisse. Beispielsweise ist die Dezimalzahl 10 binär 1010 und läßt sich 
schreiben als 

1.01 * 2 3 

Also hat 10 den binären Exponenten 3 und die binäre Mantisse 1.01. Ein anderes 
Beispiel ist die Dezimalzahl 0.375, die binär 0.011 ist und sich schreiben läßt als 

1.1 * 2 2 

Also hat die Dezimalzahl 0.375 den binären Exponenten -2 und die binäre 
Mantisse 1.1. 

Wir haben erwähnt, daß der 64 zum Speichern von Zahlen 5 Bytes benutzt. Das 
erste Byte ist der binäre Exponent plus 129. Die übrigen vier Bytes bestimmen die 
binäre Mantisse und das Vorzeichen einer Zahl. Da die erste Stelle einer binären 
Mantisse immer 1 ist, braucht sie nicht gespeichert zu werden - es reicht, die 
Stellen rechts von Dezimalpunkt der binären Mantisse zu speichern. Das erste Bit 
von Byte 2 bestimmt das Vorzeichen der Zahl, die übrigen 31 Bit der letzten vier 
Bytes enthalten die binäre Mantisse (ohne die führende 1). 

Beispielsweise würde die Dezimalzahl 10 wie folgt gespeichert: Das erste Byte ist 
129 plus den binären Exponenten 3, also 132. Das erste Bit des zweiten Byte ist 0, 
da die Zahl positiv ist. Die übrigen 31 Bit wären 

010000000000000000000000000000 

weil 10 die binäre Mantisse 1.01 hat und die führende 1 ignoriert wird. Die 32 Bit in 
den letzten vier Bytes wären also 

0010000000000000000000000000000 

Aufgeteilt in vier Blöcke zu 8 Bit ergibt das 

00100000 00000000 00000000 00000000 

was wiederum 32, 0, 0, 0 ist. Also würde die Dezimalzahl 10 in den 5 Bytes 132, 32, 
0, 0, 0 gespeichert. 

Das Verfahren läßt sich umkehren, um die Zahl zu ermitteln, die der 64 in 5 Bytes 
abgelegt hat. Nehmen wir an, die Zahl N ist in den 5 Bytes P, Q, R, S, T gespeichert. 
Folgende Programmzeiien berechnen N aus P, Q, R, S und T: 
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X = 1 : IF Q > = 128 THEN Q = Q-128 : X = - 1 
N = X * 2 p -' 29 * (1 + Q*2 7 + R*2 ' 5 + S*2~ 23 + T*2 31 ) 

Um den 64 in Aktion zu erleben, geben Sie diese Zeilen ein, und drücken Sie 
Return am Ende jeder Zeile: 

NEW : CLR 

X = 10 (beliebige Zahl eingeben) 

FOR 1=0 TO 4 : PRINT “BYTE’’ 1 + 1 “ = ” PEEK(2053+I) : NEXT I 

Sie werden allerdings feststellen, daß der Commodore 64 beim Multiplizieren 
gelegentlich kleine Fehler macht. Die Zahl 1 + 2 24 wird (korrekt) in den folgenden 5 
Bytes gespeichert: 

129, 0, 0, 0, 128 

Wird dieselbe Zahl jedoch als 1 *(1 + 2~ 24 ) geschrieben, wird sie in diesen 5 Bytes 
gespeichert: 

129, 0, 0, 0, 64 

Anders gesagt erhält der Commodore 64 das Ergebnis 

1 * (1 + 21—24) = 1 + 2|—25 

Dasselbe Problem läßt sich auch so sichtbar machen: 

X = 1 + 2|-24 
PRINT X - X, X - 1*X 

Versuchen Sie es auch mit 

X = 1 + 2|—24 

XI = 1*X 
X2 = 1 *X1 
X3 = 1*X2 
X4 = 1*X3 

PRINT X, XI, X2, X3, X4 

Probleme dieser Art scheinen nicht aufzutauchen, wenn die mittleren 3 Bytes der 
zum Speichern der Zahl verwendeten 5 Bytes nicht alle gleich 0 sind. Versuchen 
Sie, Beispiele zu finden. 
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KAPITEL 7 


TAGE UND WOCHEN 


TAGE 

Die Zellersche Kongruenz ist eine kompliziert aussehende Formel, die den 
Wochentag (Sonntag, Montag, etc.) zu jedem gegebenen Datum bestimmt. Mit 
dieser Formel können Sie z. B. herausfinden, an welchem Wochentag jemand 
geboren ist. Und falls Sie ein schlechtes Gedächtnis haben, könnten Sie berech¬ 
nen, auf welchen Wochentag ein bestimmtes Jubiläum gefallen ist. 

Die Zellersche Formel sieht so aus: 

A = INT(2.6*M - 0.1) + D + Y + INT(Y/4) + INT(C/4) -2*C 
X = A - 7*INT(A/7) 

Die Zahl X liegt zwischen 0 und 6, weil alle Vielfachen von 7, die kleiner als A sind, 
von A subtrahiert worden sind. Diese Zahlen entsprechen den 7 Wochentaqen wie 
folgt: 

0 : Sonntag 

1 : Montag 

2 : Dienstag 

3 : Mittwoch 

4 : Donnerstag 

5 : Freitag 

6 : Samstag 

Die Zahlen D, M, Y und C sind definiert als: 

D : der Monatstag. 

M : die Monatszahl - aber nicht die übliche. Januar und Februar tragen 
die Zahlen 11 und 12 des vorhergehenden Jahren (das beeinflußt das 
unten beschriebene Y und gelegentlich auch C). März trägt die Zahl 1, 
April die Zahl 2, Mai die Zahl 3, . . . und Dezember die Zahl 10. 

Y : das Jahr im Jahrhundert. 

C : die Zahl der Jahrhunderte, d. h. die ersten zwei Stellen der Jahreszahl. 

Für den 26. August 1983 ist beispielsweise die übliche Schreibweise 26/08/1983. 
Für die Zellersche Formel verwenden wir D = 26, M = 6, Y = 83 und C = 19. 
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Durch Einsetzen dieser Zahlen in die Zellersche Formel ergibt sich 


A 


und 

Y 


INT(2.6*6 - 0.1) + 26 + 83 + INT(83/4) + INT(19/4) - 2*19 
INT(15.6 - 0.1) + 26 + 83 + INT(20.75) + INT(4.75) - 38 
15 + 26 + 83 + 20 + 4 - 38 
110 


110 - 7*INT(110/7) 

110 - 7*INT(15.714285) 
110 - 105 
5. 


Wir schließen daraus, daß der Wochentag des 26. August 1983 ein Freitag ist. 
Hier sind einige Beispiele für die übliche Schreibweise des Datums und die Werte, 
die in die Zellersche Formel eingehen. 

ÜBLICHE SCHREIBWEISE SCHREIBWEISE FÜR DIE ZELLERSCHE 

FORMEL 



D 

M 

C 

Y 

03/03/1947 

3 

1 

19 

47 

01/01/2000 

1 

11 

19 

99 

26/02/1983 

26 

12 

19 

82 

29/11/1984 

29 

9 

19 

84 


Das nachstehende Programm benutzt die Zellersche Formel, um den Wochentag zu 
jedem vorgegebenen Datum auszurechnen. Ihr Commodore 64 bestimmt automa¬ 
tisch die richtigen Werte von D, M, C und Y zu einem beliebigen eingegebenen 
Datum. Beachten Sie, daß das Programm mit einigen Kontrollen sicherstellt, daß 
das eingegebene Datum sinnvoll ist. Beispielsweise wird der 30. Februar 1983 nicht 
akzeptiert. Außerdem muß das eingegebene Jahr eine ganze Zahl im Bereich von 
1752 bis 4902 sein. Die Zellersche Formel bezieht sich auf den Bereich von 1582 
bis 4902, aber der Gregorianische Kalender ist in Großbritannien, den Britischen 
Kolonien und den USA erst seit 1752 in Gebrauch. 

Schaltjahre werden automatisch vom Programm berücksichtigt. Beachten Sie, daß 
ein Jahr ein Schaltjahr ist, wenn die Jahreszahl ohne Rest durch 4 teilbar ist, es sei 
denn, sie ist durch 100, aber nicht durch 400 teilbar. Also war 1900 kein Schaltjahr, 
aber 2000 wird eines sein. 


78 



10 REM 19W0CHENTAGE 

20 DIM A< 1£),A$< 1 £ ) sFOR1 = 1TO 1£sREAD A< I) 
,A$< I):NEXTI 

30 FORI=0TO6iREADB$<I)sNEXTI 
40 POKE53280,8 s PR INTCHR$< 147)" 

*** WOCHENTAGE ***" 

50 PR INT SPRINT "DIESES PROGRAPW ZEIGT DEN 
WOCHENTAG IR- GENDEINES "; 

60 PRINT"SPEZIELLEN TAGES" 

70 PR INT!PR INT"BITTE DATUM EINGEBEN" 

30 A< £ ) =£9:D=0:PR INT! INPUT“TAG ( ZW. 1 UND 
3 1 ) " ; D 

90 IF D< 1 OR D >31 OR D<>INT<0) THENPRINT 

"NOCH MAL VERSUCHEN":GOTO80 

100 Y = 0:PR INT: INPUT"MONAT ( ZW. 1 UND 12) 

" ; M 

110 I FM< 1ORM > 1 £ORM< > INT( M ) THENPR I NT “NOCH 
MAL VERSUCHEN"!GOTO 100 
120 IF D >A< M)THENPRINT!PR INT"ZUVIEL TAGE 
IM MONAT"s FORI = 1TO2000!NEXT:GOTO40 
130 Y = 0:PRINT: INPUT"JAHR <VIERSTELL IG) 

« . Y 

140 IF Y< 1582 OR Y>490£ OR YOINT(Y) THE 

NPRINT"NOCH MAL VERSUCHEN":GOTO 130 

150 REM SCHALTJAHR 

160 L=0! IF INT<Y/4)*4=Y THENL = -1 

170 IFLANDINTC Y/1 00 ) * 100 = YTHENL =0 ! IFINT< 

Y/400 > *400 = YTHENL = -1 

175 PRINTCHR$< 147) 

180 A< £)=£8-L: IFLTHENPRINTSPRINTSPRINTY“ 
IST/WAR EIN SCHALTJAHR" 

190 IF M=£ANDD>AC2)THENPRINTSPRINTSPRINT 
" M £8 TAGE IM FEBRUAR M"sRT=l 
195 IFRT=ITHEND=D-1 

£00 PR INT! PR INTs PR INT" DER “D ”||. "M"H. " Y" I 
ST/'WAR EIN "; 

£10 M=M-£: IFM< 1THENM=M+1£ s Y=Y-1 

££0 C= INT< Y/' 100 ) ! Y = Y-C* 100 

230 REM ZELLER’SCHE KONGRUENZ 

£40 DAY=INT<£.6 +M-0. 1)+D + Y+INT( C/4) + INT< 

Y/4)-2+C 
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250 DAY=DAY-7*INT<DAY/7 >:PRINTB$< DAY> 

260 PR INT:PR INT s PR INT" NEUER START ( J/N) 

7 " 

270 GET G$:IF G$ = " " THEN270 
230 IFG$ = " J "THENRUN 

290 IFG$ = "N”THENPRINTCHR$< 147 ) "AUF WIEDE 
RSEHEN”:END 
300 GOTO270 
400 REM DATA 

410 DATA 31,JANUAR,29,FEBRUAR,31,MAERZ,3 
0,APRIL 

420 DATA 31,MAI,30, JUNI,31,JULY,31,AUGUS 
T 

430 DATA 30,SPTEMBER,31/OKTOBER,30,NOVEM 
BER,31 ,DEZEMBER 

440 DATA SONNTAG,MONTAG,DIENSTAG,MITTWOC 
H,DONNERSTAG,FREITAG,SAMSTAG 

Anmerkung: Sie haben vielleicht die Zeliersche Formel vorher schon einmal gese¬ 
hen und möglicherweise bemerkt, daß die hier verwendete Formel etwas davon 
abweicht. Häufig ist der erste Ausdruck INT(2.6*M - 0.1) durch den Ausdruck INT 
(2.6*M - 0.2) ersetzt. Daß der zweite Ausdruck hier nicht benutzt wird, hat seinen 
Grund in dem Verfahren, mit dem der Commodore den ganzen Teil INT einer Zahl 
berechnet. Für M = 7 ergibt sich zum Beispiel: 

INT(2.6*7 - 0.2) = INT(18.2 - 0.2) 

= INT (18) 

= 18 

Der Commodore 64 weist jedoch INT(2.6*7 — 0.2) den Wert 17 zu, obwohl er 
korrekt 18 ausdruckt, wenn Sie ihn nach 2.6*7 - 0.2 fragen. 


KALENDER 

Sobald wir den Wochentag für jedes Datum kennen, können wir einen Kalender 
herstellen. Das nächste Programm druckt einen Kalender für jeden Monat eines 
beliebigen Jahres. (Auf dem Bildschirm kann man nur einen Monat übersichtlich 
darstellen.) Das Programm berechnet mit der Zellerschen Formel den Wochentag, 
auf den der erste Tag des Monats fällt. Wie im obigen Programm für den Wochentag 
werden Schaltjahre automatisch berücksichtigt. 
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10 REM 20KALENDER 

20 DIMA( 12) ,A$( 12 ) s FOR 1 = 1 TO 12 s RE ADAC I ) ,A 
$( I>JNEXTI 
30 REM START 

40 POKE53280,6sP0KE53281,6 s PRINTCHR$( 147 
>" *** KALENDER ***" 

50 PRINTs PR INT"DIESES PROGRAMM 2EIGT EIN 
EN KALENDER VONJEDEM MONAT "i 
60 PRINT"EINES GEWUENSCHTEN JAHRES." 

70 PR INT s PR INT"BITTE MONAT UND JAHR EING 
EBEN" 

80 M=0:PR INT! INPUT”MONAT (ZW. 1 UND 12)" 
* M 

90 IFM< 10RM>120RM< >INT( M)THENPRINT"NOCH 
MAL VERSUCHEN"JGOTO80 

100 Y=0:PRINTsINPUT"JAHR (VIERSTELLIG) " 
; Y 

110 IFY< 1752ORY>490£ORY< >INT( Y)THENPRINT 

"NOCH MAL VERSUCHEN"sGOTO 100 

130 REM SCHALTJAHRE 

140 L =0 : IF I NT( Y/4 ) *4 = YTHENL = - 1 

150 IFLANDINT( Y/100)* 100 = YTHENL=0 s IFINT( 

Y/400)*400=YTHENL=-1 

155 PRINTCHR#( 147)CHRt( 31)sP0KE53281,7sp 
OKE53280,7 

157 PR INTs PR INTA$ ( M))Y 

160 A(2 ) =28-L: IFLTHENPRINT:PRINT“DIES IS 
T EIN SCHALTJAHR" 

170 A$ = A$(M)+STR$( Y)s MM=A< M) 

210 M=M-2: IFM< 1THENM=M+1 £ sY=Y-1 

220 C= INT( Y/100 ) : Y = Y-C*. 100 

£30 REM ZELLER'SCHE KONGRUENZ 

240 DAY= INT( 2.6*M-0.1 ) + 1 +Y+INT( C/4 ) + INT( 

Y/4)-2*C 

£50 DAY=DAY-7 * INT( DAY/7) 

260 REM KALENDER ZEIGEN 
£80 PRINTTAB( £0-LEN( A$)/£) 

290 PR INT:PR INT s PR INTCHR$ < 28)“ SO MO 

DI MI DO FR SA "CHRt(31) 

295 PRINT 
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300 FQR1=1TOMM 
3 10 DAY=DAY+1 

320 PRINTTAB<DAY*5-2+<I>9>> I; 

330 IFDAY=7THENDAY=0sPR INT:PR INT 
340 NEXTI 

350 PR INT:PR I NT:PR INTs PR INT"NEUER START 

<J/N> ?" 

360 GETG*:IFG*="“THEN3S0 

370 IFG$ = "J"THENPRINTCHR$< 154>:RUN 

380 IFG$ = "N"THENPRINTCHR$< 154>:P0KE53281 

,6:POKE53230,14:PRINTCHRSC 147) 

385 IFG*="N"THENPRINT"AUF WIEDERSEHEN":E 
ND 

330 GOTO3G0 
400 REM DATA 

410 DATA 31,JANUAR,29,FEBRUAR,31„MAERZ,3 
0,APR IL 

420 DATA 31,MAI,30,JUN1,31, JUL Y,31,AUGUS 
T 

430 DATA 30,SEPTEMBER,31,OKTOBER,30,NOVE 
PEER,31,DEZEMBER 


TERMINKALENDER 

Man braucht gelegentlich eine Liste von Terminen, die eine feste Anzahl von Tagen 
auseinanderliegen, z. B. Behandlungstermine im Krankenhaus oder Zahltage. 

Um eine solche Liste anzufertigen, verwenden wir das ‘pseudo-Julianische’ Datum. 
Das ist einfach die Zahl der Tage seit einem beliebigen festen Datum. (Zum Beispiel 
hat der erste Januar des Jahres 1 das ‘pseudo-Julianische’ Datum 1.) Eine 
verhältnismäßig einfache Formel wandelt das tatsächliche Datum in ein pseudo- 
Julianisches Datum um und umgekehrt. 

Sei D/M/Y das Datum, wobei D den Tag, M den Monat und Y das Jahr (einschließ¬ 
lich des Jahrhunderts) bezeichnet. Dann berechnet sich das pseudo-Julianische 
Datum wie folgt: 

X = INT(30.57*M) + INT(365.25*Y - 395.25) + D 
Ist M> 2 und Y ein Schaltjahr, subtrahiere 1 von X. 

Ist M> 2 und Y kein Schaltjahr, subtrahiere 2 von X. 
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Als Beispiel wollen wir das pseudo-Julianische Datum des 26. August 1983 
berechnen. Die Werte von D, M und Y sind durch D = 26, M = 8, Y = 1983 
gegeben. In die Formel eingesetzt liefern diese Werte: 

X = INT(30.57*8) + INT(365.25*1983 - 395.25) + 26 
= INT(244.56) + INT(723895.5) + 26 
= 244 + 723895 + 26 
= 724165 

Da aber die Monatszahl M größer als 2 ist und 1983 kein Schaltjahr ist, ziehen wir 2 
von X ab und erhalten das pseudo-Julianische Datum 724163. 

Um umgekehrt das Datum aus dem pseudo-Julianischen Datum zu bestimmen, 
gehen wir wie folgt vor; sei X das pseudo-Julianische Datum: 

Die erste Approximation für das Jahr ist: 

Y = INT(X/365.26) + 1 

Der Jahrestag ergibt sich als: 

D = X - INT(365.25*Y - 395.25) 

Eine Schaltjahr-Korrektur wird vorgenommen: 

Dl = 2; ist Y ein Schaltjahr, setze Dl = 1 
Ist D> 91-Dl, addiere Dl zu D. 

Monat und Tage werden berechnet: 

M = INT(D/30.57) 

D = D - INT(30.57*M) 

Monat und Jahr werden nötigenfalls korrigiert: 

Ist M> 12, setze M = 1 und addiere 1 zu Y. 

Wir möchten beispielsweise das Datum bestimmen, das dem pseudo-Julianischen 
Datum 724164 entspricht (das ist 1 höher als das vorhin berechnete pseudo- 
Julianische Datum). Folgende Berechnungen sind nötig: 

Y = INT(724164/365.26) + 1 
= INT(1982.5987) + 1 

= 1983 

D = 724164 - INT(365.25*1983 - 395.25) 

= 724164 - INT(724290.75 - 395.25) 

= 724164 - INT(723895.5) 

= 269 
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Das Jahr 1983 ist kein Schaltjahr, also ist Dl gleich 2. Weil der Wert von D größer 
als 91 - Dl ist, addieren wir Dl zu D. Es ergibt sich daraus 


D = 271 

M = INT(271/30.57) 

= INT(8.86490023) 

= 8 

D = 271 - INT(30.57*8) 

= 271 - INT(244.56) 

= 271 - 244 
= 27 

Da der Wert von M nicht größer als 12 ist, sind wir fertig: D = 27, M = 8 und Y = 
1983. Also entspricht dem pseudo-Julianischen Datum 724164 der 27. August 
1983. 

Das folgende Programm erledigt in kurzer Zeit alle diese Berechnungen und fertigt 
eine Liste von Terminen an, die eine vorgegebene Zahl von Tagen auseinanderlie¬ 
gen. Der Einfachheit halber funktioniert das Programm nur für Daten im 20. 
Jahrhundert. Die nötigen Änderungen für ein anderes Jahrhundert sollten Ihnen 
nicht schwerfallen. 

Achtung: Daten werden in der Form DD/MM/YY eingegeben; z. B. würde 12. März 
1984 als 12/03/84 oder 12/ 3/84 eingegeben, aber nicht als 12/3/84. Das Datum 
wird nach der Eingabe nicht gerade vielen Kontrollen unterzogen, daher könnten 
Sie beispielsweise 30/02/83 eingeben, was vom Programm als 2. März 1983 
interpretiert würde. (Erkennen Sie warum?) 

10 REM £ 1DATEN MANAGEMENT 

20 POKE53H80,6sPRINTCHRf< 147)" ** 

* TERMIN MANAGEMENT ***" 

30 PR INT s PR INT"DIESES PROGRAMM LIEFERT T 
ERMINE EINER BESTIMMTEN ZEITSPANNE" 

40 REM EINGABE ANFANGSDATUM 

60 PR INTSPRINTSPRINT"EINGABE DES ANFANGS 

DATUM <TT/MM/JJ>" 

80 ..PRINTsINPUT"ANFANGSDATUM s "; S* 

30 IFLEFKS$)< >8THENPRINT SPRINT"IN FORM T 
T/MM/JJ"s GOTO80 

100 IF MID*<S$,3,1 X >"/"OR MID$< S$,6,1 )< 

>"/"THENPRINT s PR INT"IN FORM TT/MM/JJ°sGO 
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1 10 T=VAL< MIDfC S$ , 1 ,2) > : M=VALC MID*C S* ,4 , 
2)> : J=VALCMID*< S$,7,2)>+1900 

120 IFT<=0ORT>31ORM<=0ORM >12THENPRINT"TA 
G, MONAT FALSCH"SGOTO80 

130 REM DATUM IST IN DER KORREKTEN FORM 
140 REM BERECHNUNG DES PSEUDO-JULIAN DAY 

150 X = INTC 30.57*M)+INTC 365.25*J-395.25) + 
T 

160 REM SCHALTJAHR 

170 IFM>2THENX=X-2sIF INTCJ/4)*4=J THEN 
K=X+1 

180 REM EINGABE ZEITSPANNE 

190 PR INT : PR INT : PR INT"EINGABE DER ZEITSP 

ANNE IN TAGEN” 

200 P=0:PRINTsINPUT"ZEITSPANNE : ";P 

2 10 IF P< =0 OR INTCPX >P THENPR INT: PR INT 
"GANZE ZAHLEN BITTE"s GOTO200 

220 PR INTsPR INTSPRINT"ANZAHL DER ZEITSPA 
NNEN" 

230 PRINT:INPUT"ZAHL : ";N 

240 IF N<=0 OR INTC NX >N THENPR INT: PR INT 

"GANZE ZAHLEN BITTE":GOTO230 

250 IF N< =0 OR N> 100 THEN PRINT"VERNUENF 

TIG BLEIBEN":GOTO230 

260 PRINTCHR$C 147)CHR$C5) 

270 FOR I=0TO N 

280 :J =INTC X/365.26) + 1 : REM JAHR 
230 : T=.X- INTC 365.25*J-395.25 ) 

300 :REM SCHALTJAHR 

310 :T1=2:IF INTCJ/4)*4*J THENT1=1 

320 :IFT>91-T1 THEN T=T+T1 

330 :M =INTCT/30.57) :REM MONAT 

340 s T=T-INTC 30.57#M>:REM TAG 

350 : IF M >12 THENM=1:J =J +1 

360 :J =J-1300 

370 :REM AUSGABE 

380 :Z =T:GOSUB 500:T$ =Z$ + "/" 

390 :Z=M:GOSUB 500: T$=T$+Z* +"/" 

400 :Z=J:GOSUB 500:Tf=T*+Z$ 
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4 10 ’.PRINTT* 

420 :X=X+P 
430 NEKT 

440 PR1NT!PRINT:PRINT:PRINTCHR$( 154)"NEU 
ER START <J/N) ?" 

450 GETGf s IFG$ = “"THEN450 
460 IFG$ = ”J"THENRUN 

470 IFG* = ' , N ,, THENPRINTCHR$< 147) "AUF WIEDE 

RSEHEN " :END 

480 GOTO450 

500 REM FORMATIERUNG 

5 10 Z$=MID$<STR*< Z),2)s IF LEN< Zf >< 2THENZ 
$= »0 "+ 2 $ 

580 RETURN 
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KAPITEL 8 

DER GRÖSSTE GEMEINSAME TEILER 


Wenn A und B ganze Zahlen sind, dann ist ein gemeinsamer Teiler oder gemeinsa¬ 
mer Faktor von A und B eine ganze Zahl, die beide Zahlen teilt. Der größte Faktor 
von A und B ist die größte solcher ganzen Zahlen. 

Zum Beispiel ist 3 ein gemeinsamer Teiler von 12 und 18, aber der größte 
gemeinsame Teiler von 12 und 18 ist 6. 

Das Berechnen des größten gemeinsamen Teilers zweier Zahlen ist nicht sonder¬ 
lich kompliziert, vor allem nicht für einen Computer. Die dabei verwendete Methode 
läßt sich gut zur Veranschaulichung eines Rechenalgorithmus benutzen. 

Der Euklidische Algorithmus ist die bekannteste und älteste (Drittes Jahrhundert v. 
Chr.) Methode zur Berechnung des größten gemeinsamen Teilers. Wenn man den 
größten gemeinsamen Teiler von A und B finden will, verfährt man wie folgt: 

1. Benenne A und B (wenn nötig) um, so daß A größer ist als B. 

2. Dividiere A und B und bestimme den Rest R v 

Rt = A - B * INT(A/B) 

Man beachte, daß jede Zahl, die A und B teilt, auch R, teilt. Und umgekehrt ist jeder 
gemeinsame Teiler von B und R, auch ein Teiler von A. Daraus ergibt sich, daß die 
gemeinsamen Teiler von A und B dieselben sind wie die von B und R,. Also ist der 
größte gemeinsame Teiler von A und B gleich dem größten gemeinsamen Teiler 
von B und R r . 

3. Teile nun B durch R, und bestimme den Rest R 2 . 

R 2 = B - Ri * INTfB/R,) 

Was oben zu den Zahlen B, R, bemerkt wurde, gilt auch für R,, R 2 . Also ist der 
größte gemeinsame Teiler von A und B gleich dem größten gemeinsamen Teiler 
von R! und R 2 . 

4. Als nächstes teile R 2 durch R, und bestimme den Rest R 3 . 

R 3 = Ri - R 2 * INT(R 1 /R 2 ) 

Das Verfahren wird in dieser Weise fortgesetzt, bis der Rest Null ist. Beachten Sie, 
daß die Reste jedesmal kleiner werden und deshalb nach einer Anzahl von Schritten 
Null erreichen. 
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R, = A - B * INT(A/B) 

R 2 = B - R, * INTtB/RO 
r 3 = Ri - R 2 * INT(R 1 /R 2 ) 

• • 

• • 

• • 

Rn-1 = Rn-3 — Rn-2 * INT(R N _3/Rn- 2 ) 
Rn = Rn-2 ~ Rn-i * INT(R n _ 2 /Rn-i) 


(0 <= Ri<B) 
(0 <= R 2 <Ri) 

(0 <= r 3 <r 2 ) 


(0 <— Rn-i Rn-2) 
(Rn = 0) 


Hat der Rest 0 erreicht, so sieht man, daß der vorhergehende Rest R N _, der größte 
gemeinsame Teiler von Rn-i und Rn -2 ist. Führt man dieses Argument aus, so 
erkennt man R N -i als größten gemeinsamen Teiler von A und B. 

Der oben skizzierte Prozeß läßt sich leicht auf den Computer übertragen. Das leistet 
das folgende Programm: 


10 REM S2GR0ESSTER GEMEINSAMER TEILER 
£0 PQKE53280,6sPRINTCHR$< 147)" *** GROES 

STER GEMEINSAMER TEILER ***" 

30 PR INT!PR INT"DIESES PROGRAMM BERECHNET 
DEN GROESSTEN GEMEINSAMEN TEILER ZWEIER 

M • 

r 

40 PRINT" GANZER ZAHLEN UNTER ZUHILFENAH 
ME DES EUCLIDSCHEN ALOGORITHMUS" 

50 PR INT:PRINT"EINGABE ZWEIER GANZER ZAH 
LEN" 

60 PR INT! INPUT"ERSTE ZAHL ";A 

70 IFA< 1 QRA< >INT< A)THENPRINT!PR INT"NOCH 

MAL VERSUCHEN"GOTO80 

SO PR INT’. INPUT"ZWE ITE ZAHL “;B 

90 IFA< 10RA< >INT<AJTHENPRINTiPR INT"NOCH 

MAL VERSUCHEN"GOTO30 

100 IFA<BTHENC=A!A=BiB=C 

110 REM EUCLIDSCHER ALOGORITHMUS 

120 R =A s S =B 

130 T=R-S*INT<R/S> 

140 IFT< >0THENR=S:S=T:GOTO 130 
150 PRINTCHR*< 147>!PR INTsPR INT!PR INT"GRO 
ESSTER GEMEINSAMER TEILER ! " ; S 
160 PR I NT s PR I NT"ERSTES GEMEINSAMES VIELF 
ACHES:"A*B/S 

170 PR INT:PR INT:PR INT:PR INT"NEUER START 
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•; j/n> 7" 

180 GETG* s IFG$ = “"THEN180 
130 IFG$ = 11 J "THENRUN 

200 1FG$ = "N"THENPRINTCHR$< 147) "AUF WIEDE 
RSEHEN"SEND 
210 GOTO 180 

Das Programm berechnet auch das kleinste gemeinsame Vielfache von A und B. 
Das kleinste gemeinsame Vielfache zweier Zahlen A und B ist die kleinste Zahl, die 
sich sowohl durch A als auch durch B teilen läßt. Das kleinste gemeinsame 
Vielfache von A und B hat den Wert 

A*B/(größter gemeinsamer Teiler). 

Bezeichnet D den größten gemeinsamen Teiler von A und B, so läßt sich D als 
Kombination von A und B schreiben: 

D = S*A + T*B. 

Dabei sind S und T ganze Zahlen. Die Werte von S und T kann man bestimmen, 
indem man den Euklidischen Algorithmus rückwärts anwendet. Ändern Sie das 
Programm zur Bestimmung des größten gemeinsamen Teilers so ab, daß es auch S 
und T berechnet. 
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KAPITEL 9 
PRIMZAHLEN 


Eine Primzahl ist eine ganze Zahl (größer als 1), die durch keine andere ganze Zahl 
außer 1 und sich selbst teilbar ist (unter Teilbarkeit verstehen wir natürlich Teilbar¬ 
keit ohne Rest). Die Zahlen 2, 3, 5, 7 und 11 sind prim, aber 4, 6, 8, 9 und 10 nicht. 
Nicht-Primzahlen werden zusammengesetzt genannt. 

Es gibt keine allgemeine Formel für Primzahlen, aber wie Euklid (etwa 300 v. Chr.) 
gezeigt hat, gibt es unendlich viele Primzahlen. Man weiß auch, daß Primzahlen 
unter großen Zahlen weniger häufig Vorkommen. 

Primzahltests sind von großer Bedeutung. In jüngster Zeit haben sie in der Krypto¬ 
graphie Beachtung gefunden. 

Eine einfache und direkte Methode, zu entscheiden, ob eine Zahl N prim ist, ist das 
sogenannte Sieb des Eratosthenes. (Eratosthenes von Kyrene war ein griechischer 
Mathematiker, ca. 276-196 v. Chr., der auch den Umfang der Erde berechnet hat.) 
Die Idee besteht darin, alle ganzen Zahlen von 1 bis N aufzuschreiben. Lasse dann 
2 stehen und streiche alle geraden Zahlen nach 2. Die nächste nicht durchgestri¬ 
chene Zahl ist prim; das ist 3. Streiche nun jede dritte Zahl nach 3. Die nächste 
verbliebene Zahl nach 3 ist 5. Diese Zahl muß prim sein. Streiche nun jede fünfte 
Zahl nach 5. Setze dieses Verfahren fort. Was übrigbleibt, sind die Primzahlen 
zwischen 1 und N. 

Die folgende Tabelle zeigt das Ergebnis eines Siebs für die Zahlen bis 100. Die 
Vielfachen von 2 sind mit / durchgestrichen, die Vielfachen von 3 mit -, die von 5 
mit \ und die von 7 mit |. 



2 

3 


5 

fr 

7 
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11 

fr 

13 

fr 

fr 

46 

17 

fr 

19 
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fr 

22 

23 

fr 

26 
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fr 

29 

M 

31 

32 

SS 

34 

* 

fr 

37 

38 

SS 

ad 

41 

fr 

43 

44 

4^- 

46 

47 

■fr 

4|9 

sä 

•ST 

52 

53 

fr 

86 

4 

S? 

56 

59 

M 

61 

52 

6fS 

64 


fr 

67 

66 

69 

fr 

71 

fr 

73 

14 

fr. 

76 

fr 

fr 

79 

90 

-84 

82 

83 

fr 

86 

56 

-87 

56 

89 

90 

91 

92 

99 

94 

96 

fr 

97 

96 

99 

t80 


Das nächste Programm verwendet das Sieb des Eratosthenes, um eine Liste von 
Primzahlen unterhalb einer gegebenen Zahl N zu erstellen. Diese Zahl wird zu 
Beginn mit INPUT eingegeben. 
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10 REM SIEB DES ERASTOSTHENES 
20 POKE53230,6:PRINTCHR$< 147 V *** SI 
EB DES ERASTOSTHENES ***" 

30 PR INT:PRINT"DIESES PROGRAMM BERECHNET 
PRIMZAHLEN BISZU EINER EINGEGEBENEN "! 

40 PRINT" ZAHL" 

50 PRINT:INPUT"ß!S ZU WELCHER ZAHL "JN 
S0 IF N< 4 OR NM5000 OR NOINTCM) THEN P 
RINT:PR INT"VERNUENFTIG BLEIBEN":GOTO50 
70 DIM AX( N):PRINTCHRi< 147 M PR INT:PR INT" 
PRIMZAHLEN VON £ EIS'WN 
75 PR INT:PRINTCHR*C 5 > 

80 FORI=£TQN 
90 IFAX<I> = 1 THEN 130 

100 IF 39-P0S< 0)<LEN( STR$< I ) JTHENPRINT 
110 PRINTI, 

120 cor j = I JON STEP I : AX< J > = 1 : NEXT 
130 NEKT:PR INT 

135 PR INT:PR INT:PR INT:PR INT:PR INTCHRf< 15 
4>"NEUER START < J/N) 7" 

150 GET G$: IFG$ = ""THEN150 
160 IFG* = "J"THENRUN 

170 IFG$ = "N"THENPRINTCHRJ< 147 ) "AUF WIEDE 
RSEHEN":END 
180 GOTO 150 

Beachten Sie, daß das Programm für große Zahlen langsam anläuft; nach kurzer 
Zeit beginnt es jedoch, sehr schnell Primzahlen auszudrucken. Es dauert etwa 7 
Minuten und 21 Sekunden, die Primzahlen zwischen 2 und 15 000 auszudrucken. 
Das Sieb des Eratosthenes ist begrifflich einfach und nützlich, wenn man an einer 
Primzahl-Liste interessiert ist. Aber es ist keine sehr praktische Methode, wenn man 
prüfen will, ob eine gegebene Zahl prim ist. 

Man kann auf einfache Weise feststellen, ob eine Zahl N prim ist, indem man 
schrittweise prüft, ob sie durch die Zahlen kleiner als N teilbar ist. Das folgende 
Programm führt diese Methode vor. Es enthält eine Uhr, die anzeigt, wie lange der 
Primzahltest dauert. 

10 REM 24PRIMZAHLEN 1 

20 POKE53280,6!PRINTCHR$< 147)" *** LÄNGS 
AMER PRIMZAHLEN TESTER I ***" 

30 PR INT:PR INT:PR INT:PR INT:INPUT"ZU TEST 
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ENDE ZAHL " ;N 

40 IFN< 4 OR N< >INT< N) THEN PR INT"VERNUEN 

FT IG BLEIBEN"JGOTO30 

45 PR INTCHRS C 147>!PR INTsPR INT 

50 A* = " EINE ":TI$="000000“SX=N-1 

60 FOR I=2 TO X 

70 T=N/I:IF INTC T> =T THEMAS = “ KEINE ": I = 

X 

80 NEXT 

90 PRINTN "IST"A$"PRIMZAHL." 

100 PR INT s PR INT s PR INT“ ZEIT UM ZAHL ZU T 
ESTEN"INT( T1/60+0.5*100>/100"SEKUNDEN" 

110 PR INT:PR INT s PR INT:PR INT" NEUER START 

<J/N> 7" 

120 GETGS:IFG*=""THEN120 
130 IFGS = "J"THENRUN 

140 IFGS = "N"THENPRINTCHRSC 147 > "AUF WIEDE 
RSEHEN"s END 
150 GOTO 120 

Das Programm funktioniert zwar, ist aber ziemlich ineffizient und langsam. Bei¬ 
spielsweise dauert der Test, ob 9001 prim ist, etwa 113 Sekunden. Mit etwas 
Nachdenken läßt sich das enorm beschleunigen. 

Zunächst einmal brauchen wir nicht alle Zahlen zwischen 2 und N zu berücksichti¬ 
gen, sondern nur die Zahlen zwischen 2 und INT(SQR(N)). Denn wenn M eine 
ganze Zahl ist, die N teilt und größer als INT(SQR(N)) ist, dann ist N/M eine ganze 
Zahl, die N teilt und kleiner als INT(SQRfN)) ist. Die nachstehende Version 2 enthält 
diese einfache Ergänzung. 

10 REM 25PRIMZAHLEN £ 

20 POKE53280,6sPRINTCHR$< 147)" *** LÄNGS 
AMER PRIMZAHLEN TESTER II ***“ 

30 PR INT!PR INTs PR INT:PR INT • INPUT"Z U TEST 
ENDE ZAHL *;N 

40 I FN< 4 OR NOINT(N) THEN PR INT "VERNUEN 

FTIG BLEIBEN"SGOTO30 

45 PRINTCHR$< 147)sPR INT:PR INT 

50 A*=" EINE “:TI$="000000":X=SQR<N> 

60 FOR 1=2 TO X 

70 T=N/I:IF INT< T> = T THEMAS = " KEINE ":I = 
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80 NEKT 

90 PRINTN "IST"A$“PRIMZAHL." 

100 PR INT:PR INT * PR INT” ZEIT UM ZAHL ZU T 
ESTEN : "INT< T1/60+0.5*100>/100"SEKUNDEN" 

110 PR INT : PR INT s PR INT:PR INT M NEUER START 

<J/N> ?" 

1£0 GETGf:IFG$=""THEN1S0 
130 IFG*="J"THENRUN 

140 I FG$ = "N " THENPR I NTCHR$( 147 ) "AUF l-JIEDE 
RSEHEN":END 
150 GOTO 1S0 

Diese Version arbeitet wesentlich schneller. Der Primzahltest für 9001 dauert jetzt 
kaum länger als eine Sekunde. Eine größere Zahl, etwa 987654323, wird in 
ungefähr 387 Sekunden getestet. Versuchen Sie nicht, eine so große Zahl mit 
Version 1 zu testen - wenn es Ihnen nicht gerade Spaß macht, Löcher in den 
Bildschirm zu starrren. 

Anmerkung: Die in Version 2 des Primzahltests enthaltene Ergänzung könnte auch 
in das Programm ‘Eratosthenes’ eingebaut werden. Mit dem so erzeugten ‘Eratos- 
thenes 2’ dauert es 6 Minuten und 21 Sekunden, die Primzahlen von 2 bis 15 000 
auszudrucken. 

10 REM 26ERAST0STHENES2 

20 POKE53280,6 * PR INTCHRJ < 147 V *** SI 

EB DES ERASTOSTHENES ***" 

30 PR INT:PR INT"DIESES PROGRAMM BERECHNET 
PRIMZAHLEN BISZU EINER EINGEGEBENEN".: 

40 PRINT" ZAHL“ 

50 PRINT:INPUT"BIS ZU WELCHER ZAHL ";N 
60 IF N<4 OR N>15000 OR N< >INT<N> THEN P 
RINT:PR INT"VERNUENFTIG BLEIBEN":GOTO50 
70 DIM A%< N >:PR INTCHRS< 147 >:PR INT:PR INT" 

PR IMZ AHLEN VON 2 BIS" ! N:S =INT<SQRC N > + 1) 

75 PRINT:PRINTCHR$<5) 

30 FORI=2T0N 

30 IFA ZS, I ) = 1 THEN 130 

100 IF 33-POSC0><LEN<STR$< I>>THENPRINT 
110 PRINTI, 

120 IF I<=S THEN FORJ = I TON STEP I : A%< J ) = 1 
: NE KT 
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130 NEXTSPRINT 

135 PRINT:PRINT:PRINT:PRINT:PRINTCHR*< 15 
4)"NEUER START < J/N> ?" 

150 GET G* ! IFG* = 11 “THEN150 
160 IFG* = “J"THENRUN 

170 IFG* = "N“THENPRINTCHR*( 147 > “AUF WIEDE 
RSEHEN": END 
180 GOTO 150 

Version 2 des Primzahltests läßt sich mit einer Idee aus dem Sieb des Eratosthenes 
weiter verbessern. Wir können alle geraden Zahlen größer als 2 und jede dritte Zahl 
nach 3 weglassen. Diese Änderungen sind in die dritte Version des Programms 
eingebunden worden. 

10 REM 27PRIMZAHLEN 3 

20 POKE53280,6:PR INTCHR*< 147}” *** P 

RIMZAHLEN TESTER III ***" 

30 PR INT:PR INT:PRINTsPRINT! INPUT”ZU TEST 
ENDE ZAHL “;N 

40 IFN< 4 OR MOINT(N) THEN PR INT"VERNUEN 

FTIG BLEIBEN"JGOTO30 

45 PR I NTCHR$( 147 >: PR INT : PR INT 

50 A* = " EINE ":TI$ = "000000":X=SQR< N) 

52 IF INT( N/2)= N/2 THENA* = " KE INE " :GOTO 
30 

53 IF INT( N/3)=N/3 THEMA* = " KE INE ":GOTO 
90 

60 FOR 1=5 TO X STEP 6 

70 T=N/I:IF INT(T)=T THENA*=" KEINE ":I= 

u 

75 T=N/< I +2 > : IF INT<T>=T THENA* = " KEINE 
" : I =X 
80 NEXT 

90 PRINTN "IST"A*"PRIMZAHL." 

100 PR INT:PR INT:PR INT" ZEIT UM ZAHL ZU T 
ESTEN :"INT<TI/60+0.5*100>/100"SEKUNDEN" 

110 PR INT:PR INT:PR INT:PR INT" NEUER START 

< J/N? 7“ 

120 GETG*: IFG* = ""THEN120 
130 IFG*="J"THENRUN 
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140 IFG$ = "N"THENPRINTCHR$( 147 ) "AUF WIEDE 
RSEHEN":END 
150 GOTO 1S0 

Damit läuft das Programm etwas schneller ab. Der Test der Primzahl 987654323 
dauert jetzt etwa 133 Sekunden, was einigermaßen akzeptabel ist. 

Eine zusammengesetzte Zahl läßt sich als Produkt von Primzahlen schreiben. Zum 
Beispiel 

6 = 3*2 

24 = 3 * 2 * 2 * 2 

81018001 = 9001 * 9001 

und so weiter. Die vorkommenden Primzahlen heißen Faktoren der jeweiligen Zahl. 
Mit einigen zusätzlichen Zeilen zu ‘Version 3' kann man alle Faktoren einer Zahl 
ausgeben lassen. 

Das nächste Programm druckt die Faktoren einer Zahl aus. 

10 REM 28PRIMFAKTOREN 

20 P = 0 : GOTO70 : REM SUBROUTINE 

30 Y=N/T: IF38-P0S< 0 >< LEN< STR$< Y> > THEN P 

R INT 

40 N=T s T=N/Y ! I FT= INT( T > THENE =E + 1 : GOTO30 
50 P=P+1:PRINTY"t"E,:S=SQR< N) + li IFP=4THE 
NPRINT:P=0 
55 RETURN 
68 REM HAUPTTEIL 

70 POKE53280,6:PR INTCHR$< 147 >" *** 

PRIMFAKTOREN ***" 

75 PR INT:PR INT:PR INT“DIESES PROGRAMM ZER 
LEGT BELIEBIGE ZAHLENIN PRIMFAKTOREN. 11 
80 PR INT: INPUT"EINGABE ZAHL ";N 
90 IF N<4 OR N< >INT< NJTHEN PR INT:PR INT"V 
ERNUEFTIG BLEIBEN":GOTO80 

100 PR INTCHR$( 147 >:PR INT s PR INT:PR INT"DIE 
PRIMFAKTOREN VON"N"SIND : "CHR$<5):PR INT 


1 10 

T I f = "000000": 

: X=SQR< N>:S = 

X+ 1 


120 

E=1:T=N/2: 

IF 

INT< TI = T 

THEN 

GOSUB30 

130 

E = 1:T = N/3: 

IF 

INT< T) = T 

THEN 

GOSUB30 
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140 FORI=5T0XSTEP6 

150 E=15 T=N/1s IF INT< T ) = T THEN GOSUB30 

1S0 E=1 :T=N/< 1+8) ! IF INT< T) = T THEN GOS 
UB30 

170 IFI>STHENI= X 
180 NEXT 

190 IF 38-POSC0><LENCSTR$<N>) THEN PRINT 

800 IF N>1 THENPRINTN"t 1" 

813 PR INT!PR INT:PR INTsPR INTCHR$ C 154)"ZEI 
T UM FAKTOREN ZU BERECHNEN : " 

880 PR INT:PR INTCHR* ( 5)INT<TI,-'60 * 100 +.5)/ 
100 "SEKUNDEN“CHR$< 154) 

830 PR INT:PR INTSPRINT:PR INT"NEUER START 

< J/'N) 7 11 

840 GETG$: IFG$ = “"THENS40 
850 IFG$ = "J"THENRUN 

860 I FG* = ”N" THENPR INTCHR$'.' 147) "AUF WIEDE 
RSEHEN":END 
870 GOTO840 


Wenn Sie eine sehr große Zahl in das Programm eingeben, rundet der Commodore 
64 allerdings die Zahl und berechnet die Faktoren der gerundeten Zahl. 


GROSSE PRIMZAHLEN 

Es gibt Leute, die zum Zeitvertreib nach großen Primzahlen suchen. Die größten 
bekannten Primzahlen sind gewöhnlich Mersennesche Primzahlen. Zahlen der 
Form 

2 P - 1 

heißen Mersennesche Zahlen nach dem französischen Mönch Marin Mersenne, der 
sich mit dem Problem beschäftigt hat, welche dieser Zahlen prim sind. 

Früher glaubten viele einschlägige Autoren, daß Mersennesche Zahlen prim sind, 
wenn der Exponent P prim ist. Für die ersten paar Exponenten ist das tatsächlich der 
Fall. 
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2 2 - 1 = 3 

2 3 - 1 = 7 
2 S — 1 =31 
2 7 — 1 = 127 

Aber für P = 11 ist die Zahl nicht prim. 

2 " - 1 =2047 = 23 * 89 

Zumindest in den letzten 100 Jahren war die größte bekannte Primzahl der Welt 
immer eine Mersennesche Primzahl (mit Ausnahme einer kurzen Phase im Jahre 
1951). Bis Januar 1983 war die größte bekannte Primzahl 


eine Mersennesche Primzahl, die im April 1979 von David Slowinski entdeckt 
wurde. Es war die 27. bekannte Mersennesche Primzahl. Ungefähr im Januar 1983 
fand Slowinski die viel größere Primzahl 

286243 _ -| 

die eine Mersennesche Primzahl mit 25 962 Dezimalstellen ist. (Vollständig ausge¬ 
schrieben würde diese Zahl mehrere Seiten dieses Buches füllen.) Der Computer, 
den Slowinski dabei benutzte, war ein Cray-1. Obwohl das ein ungeheuer schnelles 
Gerät ist, brauchte der Test 1 Stunde, 36 Minuten und 22 Sekunden Computerzeit. 
Es wäre unmöglich, eine beliebige Zahl in der Größenordnung von 

2^86243 _ -j 

mit einem der Programme aus dem vorigen Abschnitt einem Primzahltest zu 
unterziehen. Es gibt aber besondere Techniken für Mersennesche Zahlen, die in 
den letzten 100 Jahren entwickelt worden sind. Um zu prüfen, ob eine Zahl der 
Form N = 2 P - 1 prim ist, definiert man eine Folge: 

U, = 4 

U 2 = U/U, - N*INT((U 1 *U 1 -2)/N) 

• • 

• • 

• • 

U p —■] = U p -2*Up_ 2 -2- N*INT((U p _2*Up_ 2 -2)/N) 
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Dann ist N prim genau dann, wenn Up., = 0 ist. Das heißt, wir benötigen ungefähr P 
einfache Rechenschritte, um zu testen, ob N prim ist. Für den Commodore 64 
würde es ziemlich einfach, 86 243 Operationen auszuführen, wenn die betreffen¬ 
den Zahlen klein wären. Hier hat man es aber mit Zahlen zu tun, die 86 243 
Binärstelien haben, während der Commodore 64 nur 32 Binärstellen speichert. Man 
kann dieses Problem mit geeigneten Methoden umgehen, aber wir wollen hier nicht 
ins Einzelne gehen. Lesen Sie dazu das nächste Kapitel. 


STOCHASTISCHE PRIMZAHLTESTS 

Bei unbegrenzter Zeit könnten wir jede ganze Zahl daraufhin prüfen, ob sie prim ist, 
indem wir sie zu teilen versuchen. Aber die Zeit ist beschränkt, selbst für einen 
Computer. 

In den letzten Jahren hat man jedoch einen neuen Test auf der Grundlage eines 
alten Satzes von Pierre Fermat entworfen. (Fermat lebte im 17. Jahrhundert.) 
Fermat hat gezeigt, daß B p ~' - 1 durch P teilbar ist, wenn P prim ist und B irgend 
eine andere Zahl zwischen 1 und P - 1 ist. Sei z. B. P = 11 und B = 2. Dann ist 

2 1 ’-' - 1 

gleich 1023, also tatsächlich durch 11 teilbar. 

Fermat hat seinen Satz für alle Werte von B bewiesen; jedoch kannten chinesische 
Mathematiker schon im 5. Jahrhundert v. Chr. den Satz für B = 2. Außerdem 
nahmen sie irrtümlich an, daß die Umkehrung richtig sei. Mit anderen Worten, sie 
waren der Ansicht, daß P prim ist, wenn 2 P 1 durch P teilbar ist. Aber 341 teilt 

2 34 ° _ -| 

ohne daß 341 prim ist. So eine Zahl nennt man pseudoprim zur Basis 2. Allgemein 
heißt eine zusammengesetzte Zahl P, die B p -' - 1 teilt, pseudoprim zur Basis B. Die 
meisten Zahlen, die pseudoprim scheinen, sind in Wahrheit prim. Auf dieser 
Tatsache beruht der Test. 

10 REM 29PRIMZAHLEN 4 

20 POKE53280,6sPRINTCHR$< 147)" *** STOC 

HASTISCHER PRIMZAHLTEST ***« 

30 PR INT SPRINT SPRINT SPRINT! INPUT’*E INGABE 
ZAHL 11 ; N 

40 IFN< 40RN< > INT< N>THENPR INT: PR INT*'VERNU 
ENFTIG BLEIBEN"SGOTO30 
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(2tT>*X 


50 REM FAKTOR N-1 = 

60 T=0 ' X=N-1 
70 D=X,'£:IF INTCD>=D THEN T=T+1:X=D:GOTO 
70 

80 REM BASIS B 

90 B=RNDC -TI > :B= INTC RND( 1 ) *C N-3.000 1 > +8 > 


100 PRINTCHRf< 147 > s PR INT• PR INT s PR INT"TES 
T BA3IS";BCHR$C 5 > 

1S0 P=1 

130 IF H=0 THEN 170 

140 D=X/2! IFD< >INTC D) THENP = P*B!P=P-N#IN 
TC P/N> 

150 B =B *B : B =B -N+ INTC B/N > : X = I NT< D > : GOTO 1 3 
0 

160 REM KONTROLLE BtX 
170 IFP=10RP=N-1THENS00 

180 IFT< 2THENPR I NT '• PR I NT '• PR INTN" IST KEIN 
E PRIMZAHL*:GOTO250 

190 P=P*P:P=P-N*INTC P/N>:T = T-1:GOTO 170 
200 PR INTs PR INTs PR INTN" IST LIAHRSCHL ICH P 
RIMZAHL" 

210 PR INT! PR INT '• PR INTCHR$< 1 54 ) "SOLL NOCH 
MAL MIT EINER ANDEREN BASIS" 

220 PR INT"GETESTET WERDEN <J/N> ?" 

230 GETG-J ; IFG$ < >"J"ANDG$< > "N" THEN230 
240 IFG$ = "J"THEN60 

250 PR I NT! PR I NT s PR I NT ’• PR I NTCHR$< 154VNEU 
ER START CJ/N) ?" 

£60 GETG$! IFG$ = ""THEN260 
270 IFG$ = "J"THENRUN 

280 IFG$ = "N"THENPRINTCHR$< 147 ) "AUF WIEDE 
RSEHEN"!END 
290 GOTO260 

Versuchen Sie es mit Zahlen wie 341 und 561, die beide nicht prim sind. Falls sich 
herausstellt, daß sie pseudoprim zu einer bestimmten Basis sind, nehmen Sie eine 
andere Basis. 

Anmerkung: Das Programm ist schnell, aber auf Grund von Rundungsfehlern 
funktioniert es bei großen Zahlen nicht sehr gut. 
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KAPITEL 10 


VERMISCHTES 


PYTHAGOREISCHE TRIPEL 

Sie erinnern sich bestimmt an den schon erwähnten Satz des Pythagoras. Er 
besagt, daß in einem rechtwinkligen Dreieck mit den Seiten X, Y und Z, wobei Z die 
Hypotenuse (die längste Seite) ist, die folgende Beziehung gilt: 

X 2 + Y 2 = Z 2 

Das klassische Beispiel, das einem dabei sofort einfällt, ist das 3, 4, 5-Dreieck: 

3 2 + 4 2 = 5 2 

Ein anderes Beispiel ist das 5, 12, 13-Dreieck: 

5 2 + 12 2 = 13 2 

Es gibt natürlich unendlich viele verschiedene Beispiele für rechtwinklige Dreiecke. 
Aber wie viele Fälle gibt es, in denen X, Y und Z ganze Zahlen sind? Die Antwort 
lautet: unendlich viele. Solche Zahlen X, Y und Z heißen Pythagoreische Tripel. 
Wie erzeugt man eine Liste Pythagoreischer Tripel? Natürlich sollte eine solche 
Liste keine Tripel enthalten, die Vielfache eines anderen sind. Weil z. B. 3, 4, 5 ein 
Pythagoreisches Tripel ist, ist auch 6, 8, 10 eines. Pythagoreische Tripel, die keinen 
gemeinsamen Faktor haben, heißen primitiv. Also ist 3, 4, 5 ein primitives Pythago¬ 
reisches Tripel, während 6, 8, 10 nicht primitiv ist. 

Es wäre ziemlich zeitraubend, Pythagoreische Tripel zu finden, wenn es nicht 
einigen Mathematikern gelungen wäre, eine elegante Methode zur Erzeugung 
primitiver Pythagoreischer Tripel zu entwickeln. Die Technik läßt sich so be¬ 
schreiben: 

1. Wähle zwei positive ganze Zahlen A und B, so daß gilt: 

(a) A ist größer als B. 

(b) A + B ist ungerade (also ist eine der Zahlen ungerade und die andere 
gerade). 

(c) A und B haben keinen gemeinsamen Teiler außer 1. 
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2. Berechne 


X = A 2 - B 2 
Y = 2*A*B 

Z = A 2 + B 2 

Die Zahlen X, Y, Z bilden ein Pythagoreisches Tripel. Umgekehrt erhält man auf 
diese Weise jedes primitive Pythagoreische Tripel. Mit dieser Technik erzeugt das 
nächste Programm eine Liste Pythagoreischer Tripel. Das Verfahren beginnt bei 
A = 2. Das Programm bestimmt dann alle möglichen Werte von B, die die geforder¬ 
ten Bedingungen erfüllen. Der Wert von A wird dann erhöht und das Verfahren 
wiederholt. Das kann man so lange fortsetzen, wie man will. Man kann auch mit 
einem anderen Anfangswert A beginnen. 

10 REM 30TRIPEL 

£0 POKE53280,6sPRINTCHR$< 147)" *** PY 

THAGORE ISCHE TRIPEL ***" 

30 PR INT:PRINT’DIESES PROGRAMM ZEIGT EIN 
IGE PRIMITIVE PYTHAGOREISCHE TRIPEL" 

40 PR INT!PR INT:PR INT:PR INT:PR INT" 

LEER-TASTE FUER START":K=0:A=2:B=3 
50 GETGf: IFG$< >" "THEN50 

70 PR INTCHRt < 147 > " *** PYTHAGOREISCHE 

TRIPEL +**" 

80 PR INT:PR INT"ZAHL ** X ** **Yü* 

**Y**"CHR#<5):PR INT 
80 B=B-2:IFB<1THENA=A+1:B=A-1 
100 A1=A:B1=B 

110 N=INT<A1/B1):A2=A1-N*B1 

120 IFA2<>0THENA1=B1:B1=A2:G0T0110 

130 IFB1< MTHEN80 

140 K =K + 1 : PRINTK,A*A-B*B ,2*A*B,A*A+B*B 
150 IFINT< K/15X >KX15THEN90 
160 PR INT:PR INT:PRINTCHR^C 154>"WEITERE T 
RIPEL <J/N> ?“ 

170 GETGi: IFG$ = ""THEN170 
ISO IFGt = "J"THEN70 

190 IF6* = "N"THENPOKE53280,14:PR INTCHR$< 1 
47>"AUF WIEDERSEHEN":ENO 
£00 GOTO 170 
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MEHRFACHGENAUE POTENZEN 

Der Commodore 64 verwaltet Zahlen mit 9 signifikanten Stellen genau. Also werden 
ganze Zahlen, die kleiner als 999 999 999 sind, korrekt gespeichert. Größere 
Zahlen werden gerundet. Wie kann man dann große Potenzen von Zahlen genau 
berechnen? Wieviel ist z. B. 2 130 ? Das Problem läßt sich mit mehrfachgenauer 
Arithmetik lösen. 

Mehrfachgenaue Arithmetik ist u. a. möglich, wenn man die Stellen einer Zahl in 
einem Feld M(l) ablegt. Dabei speichert M(0) die letzten 4 Stellen, M(1) die 
vorhergehenden 4, und so weiter. Auf diese Weise würde 987654 gespeichert als 

M(1) = 98, M(0) = 7654. 

Die ursprüngliche Zahl erhält man zurück, indem man die Elemente des Feldes in 
Strings umwandelt und nacheinander ausdruckt. Natürlich kann man sie auch durch 
diese Formel gewinnen: 

M = M(1)*10“ + M(0) 


Angenommen, wir wollen M = 987654 mit N = 23456 multiplizieren. Als erstes 
schreiben wir diese Zahlen in die Felder M(l) und N(l). 

M(1) = 98, M(0) = 7654 
N(1) = 2, N(0) =3456 

Dann bilden wir das folgende Feld: 

C(0) = M(0)*N(0) 

C(1) = M(1)*N(0) + M(0)*N(1) 

C(2) = M(1)*N(1) 

Wir erhalten: 

C(0) = 26452224 
C(1) = 353996 
C(2) = 196 


Daraus berechnen wir das Produkt M*N als 

196* 10 8 + 353996*10“ + 26452224 
= 19600000000 + 35399960000 + 26452224 
= 23166412224 
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Natürlich geschieht die Berechnung nicht wie eben gezeigt - der Commodore 64 
würde die Zahlen einfach runden. Wir streichen vielmehr von links Stellen von C(0), 
bis nur noch 4 übrig sind. Dann addieren wir die gestrichenen Stellen zu C(1). 

C(O) 2224 

C(1) 353996 + 2645 = 356641 

Als nächstes streichen wir von links Stellen von C(1), bis noch 4 übrig sind. Die 
gestrichenen Stellen addieren wir zu C(2). 

C(1) 6641 

C(2) 196 + 35 = 231 

Damit erhalten wir: 

C(2) = 231, C(1) = 6641, C(0) =2224 
Daraus lesen wir sofort das Ergebnis ab: 

987654*23456 = 23166412224 

Dieses Verfahren gilt ganz allgemein. Bei anderen, größeren Zahlen existieren 
vielleicht Werte für die Feldelemente M(2), M(3) etc. In dem Fall würden wir das 
folgende Feld bilden: 

C(0) = M(0)*N(0) 

C(1) = M(1)*N(0) + M(0)*N(1) 

C(2) = M(2)*N(0) + M(1)*N(1) + M(0)*N(2) 
etc. 

Danach wird wie oben nach der Streich-Methode verfahren, bis schließlich das 
Resultat ausgedruckt werden kann. 

Das Programm ‘Mehrfachgenaue Potenzen’ zeigt, wie man Potenzen einer gegebe¬ 
nen Zahl mit Hilfe der eben beschriebenen Technik genau berechnet. Das Pro¬ 
gramm läuft, bis eine 40steilige Zahl erreicht ist. Um die Genauigkeit zu vergrößern, 
brauchen Sie nur den Wert von X in Zeile 80 zu ändern. Der Wert von X + 1 mal 4 ist 
der Genauigkeitsgrad. 
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10 REM 3 1 AUF 40 STELLEN 

20 PQKE53230,6:PRINTCHR$C 147J" *** E 

RHQEHTE STELLENZAHL ***" 

30 PR INT!PR IMT"DIESES PROGRAMM BERECHNET 
POTENZEN VON GANZEN ZAHLEN AUF "f 
40 PR INT"VIERZ IG STELLEN GENAU" 

50 F'RINT; PR INT "BITTE DIE ZU POTENZIEREND 
E ZAHL EINGEBEN" 

60 PRINT:INPUT"ZAHL";N 

70 IFN< 20RN>9993993990RN< >INT( N)THENPRIN 
TPRINT: m ! ! FEHLER ! ! "sGOTO60 
80 T= 10000 : X =3 s Y=X +2 : D I MMC X > , N< X > , LC Y ) : K 
= 1 : NC 0 > =N : MC 0 > =N 
30 FORI=0TO2 

! 00 I F MC I ) > =TTHENQ = I NT< M< I > /T > : MC I > =MC I > 
-Q + T: MC I + 1 > =M( I + 1 ) +Q 
110 NC I > = MC I > 

120 NEXT 

130 PRINTCHR$C 147>"ERHOEHTE STELLENZAHL 
VON";N:PR INTCHRf C5) 

140 K =K + 1 

150 FORI=0 TO Y; L< I > =0 : NEXT 
160 FOR J =0 TO 2:FOR1=0 TO X 
170 LC I + J > =L< I + J ) +NC I > J > 

180 NEXT:NEXT 

190 FOR 1=0 TO X.’NC I >=LC I) : NEXT 
200 FOR 1=0 TO X 

210 IFNC I ) > =T THENQ= INTC NC I)/T>sNC I > =NC I 

> -Q *T s NC 1 + 1 >=NC 1 + 1 >+Q 

220 NEXT 

230 FOR1=0 TO X 

240 IFNC I > >0 THEN L=I 

250 NEXT 

260 PRINTN”t"K 

270 FORI=L TO 0 STEP -1 

280 A$=MID$CSTR$CNCI)>,2> 

230 IFI< LANDLENC A$ X >4 THEN A$ = "0"+A$:GO 
TO290 

300 PRINTA$; ! NEXTSPRINT 

310 IF CMCX>*T+NCX-l>>*N>=T*TTHEN3S0 
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3£0 IFINTC K/10H >K/10THEN140 

330 PR INT!PRINTCHR$< 154)"WEITERMACHEN <J 

y'f' J ) 7 " 

340 GETGf: IFG$< >"N"ANDG*< >"J"THEN340 
350 IFG*="J"THEN130 
355 PRINTCHR$< 147> 

360 PR INTCHR$< 154>:PR INT:PR INT"FEUER STA 
RT ?" 

370 GETG* s IFG*< >"J"ANDG*< >"N"THEN370 
380 IFG$ = "J"THENRUN 

380 PR INT:PRINTCHR$< 147> "AUF WIEDERSEHEN 
M : END 
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KAPITEL 11 
MATRIZEN 


Matrizen sind rechteckige Zahlenfeider, bei denen Position und Wert jeder Zahl von 
Bedeutung sind. Eine Matrix wird gewöhnlich, aber nicht immer, in Klammern 
dargestellt: 

~ 2 —3 6~ 

-4 4 5 


Die Größe einer Matrix ist durch die Anzahl der Zeilen und Spalten gegeben. Eine 
Matrix heißt M-mal-N-Matrix, wenn sie M Zeilen und N Spalten hat. Das obige 
Beispiel ist also eine 2-mal-3-Matrix. Hier sind weitere Beispiele für Matrizen: 


8 12 0 
1 1 8 
8 12 1 


Eine 3-mal-3-Matrix 


9 12 
0 1 
8 -2 


Eine 3-mal-2-Matrix 


9 12 

3 2 


Eine 2-mal-2-Matrix 


Ist die Zeilenzahl einer Matrix gleich der Zahl ihrer Spalten, nennt man sie eine 
quadratische Matrix. 

Matrizen sind nützlich, um Informationen auf übersichtliche Weise zu speichern. 
Nehmen Sie z. B. an, Sie haben den Betrag errechnet, den Ihr Haushalt in jedem 
Quartal eines Jahres für drei Arten von Brennstoff ausgibt. Diese Information 
könnten Sie wie folgt in einer Tabelle darstellen: 


Quartal 

Gas 

Elektrizität 

Fester Brennstoff 

1 

30 

28 

5 

2 

27 

19 

4 

3 

25 

15 

0 

4 

32 

27 

3 


Diese Daten können in eine 4-mal-3-Matrix geschrieben werden. 
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30 

28 

5 

27 

19 

4 

25 

15 

0 

32 

27 

3 


In Ihrem Commodore 64 können Sie Matrizen mit Hilfe zweidimensionaler Felder 
wie etwa A(I,J) speichern. Zum Beispiel könnte die Matrix mit den Brennstoffausga¬ 
ben durch folgende Zuweisungen gespeichert werden: 

A(1,1) = 30 A(1,2) = 28 A(1,3) = 5 
A(2,1) = 27 A(2,2) = 19 A(2,3) = 4 
A(3,1) = 25 A(3,2) = 15 A(3,3) = 0 
A(4,1) = 32 A(4,2) = 27 A(4,3) = 3 

Die Informationen würden natürlich mit einem READ-Befehl eingegeben: 

FOR 1 = 1 TO 4:FOR J = 1 TO 3:READ A(I,J):NEXT:NEXT 
DATA 30,28,5,27,19,4,25,15,0,32,27,3 

Das nächste kurze Programm verdeutlicht, wie Daten eingelesen und dann ausge¬ 
druckt werden. 

5 REM 3 2 D E MO MH TR 12 E N 
10 PR INTCHRT< 147) 

20 FOR I = 1 T04 : FOR J = 1 TOS : REftDFK I , J > : NE.XT • N 
EXT 

30 FQR 1 = 1 TO4 

40 FOR J= 1T03 : PR INTFK I , J > ; : NE.XT 
50 PRINT 
60 NE.XT 

100 DHTH 30,28,5,27,13,4,25,15,0,32,27,3 

Beachten Sie, daß wir die Zähler I und J besser bei 0 statt bei 1 hätten beginnen 
lassen, um den Speicherplatz des Commodore 64 besser auszunutzen, denn 
Felder fangen bei 0 an. Die Zählung von 1 an ist allerdings begrifflich einfacher. 
Ebenso ist zu beachten, daß der Commodore 64 eine DIM-Anweisung erwartet, 
wenn die Anzahl der Zeilen oder Spalten einer Matrix 10 überschreitet. 

ADDIEREN VON MATRIZEN 

Zwei Matrizen können addiert werden, wenn sie gleich groß sind. Anders gesagt, 
zwei Matrizen können addiert werden, wenn sie die gleiche Anzahl von Zeilen und 
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die gleiche Anzahl von Spalten haben. Diese Addition wird nacheinander für je 2 
Elemente in den gleichen Positionen ausgeführt. Zum Beispiel 


2 0 5 

3 2 1 

1+2 5+0 4+5 

-2 + 3 3 + 2 1 + 1_ 

3 5 9 
1 5 2 


WOZU MATRIZEN ADDIEREN? 

Weiter oben hatten wir ein Beispiel mit Brennstoffkosten für vier Quartale eines 
Jahres. Nehmen wir an, die entsprechende Matrix sieht im nächsten Jahr so aus: 


34 

27 

4 

30 

20 

3 

20 

16 

1 

35 

29 

5 


1 

5 

4 

+ 

-2 

3 

1 



Wenn wir die beiden Matrizen addieren, erhalten wir die Ausgaben, die wir für die 
drei Brennstoffe in einem gegebenen Quartal über zwei Jahre gehabt haben. Das 
Ergebnis ist diese Matrix: 


64 

55 

9 

57 

39 

7 

45 

31 

1 

67 

56 

8 


Ihr Commodore 64 kann Matrizen oder zweidimensionale Felder sehr leicht addie¬ 
ren. Seien zwei Felder A(I,J) und B(I,J) gegeben, wobei jedesmal I von 1 bis 4 und J 
von 1 bis 3 läuft. Durch Addition der beiden Felder ergibt sich ein drittes Feld wie 
folgt: 

FOR 1 = 1 TO 4:FOR J = 1 TO 3 
C(I,J) = A(I,J) + B(I,J) 

NEXT :NEXT 

Das Subtrahieren von Matrizen geschieht nach den gleichen Regeln wie das 
Addieren. 
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MULTIPLIZIEREN VON MATRIZEN 


Das Multiplizieren von Matrizen ist begrifflich nicht ganz einfach. Um zwei Matrizen 
multiplizieren zu können, muß die Anzahl der Spalten in der ersten Matrix mit der 
Anzahl der Zeilen in der zweiten Matrix übereinstimmen. Man kann also eine 4-mal- 
3-Matrix mit einer 3-mal-2-Matrix multiplizieren, aber nicht mit einer 2-mal-3-Matrix. 
Das Produkt einer M-mal-N-Matrix mit einer N-ma!-P-Matrix ist eine M mal P-Matrix. 
Wie das Verfahren tatsächlich abläuft, läßt sich am besten an einem Beispiel zeigen, 
das anschließend erläutert wird. Angenommen, wir wollen die folgenden Matrizen 
multiplizieren: 


2-3 6 

-4 4 5 


Die erste Matrix, eine 2-mal-3-Matrix 


9 12 

0 1 

8 -2 


Die zweite Matrix, eine 3-mal-2-Matrix 


Da die Anzahl der Spalten der ersten Matrix mit der Anzahl der Zeilen der zweiten 
Matrix übereinstimmt, kann die erste mit der zweiten multipliziert werden. 


2 

-4 


-3 

4 


1 

9 

12 ’ 

* 

0 

1 

- 

8 

-2 


2*9 + -3*0 + 6*8 
-4*9 + 4*0 + 5*8 

66 9 ' 

4 -52 


2*12 + -3*1 +6*-2 
-4*12 + 4*1 + 4*-2 


Zum Multiplizieren sehen Sie sich zunächst jede Zeile der ersten Matrix und jede 
Spalte der zweiten Matrix an. Alle haben die gleiche Anzahl von Elementen. Bei 
jeder dieser Zeilen und Spalten multipliziert man zunächst jeweils die ersten 
Elemente, dann die zweiten und so fort, und zum Schluß addiert man alle so 
gewonnenen Produkte zusammen. Auf diese Weise erhält man ein Element der 
Produktmatrix. 
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Auf Ihrem Commodore 64 können Sie das Produkt der M-mal-N-Matrix A(I,J) mit 
der N-mal-P-Matrix B(I,J) wie folgt bilden. Das Ergebnis ist eine M mal P-Matrix 
C(U). 

FOR 1 = 1 TO M 
FOR J= 1 TO P 
C(I,J) = 0 

FOR K=1 TO N : C(I,J) = C(I,J) + A(I,K)*B(K,J) : NEXT 

NEXT 

NEXT 

Wenn man das Produkt A*B zweier Matrizen A und B bilden kann, braucht man 
nicht auch das Produkt B*A bilden zu können. Selbst wenn das Produkt B*A 
existiert, muß es nicht gleich A*B sein. Suchen Sie nach einfachen Beispielen! 


WOZU MATRIZEN MULTIPLIZIEREN? 

Wir wollen uns noch einmal das Beispiel mit den Brennstoffkosten ansehen. 
Vergegenwärtigen Sie sich die einzelnen Daten: 


Quartal 

Gas 

Elektrizität 

Fester Brennstoff 

1 

30 

28 

5 

2 

27 

19 

4 

3 

25 

15 

0 

4 

32 

27 

3 


Angenommen, es liegen zwei Schätzungen für den mutmaßlichen Kostenanstieg 
für diese Brennstoffe vor: 



Schätzung 1 

Schätzung 2 

Gas 

10% 

5% 

Elektrizität 

5% 

10% 

Fester Brennstoff 

10% 

10% 


Die nächste Tabelle zeigt diesen Anstieg in Dezimaldarstellung. 



Schätzung 1 

Schätzung 2 

Gas 

0.1 

0.05 

Elektrizität 

0.05 

0.1 

Fester Brennstoff 

0.1 

0.1 


in 



Wieviel mehr würden Sie in jedem Quartal für die drei Brennstoffe zahlen müssen? 
Das hängt davon ab, welche Schätzung Sie zugrundelegen. Das Ergebnis läßt sich 
wie folgt tabellieren: 


Quartal 


Schätzung 

1 


Schätzung 2 



1 

30* 

.1 

+ 

28*.05 

+ 

5*.1 

30*.05 

+ 28*.1 

+ 

5* 

.1 

2 

27* 

.1 

+ 

19*05 

+ 

4M 

27*.05 

+ 19*.1 

+ 

4* 

.1 

3 

25* 

.1 

+ 

15*.05 

+ 

0*.1 

25*.05 

+ 15*.1 

+ 

0* 

.1 

4 

32* 

.1 

+ 

27*.05 

+ 

3*,1 

32*.05 

+ 27*.1 

+ 

3* 

.1 


Nach dem Ausrechnen ergibt sich diese Tabelle: 


Quartal 

Schätzung 1 

Schätzung 2 

1 

4.9 

4.8 

2 

4.05 

3.65 

3 

3.24 

2.75 

4 

4.85 

4.6 


Wie Sie vielleicht bemerkt haben, ist die so erhaltene Matrix das Produkt der Matrix 
der Ausgaben mit der Matrix der Schätzungen. Das heißt, es ist das nachstehende 
Matrizenprodukt: 


30 

28 

5 


0.1 

0.05 

27 

19 

4 

* 

0.05 

0.1 

25 

15 

0 


0.1 

0.1 

32 

27 

3 





NULLEN UND EINSEN 

Matrizen, die nur aus Nullen bestehen, heißen Nullmatrizen. Multipliziert man eine 
Matrix mit einer Nullmatrix, so erhält man eine Nullmatrix. Zum Beispiel: 


"30 

28* 


0 

0~ 

27 

19 

* 

0 

0_ 

25 

15 





30*0 + 28*0 
27*0 + 19*0 
25*0 + 15*0 


30*0 + 28*0 
27*0 + 19*0 
25*0 + 15*0 
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0 0 

= 00 
0 0 

Das Produkt zweier Nicht-Nullmatrizen kann eine Nullmatrix ergeben, wie das 
nächste Beispiel zeigt. 

1 fl fl -l" 

1 1j -11 

1 * 1 + 1*-1 
_ri+i*-i 

= o°" 

.o 0 

Eine quadratische Matrix mit Einsen auf der Diagonalen (von links oben nach rechts 
unten) und sonst nur Nullen heißt Einheitsmatrix. Diese beiden Matrizen z. B. sind 
Einheitsmatrizen: 

1 0 0 

0 1 0 

0 0 1 


1 0 
0 1 


r-i + ri 
i*-i + ri 


Multiplikation einer Matrix mit einer Einheitsmatrix läßt die Matrix unverändert. In 
anderen Worten verhält sich die Einheitsmatrix ganz ähnlich wie die 1 bei der 
üblichen Multiplikation von Zahlen. Zum Beispiel: 


30 

28 


1 

o" 

27 

19 


_0 

1 

25 

15 





30*1 + 28*o" 


'30*0 + 28* l" 

27*1 + 19*0 


27*0 + 19*1 

25*1 + 15*0 


25*0+ 15*1 


30 

28 

27 

19 

25 

15 
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INVERTIEREN VON MATRIZEN 


Die Inverse oder Umgekehrte einer Matrix, sofern sie existiert, hat dieselbe Eigen¬ 
schaft wie das Inverse einer gewöhnlichen Zahl. Das Inverse der Zahl 4 ist 0.25, und 
es gilt 

4 * 0.25 = 1, 0.25*4 = 1 

Für eine Matrix A wird die Inverse (wenn sie existiert) mit A 1 bezeichnet, und es gilt 
A * A' 1 =1, A* 1 * A = 1 

wobei I für eine Einheitsmatrix steht. 

Nur quadratische Matrizen können Inverse besitzen. Tatsächlich haben nur 
bestimmte quadratische Matrizen eine Inverse. 

Es gibt verschiedene Methoden, die Inverse einer quadratischen Matrix zu finden. 
Für Ihren Commodore 64 ist das unten beschriebene schrittweise Verfahren gut 
geeignet. Es besteht darin, sogenannte Zeilenoperationen auszuführen. Vereinfa¬ 
chend könnte man sagen, wir setzen zunächst eine Einheitsmatrix neben die Matrix, 
die wir invertieren wollen. Die Zeilenoperationen werden dann auf beide Matrizen 
gleichzeitig angewandt, bis die ursprüngliche in eine Einheitsmatrix umgewandelt 
ist. Die Matrix, die zu Beginn eine Einheitsmatrix war, ist jetzt die Inverse. 

Ein Programm zur Berechnung der Inversen von Matrizen geben wir später an. 
Vorher beschreiben wir kurz die Methode. 

Angenommen, wir wollen die Inverse der folgenden Matrix finden: 

2 1 
— 2 4_ 

Wir setzen eine 2-mal-2-Einheitsmatrix neben diese Matrix: 


1 

C\l 

1_ 


fl 0] 

L-2 t 


1 O 

1_ 


Da die linke obere Eintragung eine 1 werden soll, teilen wir die ganze erste Zeile 
durch 2: 


r 

LO 

ö 

i _ 


0.5 Ol 

-2 4 


—k 

1 _ 


Als nächstes soll die linke untere Eintragung 0 werden; also addieren wir das 
Zweifache der ersten Zeile zur zweiten Zeile: 
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’l 0.5~ 


"0.5 0" 

0 5 


1 1 


Nun soll aus der rechten unteren Eintragung (der linken Matrix) eine 1 werden. Dazu 
teilen wir die zweite Zeile durch 5: 


1 0.5 


" 0.5 0 

0 1 


0.2 0.2 


Um aus der linken Matrix eine Einheitsmatrix zu machen, subtrahieren wir das 
Einhalbfache der zweiten Zeile von der ersten: 


1 


L° ij 



0.4 

0.2 


- 0.1 

0.2 


Die rechte Matrix ist jetzt die Inverse der Ausgangsmatrix. Folgende Rechnung 
bestätigt das: 


2 1 


f0.4 

-0.1 

-2 4 

* 

1 0.2 

0.2 


2*.4 + 1*.2 2 * —.1 + 1 *.2 

—2*.4 + 4*.2 —2*-.1 + 4*.2 


0 

1 


Das unten angeführte Programm berechnet die Inverse quadratischer Matrizen, 
sofern überhaupt eine Inverse existiert. 

10 REM 33 INVERSE MATRIZEN 

20 PÜKE53280,S:PR INTCHR$< 147)" *** 

INVERSE MATRIZEN ***" 

30 PR INT!PR INT“DIESES PROGRAMM BERECHNET 
DIE INVERSE EINER N MAL N MATRIX" 

40 PR INT:PR INT"BITTE DIE GROESSE DER MAT 
RIX EINGEBEN" 

50 PR INT: INPUT"GROESSE VON N " ; N 

60 IFN< 1ORN< >INT< N > THENPRINT"! ! UNMOEGLI 

CH !!":GOTO50 

70 REM EINGABE DER MATRIX 


115 



20 D IM A(N,N),B( N, N > 

30 FOR1 = 1TON!FORJ =1 TON 

100 PR INTCHR$( 147 > s PR INT! PR INT"E INGABE D 
ER MATRIX < TERM FUER TERM>“ 

110 PR INTCHR$ C 13>:PR INT:PR INTsPR INT:PR IN 
T"R£IHE"; I , 

130 PR I NT" SPALTE " } J f s INPUT A<I,J) 

140 NEXTJNEXT 
150 REM RECHNUNG 

160 PRINTCHR*< 147)"DIE MATRIX"iPR INT 
165 P =7: IFN>STHENP=INT< 40/<N+1 > > 

170 FORI = 1TON!O=0sFORJ=lTON:O =0 +P s PRINTA 
< I,J>TAB<0>sNEXTSPRINTsNEXT 
180 PRINT:PR INT"RECHNUNG" 

190 FOR1 = 1 TON:B< I,I> = 1sNEXT 
200 X =0:GOSUB300 
210 REM SCHLUSS 

220 PR INT"NEUER START ( J/N) 7" 

230 GETG*: IFG*< >”N"ANOG$< >"J"THEN230 
240 IFG$ = "J"THENRUN 

250 PRINTCHR$<147>"AUF WIEDERSEHEN":END 

300 REM RECHNUNG 

310 X=X+1:Z=X:PRINT"*"; 

320 IF A< Z , X )=0 THEN Z=Z + 1:IFZ<=N THEN32 
0 

330 IFZ>NTHENPRINTSPRINTSPRINT"KE INE INV 
ERSE MATRIX"sRETURN 

340 IFZ< >XTHENR=1/A( 2,X)s I=X:K*ZsGOSUB50 
0 

350 I FA<X,X X >1THENR =1/A< X,X > s I =X:G0SUB5 
50 

360 FOR1=1TON 

370 IFI=XTHEN1=1+1sIFI>NTHEN390 

380 IF A< I , X X >0THENR = -A< I , X X K =X s GOSUB50 

0 

390 NEXT 

400 IFX< NTHEN310 

4 10 PR INT:PRINT"DIE INVERSE MATRIX"CHR*< 

5 >:PR INT 

4 15 P =7 s IFN>6THENP = I NT< 40X< N+ 1 > > 
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4H0 FOR1 = 1TONs 0=0!FORJ =1TONs 0=0 +P s PR INTB 
< I,J)TAB< O >:NEXT:PR I NT:NEXT 
430 PR I NTCHR-f < 154 >s RETURN 
5Q0 REM RECHNUNG 
5 10 FORJ =1TON 

5£0 I , J ) = RK I , J > +R*A< K , J > • BC I , J ) =B< I , J > 

+R*B<K,J> 

530 NEXT 
540 RETURN 
550 REM RECHNUNG 
560 FORJ=1T0N 

570 f\<. I ,J>=R*A< I ,J> :B< I , J>=R*B< I ,J) 

580 NEXT 
590 RETURN 


GLEICHUNGSSYSTEME 

Matrizen lassen sich zum Lösen von Gleichungssystemen verwenden. Das nächste 
Beispiel zeigt ein System zweier Gleichungen mit zwei Unbekannten: 

3*X + 1*Y = 7 
5*X + 2*Y = 9 

Dasselbe sieht in Matrixschreibweise so aus: 


3 l" 


V 


7 ~ 

5 2 


Y 


_9 


Eine mögliche Lösung der Aufgabe besteht darin, die Inverse der 2-mal-2-Matrix 
auf der linken Seite zu finden und dann die Gleichung mit dieser Inversen zu 
multiplizieren. Im Beispiel ist die Inverse die Matrix 

2 -1 
-5 3 
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Wir erhalten also die Gleichungskette 

X 
Y 

= 1 o 

0 1 _ 

2 -l" 

-5 3 

= 2 ~l" 

-5 3 

= 2*7-1*9 

—5*7+3*9 

5 

-8 

Also lösen X = 5 und Y = -8 das Gleichungssystem. 

Das unten abgedruckte Programm ahmt im wesentlichen diese Schritte nach, um 
ein System von N Gleichungen in N Unbekannten zu lösen. Zu Anfang des 
Verfahrens wird wie oben das Gleichungssystem in zwei Matrizen A(I,J) und B(l) 
geschrieben. Das obige Beispiel würde wie folgt dagestellt: 




Dann werden Zeilenoperationen ausgeführt, bis die Matrix auf der linken Seite eine 
Einheitsmatrix geworden ist. Die rechte Matrix geht dabei in die Lösung des 
Gleichungssystems über. 

10 REM 34GLEICHUNGSSYSTEME 

20 POKE53280,6sPRINTCHRS(14?)" *** G 

LE ICHUNGSSYSTEME *#*" 

30 PR INT:PR INT"DIESES PROGRAMM LOESST GL 
EICHUNGSSYSTEMEMIT N UNBEKANNTEN"; 

40 PR INT:PR INT"ANZAHL DER UNBEKANNTEN" 

50 PRINT:INPUT "ZAHLEN 

60 IFN<20RN< >INTC N>THENPRINT:PRINT"NOCHM 
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AL VERSUCHEN":GOTO50 
70 REM EINGABE KOEFFIZENTEN 
80 DIM ACN,N),B<N) 

90 FOR I = 1 TON!FORJ = 1TON 

100 PRINTCHRfC 147)"EINGABE MATRIX DER KO 
EFF 12 ENTEN C TERM FUER TERM)" 

1 10 PR I NT s PR I NT 11 RE I HE " ) I , "SPALTE")J) : INP 
UTAC I , J> 

120 NEKT 

130 PR INT: INPUT"RECHTE SEITE")BC I> 

140 NEKT 

180 REM RECHNUNG 

190 PRINTCHR$C 147)"MATRIX DER KOEFFI2ENT 
EN":PRINT 

195 P=INTC 40/CN+l .3) > 

£00 FORI = 1TON!FORJ =1 TON:PRINTAC I,J)TAB<P 
* J ) : NEXT: PR INTTABC P *C N+ . 5 ) ) BC I ) : NEXT 
£ 10 P R I NT: P R I NT " R E C H NU NG " 

££0 X=0:GOSUB300 
£30 REM SCHLUSS 

£40 PR INT:PRINT"NEUER START C J/N) ?" 

£50 GETGJ: IFG$< >"J"ANDG$< >"N"THEN£50 
£60 I FG$ = "J"THENRUN 

£70 PRINTCHRfC 147) "AUF WIEDERSEHEN":END 

300 REM RECHNUNG 

310 X=X +1:Z=X:PRINT"*" ) 

3£0 IFAC2,X ) =0THENZ=2 + 1: IF2< =NTHEN3£0 
330 IF2>NTHENPRINTSPRINT:PR INT"KE INE LOE 
SUNG":RETURN 

340 IF2 < >XTHENR=1/A<Z,X): I=X:K=2:GOSUB50 
0 

350 IFAC X , X X >1 THE NR = 1 / A< X X ) : I =X : G0SUB5 
50 

3S0 FOR1 = 1 TON 

370 IFI=XTHENI=I+1:IFI>NTHEN390 
380 IFAC I ,X)< >0THENR = -AC I ,X):K=X:GOSUB50 
0 

390 NEXT 

400 IFX< NTHEN310 

4 10 PR INT:PR INT"DIE LOESUNG":PR INT 



420 FOR1 = 1TON 

430 IF POS< 1)+LEN<STR$< B< I)> >>38THENPRIN 
T 

440 PRINTB< I> 

450 NEMT 

460 PR INT:RETURN 

500 : 

5 10 FORJ =1 TON 

520 fl< I,J)=AC I,J)+R*WK,J) 

530 NEMT 

540 B< I>=B< I >+R*B<K>:RETURN 
550 : 

560 FORJ =1 TON 

570 (=K I ,J>=R*A< I ,J> 

580 NEMT 

590 BC I > = R*B< I ) : RETURN 
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KAPITEL 12 


CODES 

Der Bedarf an gesicherter Datenübertragung wird ständig größer. Dadurch wurde 
die Forschung auf dem Gebiet der Kryptographie - der Kunst des Verschlüsselns 
- sehr angeregt. Kryptographie besteht aus zwei Schritten: die Nachricht bzw. die 
Daten müssen zuerst verschlüsselt und dann wieder entschlüsselt werden. 

TAUSCHALPHABETE 

Die einfachsten Schlüssel sind diejenigen, in denen jeder Buchstabe durch etwas 
anderes ersetzt wird, in der Regel durch einen Buchstaben desselben Alphabets. 
Hier ist ein Beispiel: 

verschlüsselnd verschlüsselnd 
•«—entschlüsseln «^-entschlüsseln 
A - A B - F 

C - K D - P 

E - U F - Z 

G - E H - J 

1-0 J - T 

K - Y L - D 

M - I N - N 

O - S P - X 

Q - C R - H 

S - M T - R 

U - W V - B 

W - G X - L 

Y - Q Z - V 

Die Nachricht MESSAGE würde in IUMMAEU übersetzt. 

In diesem Schlüssel ist das Tauschalphabet nicht zufällig gewählt worden. Werden 
die Buchstaben des Alphabets von 0 bis 25 durchnumeriert, dann wird der Buch¬ 
stabe mit der Nummer NUM durch den Buchstaben mit folgender Nummer ersetzt: 

NUM*5 - 26*INT(NUM*5/26) 

Diese Zuordnung wird rückgängig gemacht durch die Formel 

NUM*21 - 26*INT(NUM*21/26) 
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Dieses Tauschalphabet kann mit einem kurzen Programm aufgelistet werden: 


10 REM 35EINFACHCODE 

£0 PRINTCHR*< 147 >"CODIERUNG"SPRINT 

30 FORI=0TO£5 

40 J =I*5-26 * INT( I*5/26 > 

50 PRINTCHR*<65+I >" -"CHR*< 65+J ) , 

60 NEXT 

70 PR INT!PR INT!PR INT"ENTSCHLUESSELUNG":P 
RINT 

80 FORI= 0TO25 

100 J=I*£1-26*INT< I*£l/£6> 

110 PRINT CHR$< 65+IV -"CHR$( 65+J >, 

120 NEXT 

Um einen Buchstaben zu verschlüsseln, multiplizieren wir die Nummer des Buch¬ 
stabens mit 5 und lassen Vielfache von 26 weg. Zum Entschlüsseln des Buchstaben 
dividieren wir durch 5 und lassen Vielfache von 26 weg. Ebenso könnten wir mit 21 
multiplizieren und Vielfache von 26 vernachlässigen, weil gilt: 

5 * 21 = 105 

= 1 + 4 * 26 

Mit anderen Worten: Wenn wir Vielfache von 26 weglassen, ist 21 reziprok zu 5. 
Man sagt auch, daß 21 das Inverse von 5 modulo 26 ist. 

Dieser Schlüssel hat den offensichtlichen Nachteil, daß eine Nachricht wie PLEASE 
COME QUICKLY zu XDUAMU KSIU CWIKYDQ verschlüsselt würde: Leerstellen 
bleiben Leerstellen. Um diesem Mangel abzuhelfen, sollten wir Leerstellen, Punkt, 
Komma, Ziffern und vielleicht das Fragezeichen in unser Tauschalphabet mit 
aufnehmen. 

Da uns ein Computer zur Verfügung steht, sollten wir die ASCII-Zeichen verwen¬ 
den. Es bieten sich die 59 ASCII-Zeichen von 31 bis 90 an; sie umfassen alle 
gewünschten Buchstaben, und außerdem ist 59 eine Primzahl, was sich später 
noch als nützlich erweisen wird. CHR$(31) ist ein Farbcode und wird nicht wirklich 
gebraucht. 

Verschlüsselt wird in erster Linie durch Multiplikation mit 5. Der Ausdruck für die 
Verschlüsselung des einzelnen Zeichens A$ zu C$ lautet präzise: 

N = 5 * (ASC(A$)-31) 

C$ + CHR$(31 + N - 59*INT(N/59)) 

Entschlüsselt wird durch Multiplikation mit 12 (das Produkt von 5 und 12 ist 60, also 
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gleich 1 ohne Vielfache von 59). Um ein einzelnes Zeichen C$ zu entschlüsseln, 
gehen wir demnach wie folgt vor: 

M = 12 * (ASC(C$)-31) 

A$ = CHR$(31 + M - 59*INT(M/59)) 

Hier ist ein kurzes Programm, das auf dem genannten Schlüssel beruht und eine 
Nachricht verschlüsselt oder entschlüsselt. Als zusätzliche Besonderheit müssen 
Sie eine der Zahlen von 2 bis 58 zum Verschlüsseln und Entschlüsseln eingeben. 
Natürlich muß beim Entschlüsseln dieselbe Zahl wie beim Verschlüsseln gewählt 
werden. Zum Verschlüsseln multiplizieren wir im wesentlichen mit der gewählten 
Zahl, sagen wir N. Zum Entschlüsseln multiplizieren wir mit dem Inversen von N 
modulo 59, also mit einer Zahl M, für die N*M bis auf Vielfache von 59 gleich 1 ist. 

10 REM 36TAUSCHC0DE 

£0 POKE53280,6sPRINTCHR$( 147)" *** 

TAUSCHCODE ***" 

30 PR INT!PRINT"DIESES PROGRAMM VER- UND 
ENTSCHLUESSELT NACHRICHTEN" 

40 PR INT!PR INT"BITTE CODE SCHLUESSEL EIN 
GEBEN" 

50 PR INT:INPUT"ZAHL ZWISCHEN 2 UND 58 
C 

60 IFC<MNT<O0RC<20RO58THEN PRINT:PRIN 
T"FALSCHE EINGABE":GOTO50 
70 PR INT:PR INT"SOLL VER- ODER ENTSCHLUES 
SELT WERDEN 7" 

80 PR INT: INPUT"V ODER E " } A$ 

90 IFA$ < >"E"ANDA* < >"V"THEN80 

100 D = 1 : I =C : B$ = "ENT" : IFA$ = 11 V "THEND=C : B$ = 

"VER":GOTO 130 

1 10 C=C+I : D =D+ 1 : IFC-59*INT< C/59X MTHEN1 
10 

130 PRINTCHR*< 147>"NACHRICHT EINGEBEN <B 

IS ZU 254 ZEICHEN)":PR INT 

140 M$="":N$="":L=0 

150 GETG$:IFG*=""THEN150 

170 G=ASC< G $) 

180 IFG>31ANDG<91THENM*=M$+G$:PRINTG$ ; :L 
=L +1 
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130 IFG=20ANDL>0THENPRINTG$; :L=L-1 : M$ =LE 
FT#<M$,L > 

200 IFG< >13ANDL<255THEN150 
210 PRINT" " 

220 REM VER- ENTSCHLUESSELN 
230 FOR I = 1TOL 

240 N=ASC<MID#<M$,I,1>>-31:N=D*N 
250 N#=N# +CHR#< 31+N-59*INT<N/59>) 

260 NEKT 

270 PRINTCHR#< 147> "DIE "B$"SCHLUESSELTE 
NACHRICHT LAUTET":PR INT:PRINTN$ 

280 PR INT! PR INT: PR INT;PR INT"NEUER START 

<J/N> 7" 

230 GETG#: IFG#< >"J"ANDG$< >"N"THEN£90 
300 IFG# = " J"THENRUN 

310 PRINTCHR#( 147> "AUF WIEDERSEHEN":END 

All unseren Bemühungen zum Trotz ist es (für Experten) leicht, unseren Code zu 
knacken, d. h. zu entschlüsseln. Das Problem beim Tauschalphabet ist, daß jedes 
Zeichen prinzipiell durch ein festgelegtes anderes Zeichen dargestellt wird. Manche 
Buchstaben und Buchstabenpaare kommen jedoch häufiger vor als andere. Bei¬ 
spielsweise macht der Buchstabe E in einem normalen englischen Text ungefähr 
13% der Zeichen aus, T etwa 9%, P etwa 2% und Q etwa 0.2%. Mit solchen 
Informationen kann man ein Tauschalphabet dechiffrieren, obwohl man mit 59 
Zeichen (ASCII-Codes 31 bis 90) ungefähr 1.4*10 80 Tauschalphabete bilden kann. 
Daran zeigt sich, wie sehr die große Fülle von Möglichkeiten täuschen kann. 

Im nächsten Abschnitt werden Codes vorgestellt, die ein gegebenes Zeichen nicht 
immer durch dasselbe Zeichen ersetzen. 


MATRIXSCHLÜSSEL 

Zum Verschlüsseln von Nachrichten kann man Matrizen verwenden. Dieses Verfah¬ 
ren erläutern wir an einem Beispiel: Angenommen, wir wollen die Nachricht 
PLEASE COME QUICKLY verschlüsseln. Dazu ordnen wir sie zunächst in zwei 
Zeilen an: 

PES OEQIKY 
LAECM UC L . 
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Daraus bilden wir dann eine zweizeilige Matrix, indem wir die Buchstaben in ihre 
ASCII-Codes minus 31 verwandeln: 

49 38 52 1 48 38 50 42 44 58~ 

45 34 38 36 46 1 54 36 45 15 

Nun multiplizieren wir die Matrix von links mit der folgenden Matrix: 

2 - 1 ~ 

-5 3_ 

Das Ergebnis ist: 


2 

-1| 


49 38 

52 

1 48 38 50 

42 

44 58 

5 

3 


45 34 

38 

36 46 

1 54 

36 

45 15 

- 

53 

42 

66 

-34 

10 

75 

46 

48 


-110 

-88 

-146 

103 

-102 -187 - 

88 

-102 - 


Dann wandeln wir die Elemente der Matrix in Zahlen zwischen 0 und 58 um, indem 
wir Vielfache von 59 addieren oder subtrahieren, und zwar nach der Regel: 

N = N - 59*INT(N/59) 

Es ergibt sich diese Matrix: 

53 42 7 25 50 16 46 48 43 42~ 

8 30 31 44 16 49 30 16 33 50 

Schließlich addieren wir 31 zu allen Zahlen und fragen mit PRINT CHR$(X) nach den 
zugehörigen Zeichen. Das Ergebnis sieht so aus: 

TI & 8 Q / M 0 J I 

’ = > K / P = / @ Q 

Die endgültige verschlüsselte Botschaft ist 

TT=&>8KQ//PM=0/J@IQ 
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Beachten Sie, daß in der ursprünglichen Form der Buchstabe E dreimal vorkommt. 
Diese drei E sind zu I, > und / verschlüsselt. Der Schlüssel ist also subtiler als ein 
einfaches Tauschalphabet. 

Ein Code taugt nichts, wenn Nachrichten nicht wieder entschlüsselt werden kön¬ 
nen. Der Trick bei dem gerade beschriebenen Code besteht darin, das Verfahren 
umzukehren. Auf den ersten Blick sieht das schwierig aus, aber wir können etwas 
Mathematik zu Hilfe nehmen. 

Wir wollen folgende Nachricht entziffern: 

l = M?"S?$:8M@>’<M 

Zunächst schreiben wir die verschlüsselte Nachricht in zwei Zeilen: 


IM” ? M > < 
= ? S $ 8 @ ’ M 


Dann bestimmen wir die ASII-Codes minus 31, um eine Matrix zu erzeugen: 

42 46 3 32 27 46 31 29 

30 32 52 25 5 33 8 46_ 

Diese multiplizieren wir von links mit der Matrix 


3 1 

_5 2 _ 

Beachten Sie, daß diese Matrix nicht dieselbe ist, die zum Verschlüsseln benutzt 
wurde. Die Ehtschlüsselungsmatrix ist vielmehr die Inverse der Verschlüsselungs¬ 
matrix, wie die nachstehende Rechnung zeigt. 


3 

5 


1 

2 


2 

-5 


-1 

3 


3*2 + 1*—5 
5*2 + 2*-5 


3* —1 + 1*3 
5*—1 + 2*3 


1 0 
0 1 
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Wir multiplizieren also die Entschlüsselungsmatrix mit der Matrix, in der die Nach¬ 
richt verschlüsselt ist: 


3 

1 


42 

46 

3 

32 

27 

46 

31 

_5 

2_ 


30 

32 

52 

5 

25 

33 

8 


r 156 

170 

61 

101 

106 

171 

101 

133] 



[270 294 119 170 185 296 171 237j 

Dann addieren (oder subtrahieren) wir Vielfache von 59, um Zahlen zwischen 0 und 
58 zu erhalten: 

"38 52 2 42 47 53 42 15* 

34 58 1 52 8 1 53 1 

Schließlich addieren wir 31 zu diesen Zahlen und suchen die zugehörigen Zeichen: 

ES! I N T I 

A Y S ’ T 

Auf diese Weise erhalten wir die Mitteilung EASY! ISN’T IT, 

Das nächste Programm verwendet solche Matrizen, um Nachrichten zu verschlüs¬ 
seln und/oder zu entschlüsseln. 

10 REM 37MATRIXCODE 

20 PQKE53280 ,6 : PR INTCHR$ C 147 V * 

:f * MATRIX CODE * * * ” 

30 PR INT:PR INT"DIESES PROGRAMM VER- UND 
ENTSCHLUESSELT NACHRICHTEN 

40 S =2 : FOP. I = 1 TOS : FOR J = 1 TOS : READAC I,J>:NE 
XTs NEXT 

50 PR INT: INPUT“ENT- ODER VERSCHLUESSELN 

* £ /W ) h ; 

70 I FA$ < ;> " E " ANDA$ < > " V " THEN50 
80 ET = "VER" 

90 IFAf = "E"THENBf="ENT":FORI = 1TOS:FORJ = l 
TOS:READAC I,J):NEXT:NEXT 

100 PRINTCHR$< 147>"NACHRICHT EINGEBEN CM 


AH 254 2EICHEN>":PR INT 
110 M$="":L=0 
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120 GETGf:IFG*=""THEN1E0 

1 30 G = ASCC Gf > 

140 I FG >3 1 ANDG< 3 1 THENMf = M$ + Gf:PR I NTGf : L 
= L+ 1 

150 IFG=£0ANDL >0THENPR INTGf ' :L=L-1 : Mf =LE 
FTf < Mf ,L > 

160 I FG< > 1 3ANDL< S25THEN1 20 

170 PRINT" "sM=INT<L/S+0.9)!DIMB<S,M>,C< 

S ,M> 

180 REM NACHRICHT VERXENTSCHLUESSELN 
130 FORI = 1TOS s FORJ =1 TOM 
200 K = S * J +I -S 

2 10 IFK<=LTHENN=ASC< MIDf < M$ , K , 1 > > -3 1 : B< I 
,J >=N-59 * INT<N/59> 

220 NEXT:NENT 

230 FOR I =2TOS: IFB< I , MX 1 THENBt I,M) = 1 
£40 NEXT 

£50 FOR1 = 1 TOS!FORJ = 1 TOM 

£60 F ORK = 1 TOS : C< I,J>=C< I,J>+A< I,K>*B<K,J 
>:NEKT 

£70 C< l,J>=C-< I , J >-59*INT< C( I,J>/59> 

£80 NEXT s NEKT 

£30 PR INTCHRfC 147 > "D IE "Bf "SCHLUESSELTE 
NACHRICHT LAUTET":PRINT 
300 FORJ =1 TOM•FOR1 = 1 TOS 

3 10 PR INTCHRf< 3 1 +C< I , J > >; 

320 NEXTJNEXT sPR INT 

330 PR INT:PR INT:PR INT"NEUER START <J/N> 

340 GETGf:IF Gf<>"J"ANDGf<>"N"THEN340 
350 IF Gf="J"THENRUN 

360 PR INTCHRf< 147J “AUF WIEDERSEHEN":END 
400 DATA £,- 1 ,-5,3 
4!0 DATA 3,1,5,2 

Für diesen Code haben wir eine 2-mal-2-Matrix zum Verschlüsseln und die Inverse 
der Matrix zum Entschlüsseln gebraucht. Grundsätzlich kann man eine beliebige 
Matrix und ihre Inverse verwenden, sofern beide Matrizen nur ganze Zahlen als 
Elemente besitzen. Hier ist eine weitere Matrix, die zum Verschlüsseln verwendet 
werden könnte: 
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0 1 

1 -2 

Die zugehörige Entschlüsselungsmatrix sieht so aus: 

2 1 

_1 0 _ 

Dieses Verschlüsselungsverfahren könnte man verfeinern, indem man 3-mal-3- 
Matrizen nimmt. Die Botschaft müßte dann in drei Zeilen geschrieben werden und 
das oben beschriebene Verfahren entsprechend übertragen werden. Hier sind zwei 
Matrizen, die man zum Ver- und Entschlüsseln verwenden kann: 


1 

o 


i 2 —r 

2 1 3 


-2 -9 5 

4 2 5 


1 

r— 

1 

C\l 

O 

_1 


Das untenstehende Programm benutzt diese Matrizen zum Ver- und Entschlüsseln. 
Wenn Sie wollen, können Sie Ihre eigenen Matrizen einsetzen. Prüfen Sie aber 
nach, ob sie invers zueinander sind und nur ganzzahlige Elemente enthalten. 

10 REM 38MRTRIXCODE 3 

20 POKE53280,6:PR INTCHRf( 147 >" * 

*+ MATRIX CODE ***" 

30 PR INT SPRINT"DIESES PROGRAMM VER- UND 
ENTSCHLUESSELT NACHRICHTEN 

40 S =3;FDP 1 = 1 TOS:FORJ =1 TOS sREADA< I,J>:NE 
XTsNEMT 

50 PR INT: INPUT"ENT- ODER VERSCHLUESSELN 
(. E/v > " ;a$ 

70 I FAf < >"E"ANDAf < >"V"THEN50 
80 Bf = "VER" 

90 IFAf = "E"THENBf = "ENT" :FOR I = 1 TOS:FORJ = 1 
TOS:REAOA<I,J):NEHT:NEXT 

100 PRINTCHRfC 147>"NACHRICHT EINGEBEN <M 
AM 254 ZEICHEN>":PR I NT 

1 10 Mf ... =0 

120 GETGf: IFGf = " "THEN120 
136 G = ASC< Gf > 

140 IF6>31ANDG<9ITHENMf=Mf+Gf:PRINTGf;:L 
= L+ 1 


129 



150 IFG=20ANDL>0THENPRINTGf ;:L=L-1:M$=LE 
FTt( Mf ,L > 

160 IF G < >13ANDL<££5THEN1£0 

170 PRINT..L/S+0.9).S,M) ,C< 

S , M > 

180 REM NACHRICHT VER/ENTSCHLUESSELN 
190 FOR I = 1TOS:FORJ =1 TOM 
£00 K“S*J+I-S 

£10 I FK< =LTHENN=ASC< MI Df- < M$ , K , 1 > > -3 1 : BC I 
,J>=N-59*INT<N/53> 

££0 NEHT:NEHT 

£30 F0RI=2T0S: IFB< I .. M ) < 1THENBC I,M>=1 
£40 NEKT 

£50 FOR1 = 1 TOS:FORJ =1 TOM 

£60 FORK = 1 TOS : C( I , J > =C< I , J > +A< I ,K ) *B< K , J 
> : NE.KT 

£70 C< I , J > =C< I ,J>-59*I NT< C< I , -T > 759 > 

£30 NEKT:NEHT 

£90 PRINTCHRf< 147>"DIE "Bf"SCHLUESSELTE 
NACHRICHT LAUTETPR INT 
300 FORJ=1TOM:FORI=1TOS 
310 PR INTCHRf < 3 1 +C'l I ,J >> J 
3 £0 NEKT : NE XT : PR I NT 

330 PR INT:PR INT:PR INT"NEUER START <J/N) 

9 II 

343 GETGf: IF Gf< >"J"ANDGf < >"N"THEN340 
350 IF Gf = “J”THENRUN 

360 PR INTCHRf < 147 > "AUF LIIEDERSEHEN" : END 

400 DATA 1 ,0,-1 ,£, 1 ,3,4,£,5 
410 DATA 1,2,-1,-£,-9,5,0,£,-l 


OFFENE SCHLÜSSEL 

Die im vorigen Abschnitt beschriebenen Codes haben einen Nachteil: Wenn man 
einmal weiß, wie eine Nachricht zu verschlüsseln ist, weiß man auch, wie sie zu 
entschlüsseln ist. Das gilt nicht für offene Schlüssel. Sie bestehen aus zwei Teilen: 
dem Verschlüsselungsverfahren, das veröffentlicht wird und jedermann das Ver¬ 
schlüsseln von Nachrichten ermöglicht, und dem Entschlüsselungsverfahren, das 
geheimgehalten wird und nur dem Hersteller des Codes das Entschlüsseln von 
Nachrichten erlaubt. 


130 



Im folgenden beschreiben wir einen solchen offenen Schlüssel. Wählen Sie 
zunächst zwei sehr große Primzahlen P und Q. Jede sollte ungefähr 50 Dezimalstel¬ 
len haben; deshalb läßt sich die Suche danach kaum mit dem Commodore 64 
realisieren. Sei N das Produkt von P und Q. Wählen Sie nun eine ganze Zahl A, die 
kleiner als N ist und keinen gemeinsamen Faktor mit (P— 1)*(Q — 1) hat. Die Zahlen 
N und A dürfen dann jedem zugänglich sein. 

Wie wird eine Nachricht mit den Zahlen N und A verschlüsselt? Man geht dabei so 
vor: 

1. Übersetze die Nachricht in Zahlen (Leerstelle = 01, A = 34, etc.); sie ist dann 
eine große Zahl. 

2. Zerlege die in eine Zahl übersetzte Nachricht in Blöcke übersichtlicher Länge. 

3. Verschlüssele jeden Block B wie folgt: 

C = BfA—N*INT((BfA)/N) 

Die Zahl B wird zur Potenz A erhoben, und Vielfache von N werden weggelassen, 
um eine Zahl zwischen 0 und N zu erhalten. 

Wie wird die entstandene Nachricht wieder entschlüsselt? Da der größte gemein¬ 
same Teiler von A und (P —1 )*(Q —1) die Zahl 1 ist, gibt es zwei Zahlen X und Y, so 
daß A*X + (P —1 )*(Q —1 )* Y = 1 ist. Mit Hilfe von X kann man die Nachricht wie folgt 
entschlüsseln: 

1. Zerlege die verschlüsselte Nachricht in Blöcke. 

2. Berechne für jeden Block C: 

C|X —N*INT((CtX)/N). 

3. Füge die entstandenen Blöcke wieder zusammen und übersetze die Zahlen 
wieder in Zeichen zurück. 

Das Verfahren funktioniert, weil 

(BjA)fX = Bj(A*X) 

= B|(1—(P—1)*(Q-1)*Y) 

= B + Vielfache von N. 

Die letzte Aussage folgt aus einem Satz, der von dem Mathematiker Fermat 
bewiesen worden ist. 

Warum ist dieser Code schwer zu knacken? Beachten Sie, daß wir zum Entschlüs¬ 
seln der Nachricht den Wert von X kennen müssen. Dieser Wert kann aus dem Wert 
von (P—1)*(Q-1) berechnet werden. Um (P-1)*(Q-1) zu kennen, müssen wir P 
und Q kennen. Bekannt sind nur A und N. Wenn man N kennt, ist es theoretisch 
möglich, N zu faktorisieren und P und Q zu bestimmen. Die Zahl N hat aber 
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ungefähr 100 Dezimalstellen, und es sind ungeheure Rechenzeiten (mehrere 
Millionen Jahre) erforderlich, um solche Zahlen zu faktorisieren. Aus diesem Grund 
ist der Code sicher. 

Um einen sehr sicheren Code herzustellen, brauchen Sie jetzt nur noch ein 
Programm für Ihren Commodore 64 zu schreiben, das mit sehr großen Zahlen 
fehlerfrei umgehen kann. Siehe dazu das Kapitel ‘Vermischtes’. 
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KAPITEL 13 


ZUFALL! 

KOPF ODER ZAHL 

Wirft man eine Münze, so sind Kopf oder Zahl als Ergebnis gleich wahrscheinlich, 
zumindest, wenn die Münze nicht verfälscht ist. Ihr Commodore 64 kann das 
Werfen von Münzen mit der eingebauten Funktion RND simulieren. 

Für jede positive Zahl X liefert RND(X) eine (Pseudo-)Zufallszahl zwischen 0 und 1 
(die 1 nicht eingeschlossen). Wenn Sie Ihren Commodore 64 einschalten und 
folgendes tippen: 

FOR 1 = 1 TO 5:PRINT RND(1 ):NEXT I 

erhalten Sie diese Zahlenfolge: 

.185564016 

.0468986348 

.827743801 

.554749226 

.897233831 

Falls Sie Ihren 64er aus- und wieder einschalten (nicht empfehlenswert) und die 
Anweisung wiederholen: 

FOR 1 = 1 TO 5:PRINT RND(1 ):NEXT 1 

so erhalten Sie dieselbe Zahlenfolge. Um das zu vermeiden, beginnt man gewöhn¬ 
lich ein Programm, in dem Zufallszahlen Vorkommen, mit einer Zeile wie 

Y = RND(-TI) 

die dazu führt, daß eine neue Folge von Zufallszahlen ausgelöst wird. Allgemein gilt. 
Für negatives X beginnt RND(X) eine neue Folge von Zufallszahlen. Ist X — 0, so 
ergibt sich dieselbe Zahl wie vorher. 

Das folgende kurze Programm simuliert das Werfen einer unverfälschten Münze. Es 
wird eine Tabelle ausgedruckt, die zeigt, ob Kopf (K für Kopf) oder Zahl (Z für Zahl) 
erscheint. Nachdem 100 solcher Buchstaben ausgedruckt sind, wird angezeigt, wie 
häufig Kopf und Zahl erschienen sind. 
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10 REM 39KQPF ODER ZAHL 

20 POKE53280,6sPRINTCHRJY 147)" *** 

KOPF ODER ZAHL ***“ sK=0 
30 PR INTsPR INT"DIESES PROGRAMM SIMULIER 
T DAS WERFEN EINER MUENZE "} 

40 PRINTsPRINT!INPUT"ANZAHL DER WUERFE " 

.sN 

50 I FN< 1 ORN< > INT< N> THEN40 
60 Y = RNDC -TI) 

70 K=K+1s J=0:PRINTCHR$< 147 >STR$< K >". DUR 
CHLAUF MIT "N"WUERFEN"CHR*< 158>sPR INT 
80 FOR1 = 1 TON 

90 A*="K"!IF RNDC 1>>=.5THENA*="Z":J=J+1 
100 PRINT" "A*; 

1 10 WEHT 

120 PR INT s PR INT:PRINTCHRiT 154>N "WUERFE" 

130 PRINTsPRINT” DAVON "N-J" KOEPFE, "J 
H ZAHLEN" 

140 PRINTSPRINTSPRINT" NEUER DURCHLAUF < 

J/N> ?" 

150 GETG* s I FG$< >"N"ANDG$ < >"J"THEN150 
160 I FG$ = " J 11 THEN200 

170 PR I NTCHR* < 147 > " AUF WIEDERSEHEN" s END 
200 PRINTCHR*< 147>"ANZAHL DER WUERFE VE 
RAENDERN ( J/N) ?" 

210 GETG*: IFG*< >"J"ANDG$< >"N"THEN210 
220 IFG* = "J »THENRUN 
230 GOTO60 

Geben Sie statt Zeile 60 diese Zeile ein: 

60 Y = RND(-I) 

Nun müßten Sie feststellen, daß Sie bei jedem Programmdurchlauf dieselbe Folge 
von Kopf und Zahl erhalten. Deshalb verwenden wir 

60 Y = RND(-TI) 

um die Folge von Zufallszahlen zufällig zu machen. 
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WÜRFELN 


Wird ein unverfälschter sechsseitiger Würfel geworfen, so erscheint eine der sechs 
Zahlen 1, 2, 3, 4, 5 und 6, keine häufiger als eine andere. Wiederum kann der 
Commodore 64 benutzt werden, um das Würfeln zu simulieren. Das nächste 
Programm listet das Ergebnis von 120 Würfen auf. Am Schluß wird ausgedruckt, wie 
häufig jede Augenzahl vorgekommen ist. 

10 REM 40WUERFEL 

20 POKE53280,6 s PRINTCHR$( 147)" 

*** WUERFEL ***":K=0 

30 PR INT!PRINT"DIESES PROGRAMM SIMULIER 
T DAS WERFEN EINES WUERFELS" 

40 PR INT:PR INT! INPUT"ANZAHL DER WUERFE " 

;n 

50 IF NC 1OR N< > I NT< N > TH E N4 0 
60 Y = RNOC -TI) 

70 K =K + 1:FOR1 = 1T06:A< I> =0:NEXT 
30 PR INTCHRf ■: 147 >STR$< K > " . DURCHLAUF MIT 
”N”WUERFEN"CHR$< 158>:PR INT 
30 FOR1=1TON 

100 L= INT< RND< 1 > *6 > + 1 : PR INTL " ";:A<L>=A< 

LD + l 

1 10 NEXT 

120 PR INT:PR INT:PR INTCHR4 C 154 >" BEI "N"W 
UERFEN FIEL" 

130 PR INT:FORI = 1T06:PRINTAC I)"MAL DIE“I, 

: NEXT 

140 PR INT:PR INT:PR INT" NEUER DURCHLAUF < 

J/N) ?" 

150 GETG$ : IFG$< >"N"ANDG*< >"J"THEN150 
160 IFG# = "J"THEN200 

170 PRINTCHR4< 147>"AUF WIEDERSEHEN":END 
200 PR INT:PR INT" ANZAHL DER WUERFE AENDE 

RN <J/N> 7" 

2 10 GETGf: IFG*< >"N"ANDG$ < >"J"THEN210 
220 IFGi="J"THENRUN 

230 IFG$ = "N"THEN60 
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Als nächstes folgt eine andere Version des Programms. Es liefert auch eine 
graphische Darstellung der Würfelseiten. Die Gestaltung der Seiten ist im Feld 
D$(I,J) enthalten. 

10 REM 41UIUERFEL MIT 

20 PQKE53880,6 : PRINTCHRIY 147 > " *** 

FALLENDER WUERFEL ***" 

30 PRINT:PR INT"DIESES PROGRAMM SIMULIERT 
DAS FALLEN EINES WUERFELS" 

40 FOR I =0TO6 : FOR J =0TOS : READA ’■ A$ < I > =A4( I ) 

+ CHR* (A):NEKT:NEKT 

50 FORI=1TOS:FORJ=0TO6:READA:D$C I,J>=A$< 

A $ :NEKT:NEKT 

SO PR INT:PR INT:PR INT:PR INT"BITTE LEER-TA 
STE DRIJECKEN" 

70 GETGt: IFG$< >" "THEN70 
80 Y=RND<-TI) 

90 K= K + 1:PR INTCHR*< 147)K " . WURF"CHR$< 158 
>:PRINT 

100 L -• I NTt RND< 1 > *8 > + 1 : Ar L > = A( L > + 1 

110 FOR J =0TO7 : PR INTTAB't 16) D$ < L , J > : NEKT 

115 PR INT:PR I NT:PR INT:PR INTCHRf< 154 >"IN" 

K" WUERFEN FIEL" 

ISO FR INT:FORI = 1T02:PR INT"DIE"I;AC I ) "MAL 
: NEKT 

185 PR I NT:FOR I =3T04:PR INT"DIE"I;AC I:■ "MAL 
" : NEKT 

130 PR INT:FORI=5T06:PR INT" D I E " I ' AC I )"MAL 
: NEKT 

140 PR INT:PR I NT:PR INT:PR INT"NEUER WURF ( 

J/N) 7 " 

150 GETGf: IFG$< >"J"ANDGS< >"N"THEN150 
160 IFG4="J"THEN80 

170 PRINTCHR$C 147VAUF WIEDERSEHEN " : END 

809 REM DATEN FUER ' S AUSSEHEN 

£10 DATA 213> 192,192,192,192,192,20 1 

820 DATA 881,32,32,32,32,32,221 

230 DATA 221,113,38,38,32,32,221 

240 DATA 221,32,32,113,32,32,221 

250 DATA 221,32,32,32,32,113,221 

2S0 DATA 221,113,32,32,32,113,221 
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£70 DAT« 202,192,192,192,192,192,203 

300 REM DATEN FUER MUERFEL SEITE 

3 10 DATA 0,1, 1,3,1,1,6 

320 DATA 0,4,1 , 1 , 1 ,2,6 

336 DATA 0,4,1 ,3,1 ,2,6 

340 DATA 0,5,1,1, 1,5,6 

350 DATA 0,5,1 ,3,1 ,5,6 

3S0 DATA 0,5,1,5,1,5,6 

Gleichzeitiges Werfen zweier Würfel kann ebenso leicht simuliert werden. Die 
mögliche Augenzahl bei jedem Wurf ist eine Zahl von 2 bis 12. Wie Sie sicher 
wissen, kommen manche Augenzahien häufiger als andere vor. Das sollte auch mit 
jedem der beiden nächsten Programme deutlich werden. Werden zwei Würfe! 
mehrmals geworfen, so ergibt sich die erwartete Häufigkeit (oder die Wahrschein¬ 
lichkeit) jeder Augenzahl aus dieser Tabelle: 

Augenzahl Wahrscheinlichkeit 

2 1/36 

3 2/36 

4 3/36 

5 4/36 

6 5/36 

7 6/36 

8 5/36 

9 4/36 

10 3/36 

11 2/36 

12 1/36 

10 REM 42 2MUERFEL 

20 P0KE53280,6 ! PR I NTCHR4( 147 V *** 

2 ME I MUERFEL ***" 

30 PR INT:PR INT"DIESES PROGRAMM SIMULIER 
T DAS MERFEN VON 2MEI MUERFELN":K=0 
40 PR INT:PR INT: INPUT"MIEVIELE MUERFE "'N 

50 IF N< 10RN< >INT < N>THEN40 
60 Y=RND<-TI> 

70 K=K+ 1 : PR INTCHR$< 147 >K. DURCHLAUF MIT 
"N"WUERFEN":PRINT 
80 FOR1=1TON 
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30 L = I NT< RND< 1 > *6 ) + 1 : M= I NT( RNDC 1 ) *6 > + 1 : ft 
< L +M-2)= AC L+M-2> + 1 

100 PR I NTCHR$ (. 150 > " + " CHR4 < 158 > L ; MCHR* < 15 

0 > " * " j 

110 NE X T 

115 PR INT:PR INTCHR* < 154 > 

120 PRINT" IN"N"WUERFEN FIEL"SPRINT 

130 FORI=0TO 1 0 

140 PR I NT PK I > " X " < I +2 > , 

150 NEXT 

160 PR INTs PR INTsPR INT" NEUER WURF ( J/N> 

n II 

170 GETSf : I FGt < > " N " ANDG-f < > " J " THEN170 
180 IFG$ = "N"THENPRINTCHR$< 147 ) "AUF WIEDE 
RSEHEN":END 

130 PR INTSPRINT" ANZAHL DER WUERFE VERAE 
NDERN <J/N> 7" 

200 GETG$ s I FGf < > "N 11 ANDG4 < > " J " THEN200 
2 10 I F G -f = "N"THEN60 
220 IFG$="J"THENRUN 


10 REM 43 2WUERFEL MIT 

20 POKE53230,6 sPR I NTCHR4- < 147 ) " RMU* * * 2 WE 
I FALLENDE WUERFEL *** " 

30 PRINTS PRINT"DIESES PROGRAMM SIMULIER 
T DAS FALLEN VON ZWEI WUERFELN" 

40 FORI=0TO6 s FORJ=0TO6 s READA s A#< I)=A$< I ) 
+ CHR$ CA >sNEXTsNEXT 

50 FOR 1 = 1 TOS s FOR J =0TO6 s READA s D$ ( I,J)=A'I( 
A >:NEXTs NEXT 

60 PR INTSPRINTsPR INTSPRINT"BITTE LEER-TA 
STE DRUECKEN" 

70 GETGJ6IFG$<>" "THEN70 

SO Y=RND<-TI> 

90 k =K+1 : PR INTCHRf< 147 >K " . WURF"5 PR INT sP 
RINT:PR INTCHR*< 153) 

100 L = I NTC RNDC 1 > +6 5 + 1 : M= I NT< RND( 1 J *6 > + 1 s 
A<L +M-2)=A<L + M-2> + 1 

110 FORT = 0TO7:PR INTTABt 10>D$< L , J>SPC<6 ) D 
$< M,J)s NEXT 
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120 PR I NT:PR INT:PR INTCHRt < 154 > "BE I"K"WUE 

RFEN FIEL"SPRINT 

130 FORI=0TO10 

140 PR INTA< I > " * I +2 > , 

150 NE KT 

160 PR INT SPRINTSPRINT"NEUER WURF < J/N) ? 

II 

170 GETG*:IFG*<>"N"ANDG*<>"J"THEN170 
180 IFG$ = "J"THEN80 

180 PRINTCHR$< 147)"RUF WIEDERSEHEN":END 

£00 REM DATEN FUER'S AUSSEHEN 

£10 DATA £ 13,192,192,19£ , 19£ , 19£ ,20 1 

220 DATA 221 ,32,32,32,32,32,221 

230 DATA 2£1,113,32,32,32,32,221 

240 DATA 22 1 32,32,113,32,32 , £2 1 

250 DATA 221 - 32,32,32,32,1 13,221 

260 DATA £2 1,1 13,32,32,32,1 13,221 

270 DATA 202,192,192,192,192,192,203 

300 REM DATEN FUER UUERFEL SEITE 

3 10 DATA 0,1, 1,3,1,1,6 

320 DATA 0,4,1,1, 1,2,6 

330 DATA 0,4,1 ,3,1 ,£,6 

340 DATA 0,5,1,1,1,5,6 

350 DATA 0,5,1,3,1,5,6 

360 DATA 0,5,1,5,1,5,6 


KARTENSPIELE 

Ein normales Kartenspiel hat 52 Karten. In einem gut gemischten Spiel kann jede 
Karte mit gleicher Wahrscheinlichkeit oben liegen. Das nächste Programm zeigt, 
wie der Commodore 64 das Ziehen einer Karte aus einem gut gemischten Karten¬ 
spiel simulieren kann. Bei jedem neuen Ziehen wird angenommen, daß die vorher 
gezogene Karte zurückgelegt und das Spiel wieder gut gemischt worden ist. 

10 REM 44 KARTEN 

£0 POKE53H80,15:P0KE53£81,15!PR INTCHR*< 1 
47)CHR*(31)" *** KARTEN ***" 

30 PR INT:PRINT"DIESES PROGRAMM SIMULIER 
T DAS ZIEHEN EINER KARTE AUS EINEM "; 

35 PR INT"GUT GEMISCHTEM HAUFEN" 
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50 DIM A$ < 12)s FORI= 0TO12 : READA$ < I } s NEXT 
60 FORI =0TO3 ! READft s B$C I } =CHR* < A } : NEXT 
70 FOR I =0TO3 : READA : CIY I) =CHR*< 1=0 s NEXT 
80 PR INT!PR INT:PR INTSPRINTSPRINT"BITTE L 
EER-TASTE DRUECKEN" 

90 GETG*:IFG*<>" "THEN90 
100 Y = RND(-TI ) 

110 K =K + 1 : PR INTCHR$ C 147 > K "2LIG " : PR I NT: PR I 
NT 

120 L = I NT< RND't 1}* 13 >:M=I NTC RND< 1 > *4 > 

130 PRINTTABC 16 + ( L=9>)C$( M) BIT M>" "A* 

< L } 

140 PR INT:PR INT:PR INT:PR INT s PR INTCHRf ( 31 
>" NE U E S MISCHEN CJ/N) ?" 

150 GETGfs IFGf < >"J"ANDGS< >"N"THENI50 
160 IFG$ = “1"THEN100 

170 P0KE53231,6:POKE53280,14:PRINTCHRtC 1 
4 7> CHR$ < 154} "AUF WIEDERSEHEN":END 
200 REM DATEN 

210 DATA A,2,3,4,5,S,7,8,9,10,B,D,K 
220 DATA 120,122,115,97 
230 DATA 144,28,28,144 

Im Programm ‘Karten’ wählt der Computer erst eine Zahl zwischen 0 und 12, die 
(um 1 erhöht) den Wert der Karte bestimmt. Dann wählt er eine Zahl von 0 bis 3, die 
für eine der vier möglichen Farben der Karte steht. Statt dessen wäre es auch 
möglich, eine Zahl zwischen 0 und 51 zu wählen: 

K = INT(RND(1 )*52) 

und Wert und Farbe der Karte aus dieser Zahl zu bestimmen. Das leistet diese 
Zeile: 

L = INT(K/4):M = K - 4*L 

Jetzt bestimmt die Zahl L den Wert der Karte, während M die Farbe angibt. 

Das Programm ‘Karten’ zieht eine Karte aus einem gut gemischten Spiel. Vor jedem 
neuen Ziehen wird die Karte zurückgelegt und das Spiel wieder gemischt. Was 
wäre, wenn wir nun das Spiel einmal mischen und dann die Karten der Reihenfolge 
nach von oben nach unten auflisten wollten? Dazu wird eine andere Routine 
benötigt. Im wesentlichen numerieren wir die Karten von 0 bis 51 und ordnen diese 
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52 Zahlen dann zufällig um. Dieses Umordnen geschieht systematisch. Zunächst 
wird die erste Zahl zufällig mit einer der 51 anderen vertauscht. Dann wird die zweite 
Zahl mit einer der 50 übrigen Zahlen vertauscht, und so weiter. Das Programm 
'Kartenmischer’ veranschaulicht die Methode. 

10 REM 45KARTEN MISCHER 

20 P0KE53281,15:POKE532S0,15:PRINTCHRSC 1 
47 ) CHRJ(31) " *** KARTEN MISCHER ** + 

M 

30 PR INT SPRINT"DIESES PROGRAMM ILLUSTRIE 
RT DAS MISCHEN VON KARTEN" 

50 DIMfi$< 12 >:FORI=0TO12:READAf< I>:NEXT 
60 FORI=0TO3s READA s B*< I>=CHR$ < A):NEXT 
70 FORI=0TO3:READA:C$< I>=CHR*C A>:NEXT 
90 D I MD MC 51)! FOR I =0TO5 1 : DYS I ) = I s NEXT 
30 PR INTsPRlNT «PR INT:PR INT"BITTE LEER-TA 
STE DRUECKEN" 

100 GETG4 : IFG *<f " "THEN100 

119 Y=RNDC-TI> 

120 K=K+1:PRINTCHR$< 147>:PRINTK". MISCHE 
M":PR INT:PR INT 

130 REM MISCHEN 

140 FORI=0TO50 

15 0 L = I NT < RND( 1 > *< 52 - I > ) + I 

1 70 T = D YS I > : OYS I ) =D/Ä L ) t OYS L > =T 
180 NEXT 

130 FORI=1T052 

£00 L = INTC DYS I - 1 > /4 > : M=DX< I - 1 ) -4 *L 

2 10 PRINTTABC O >C$< M)B$< MD" "A*<L>;:0=0 

+ 10! IF I NT< I Y4 ) = < I /4 > TH E NP R I NT: 0 = 0 

220 NEXT 

£30 PR INT:PR INT:PR INTCHR* ( 3 1)" NEUES MIS 
CHEN < J./'N) ?" 

£40 GETGf:IFGf<>"J"ANDGf<>"N"THEN240 
250 lFGf="J"THEN120 

£60 P0KE532S1,S:POKE532S0,14:PRINTCHR*< 1 
47 >CHR$< 154)"AUF WIEDERSEHEN":END 
£70 REM DATEN 

30© DATA A , £ , 3,4,5,6,7,3,9,10 , B , D , K 
310 DATA 120,122,115,97 
320 DATR 144,£3,28,144 
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UNTERSCHIEDLICH WAHRSCHEINLICHE EREIGNISSE 


ln den meisten bisher behandelten Beispielen ist jedes mögliche Ereignis ebenso 
wahrscheinlich wie ein anderes. Beim nächsten Beispiel ist das nicht der Fall. 

Ein Topf enthält 100 farbige Knöpfe, und zwar 6 rote, 54 blaue und 40 grüne. 
Nachstehende Zeilen simulieren die Auswahl eines Knopfes aus dem Topf: 

X = RND(1) 

R$ = “ROT” 

IF X> = 0.06 THEN R$ = “BLAU” 

IR X> = 0.60 THEN R$ = “GRUEN” 

PRINT R$ 

Das Programm ‘Knöpfe’ simuliert diese Entnahme von Knöpfen aus dem Topf. 100 
mal wird ein Knopf gewählt, und nach jeder Ziehung wird er wieder zurückgelegt. 

10 REM 4SKNOEPFE 

£0 P0KE5323 1,15 : POKE53S80,15 : F'R INTCHR* ( 1 
47 )CHR$'t 144 > " *** KNOEPFE ***" 

30 PR INT:PR INT"IN DIESEM PROGRAMM WIRD E 
IN KNOPF AUS EINEM SACK MIT 6 ROTEN"; 

40 PRINT”, 40 GRUENEN UND 54 BLAUEN KNO 
EPFEN GEZOGEN" 

SO PR INT:PR INT:PR INT:PR INT"BITTE LEER-TA 
STE DRUECKEN" 

70 GETGf:IF6*<>" "THEN70 

80 V=RND<-TI> 

30 K =K+ 1 :PR INTCHRS< 147)K". DURCHLAUF" 

100 X =RND( 1):R4 =CHR$(£8)+ CHR4 C 203 > + " RO 
T" 

1 10 IFX >=0.06THENR4=CHR*< 31 1+CHR$<£09 > + " 

BLAU " 

1£0 IFX > = 0.S0THENR* =CHR$<30>-+CHR$< £09) + " 

GRUEN" 

130 PR INT:PR INT:PR INTTAB< 16 >R$:PR INT:PR I 
NT : PR INTCHRP't 144)" BE I "K "Z IEHUNGEN" ; 

1 35 PR INT" WURDE":PR INT 

140 S =ASC( R$ )-£8:A<S)=A<S > + 1 

150 PR INTCHR4-: £8 ) A( 0 ) " * ROT " ,CHRf < 3 1 >TAB 

< 1 £ ) A < 3) "+• BLAU”; 

155 D R INTTAB't £6 )CHR$( 30 ) A< £ ) 11 * GRUEN” 
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160 PRINTiPRINTCHR*< 144>" GEZOGEN" 

170 PR INT:PR INT:PR INT" NEUE ZIEHUNG ( J/N 
^ 

i80 GETG*: I FG $< >"N"ANDG$ < >"J"THEN180 
190 I FG$ = "J”THEN90 

200 POKE53280,14:P0KE53281,6:PRINTCHR*< 1 
47 > CHR# C i54 > "AUF WIEDERSEHEN":END 
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KAPITEL 14 


STATISTIK 

Was tut man, wenn man mit einer großen Menge numerischer Daten konfrontiert 
ist? Hier ist etwa eine Reihe von Zahlen, die z. B. die Ergebnisse eines Experiments 
oder einer Prüfung darstellen könnten: 

23 67 89 45 40 10 5 
19 99 40 23 9 11 21 
34 34 56 41 42 27 80 

Man kann einiges über diese Zahlen erfahren, wenn man sich deskriptive Maßzah¬ 
len wie das Mitte!, die Varianz und die Standardabweichung ansieht. Aber zunächst 
muß man die Daten in den Computer eingeben. 

Im Commodore 64 kann man die Daten mit Hilfe eines Feldes X(l) abspeichern, 
wobei X(0) = 23, X(1) = 67, etc. zugewiesen wird. Die Eingabe kann im Dialog mit 
einem einfachen Programm vor sich gehen. Das veranschaulichen die Programme 
'Daten-Eingabe 1 ’ und 'Daten-Eingabe 2'. 

200 REM DATEN-EINGABE 1 

210 POKE53E80,6JPRINTCHRiK 147> " * 

*■* DATEN-EINGABE" 

220 PR INT:PR INT"MIT DIESEM PROGRAMM KOEN 
NEN DATEN IM COMPUTER GESPEICHERT"; 

£25 PRINT" WERDEN" 

230 PR INT:PRINT"WIEVIELE DATEN SOLLEN EI 
NGEGEBEN WERDEN?" 

£40 »RINT: INPUT"WIEVIEL " ;N 

£50 IFN< 20RN< >I NT<N > THENPRINT:PR INT"UNMO 

EGLICH":GOTOS40 

£60 N=N-1:DIMX<Ni> 

270 FDP 1 = 1 TON 

£30 PRINTCHRfC 147) "BITTE -JETZT DIE "N+l" 

DATEN E INGEBEN” : PP. INT 
230 PR INT"BITTE "1 + 1“. WERT EINGEBEN 
INPUTNC I) 

3O0 NEXTI 

3 10 PR INT:PRINT"WEITER MIT LEERTASTE " 

320 GETGf: IFG$ < >" "THEN320 
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‘Data Entry I’ fragt zu Beginn nach der Anzahl der einzugebenden Daten. Das Feld 
X(l) wird dann entsprechend dimensioniert. 'Data Entry II’ nimmt hingegen an, daß 
nicht mehr als 100 Zahlen eingegeben werden. Sie geben die Daten auf Anforde¬ 
rung ein; wenn Sie fertig sind, geben Sie -99999 ein. Falls nötig, kann die Zahl M in 
Zeile 240 von 100 in eine beliebige andere Zahl geändert werden. 

Dieses Kapitel enthält mehrere kurze Routinen, die an die 'Data Entry’-Programme 
angehängt werden können. Insgesamt ergibt sich damit ein nützliches Programm 
zur Auswertung Ihrer Daten. 

200 REM DATEN-EINGABE 2 

210 POK.E53280 , S s PR I NTCHR$ < 147>" * 

*+ DATEN-EINGABE" 

220 PR INT SPRINT"MIT DIESEM PROGRAMM KOEN 
NEN DATEN IM COMPUTER GESPEICHERT"; 

225 PRINT" WERDEN" 

230 PRINTsPRINT"WEITER MIT LEERTASTE " 

240 GETG* : IFG*< >" "THEN240 

250 M= 100 : D IMX< Mi' 

260 FORI= 0TOM 

270 PRINTCHR$< 147)"BITTE DATEN EINGEBEN" 

2S0 PR INTsPR INT"EINGABE MUSS MIT -93999 
BEENDET WERDEN" 

290 PR INT:PRINT"BITTE "I + 1 ". WERT EINGEB 
EN " ; s INPUTXt I > 

300 I FX< I >=-93999AND I >1THENN=1 : I =M 
305 NE 7 T 

3 10 PR INT SPRINTSPRINTSPRINT"WEITER MIT L 
EERTASTE " 

320 GETG$ : IFG$< >" "THEN320 


MITTELWERT 

Der Mittelwert oder Durchschnitt ist eine wichtige statistische Maßzahl. Man erhält 
ihn, indem man alle Zahlen addiert und die Summe durch die Anzahl der Zahlen 
teilt. 

Mittelwert = Summe der Daten 
Anzahl der Daten 
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Sind die Zahlen im Feld X(l) für I = 0 bis N — 1 gespeichert, so kann der Mittelwert 
mit Hilfe der folgenden Programmzeilen berechnet werden: 


X=0 

FOR 1=0 TO N —1 : X=X+X(I) : NEXT 
XM=X/N 

Das unten angeführte kurze Programm berechnet den Mittelwert und druckt ihn aus. 
Es kann in die ‘Daten-Eingabe’- Programme eingebaut werden. 

400 REM MITTELEWERTE 

410 ROKE53S30,6 ■ RR INTCHR'Tf 147 >" *** D 

ATEN-ANALYSE SEITE I * * *" SPRINT 
4 15 PR INT: RR INT"DATEN-ANALYSE VON"N+l"WE 
RTEN " 

420 X=0:FORI=0TON:K=X+K( I ) : NEXT: XM=X/T N+ 

1 > 

430 PR INT:PR INT:PR INT"MITTELWERT 
= "XM 


MAXIMUM, MINIMUM, SPANNWEITE 

Es ist oft sinnvoll, den größten und kleinsten Wert der Daten zu kennen. Sie können 
Ihren Commodore 64 danach suchen lassen. Ein solches Suchprogramm ist unten 
abgedruckt. Es berechnet zusätzlich die Streubreite oder Spannweite der Daten. 
Das ist einfach die Differenz zwischen dem größten und kleinsten Wert der Daten. 

500 REM MAXIMUM,MINIMUM UND SPANNWEITE 
510 MAX = -1E +37:MIN=1E+37 
520 FORI=0TON 

530 IF X< I ? >MA.X TH EN MAX=X< I) 

540 IF X< I)< MIN THEN MIN=X< I> 

350 NEXT 

5S0 PRINT:PR INT;P RINT"MINIMUM 
= " M IN 

570 P RINT:P RINT"MAXIMUM = " MA 

580 PR INT:PR INT"SPANNWEITE = "MA 

X -M I N 
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Vielleicht möchten Sie Ihre Daten auch in ansteigender (oder absteigender) Reihen¬ 
folge ordnen. Dafür gibt es verschiedene Verfahren (wie Bubble-Sort, Quick-Sort, 
Shell-Sort). Wir gehen an dieser Stelle jedoch nicht ins Detail. 


STANDARDABWEICHUNG UND VARIANZ 

Der Mittelwert ist eine einfache, nützliche und aussagekräftige Maßzahl. Aber er 
sagt uns nicht alles, was wir wissen müssen, Sehen Sie sich etwa diese beiden 
Datensätze an: 

DATEN für X(l) 20, 21, 20, 19 
DATEN für Y(l) 38, 26, 14, 2 

Die Mittelwerte XM und YM sind beide gleich 20. Aber die Daten für Y(l) streuen viel 
stärker als die Daten für X(l). Diese Streuung kann man mit der Standardabwei¬ 
chung der Daten messen. 

Die Standardabweichung eines Datensatzes ist durch die Formel gegeben: 

Standardabweichung = Quadratwurzel (Summe der Quadrate der Differen¬ 
zen zwischen Daten und Mittelwert / (Datenzahl 
-1». 

Die Varianz ist das Quadrat der Standardabweichung. Die Standardabweichung der 
im Feld X(l) gespeicherten Daten wird nach diesem Verfahren berechnet: 

1. Berechne den Mittelwert XM. 

2. Bestimme die Abweichungen vom Mittelwert, das-heißt die Werte X(l) -XM. 

3. Quadriere jede Abweichung, das heißt berechne (X(I)-XM) 2 . 

4. Addiere die Quadrate der Abweichungen. 

5. Teile durch die Datenanzahl minus 1. Das ergibt die Varianz XV der Daten. 

6. Ziehe die Quadratwurzel. Das ergibt die Standardabweichung XD der Daten. 

Die Standardabweichung gibt eine Vorstellung davon, wie weit die Daten um den 
Mittelwert streuen. Sehen wir uns noch einmal die Beispiele an: 

DATEN für X(l) 20, 21, 20, 19 
DATEN für Y(l) 38, 26, 14, 2 
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Die Standardabweichungen ergeben sich aus diesen Berechnungen: 

XD = SQR((0*0 + 1*1 + 0*0 + (— 1 )*( — 1 ))/3) 

= SQR(2/3) 

0.816496581 

YD = SQR((18*18 + 6*6 + (-6)*(-6) + (—18)*( —18))/3) 

= SQR(240) 

= 15.4919334 

Die Ergebnisse spiegeln deutlich den Unterschied in der Streuung der Daten wider. 
Das nächste kurze Programm berechnet die Standardabweichung der im Feld X(l) 
abgelegten Daten. 

600 REM STANDARDABWEICHUNG 

S10 v= 0 :FQRI=0TON:Y=X< I>-XM:X=X+Y*Y:NEXT 

620 XD = SQR< X. -'N > 

630 PR INT:PR INT:PR INT"STANDARDABWEICHUNG 
= "XD 


KONFIDENZINTERVALLE 

Die Standardabweichung ist nützlich, weil sie anzeigt, wie stark die Daten um den 
Mittelwert streuen. In vielen Massenherstellungsverfahren variieren die erzeugten 
Produkte leicht in Größe, Qualität, Länge etc. Man nennt die zu messenden Objekte 
die Population. Die Abweichungen einer Population sind häufig normalverteilt. 
Statistiker haben festgestellt, daß die Normalverteilung in vielen Fällen die tatsächli¬ 
chen Daten gut approximiert. Wenn die Anzahl der Daten groß ist (mehr als 
ungefähr 30), dann nimmt man für die Berechnungen oft an, daß die Population 
normalverteilt ist, auch wenn es sich in Wahrheit nicht ganz so verhält. 

Vereinfacht gesagt ist eine Normalverteilung symmetrisch um den Mittelwert, wobei 
der größte Anteil der Population nahe beim Mittelwert und nur ein sehr geringer 
Anteil weit davon entfernt liegt. 

Etwas genauer: Ungefähr 80% einer normalverteilten Population sind weniger als 
die Standardabweichung vom Mittelwert entfernt, und ungefähr 96% sind weniger 
als das Doppelte der Standardabweichung vom Mittelwert entfernt. Die folgende 
Tabelle enthält - noch ausführlicher - die zu verschiedenen Vielfachen der 
Standardabweichung gehörigen Prozentsätze. 
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% der Population Vielfaches der Standardabweichung 
50% 0.6745 

68.27% 1 

80% 1.28 

90% 1.645 

95% 1.96 

95.45% 2 

99% 2.575 

99.73% 3 

Diese Tabelle zeigt, daß 95% der Population weniger als das 1,96fache der 
Standardabweichung vom Mittelwert entfernt sind. Anders gesagt liegen 95% der 
Population im Bereich von 

XM - 1.96*XD bis XM + 1.96*XD, 

wobei XM der Mittelwert und XD die Standardabweichung ist. Man nennt dieses 
Intervall das 95%-Konfidenzintervall für die Population. Ähnlich reicht das 99%- 
Konfidenzintervall von 

XM - 2.575*XD bis XM + 2.575*XD. 

Wir wollen jetzt an einem Beispiel erläutern, wie Konfidenzintervalle angewandt 
werden können. Angenommen, Sie haben den Verdacht, daß Ihr Händler Ihnen 
nicht ganz volle Zweiliterflaschen Limonade verkauft. Sie kaufen zehn Flaschen und 
messen sorgfältig den Inhalt. Das Ergebnis in Litern ist: 

2.001, 2.040, 2.020, 2.000, 2.015 
2.006, 2.005, 2.031, 2.008, 2.018 

Alle Flaschen enthalten mindestens 2 Liter. Lassen Sie uns jedoch Mittelwert und 
Standardabweichung dieser Daten berechnen. Das Ergebnis: 

XM = 2.0144 
XD = 0.0132 

Unter der Annahme, daß unsere Stichprobe aus einer normalverteilten Population 
stammt, können wir einige Konfidenzintervalle berechnen. Das 95%-Konfidenzin- 
tervall reicht von 

2.0144 - 1.96*0.0132 bis 2.0144 + 1.96*0.0132, 
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also von 


1.989 bis 2.040. 

Wir erwarten also, daß 95% der Flaschen zwischen 1.989 und 2.040 Liter enthalten. 
Das heißt, wir erwarten, daß 2,5% mehr als 2.040 Liter enthalten, während 2,5% 
weniger als 1.989 Liter enthalten. Wir können auch das 90%-Konfidenzintervall 
ausrechnen. Es reicht von 

2.0144 - 1.645*0.0132 bis 2.0144 + 1.645*0.0132, 
oder von 
1.993 bis 2.036. 

Also erwarten wir, daß 90% der Flaschen zwischen 1.993 und 2.036 Liter enthalten. 
Folglich würden mindestens 5% der Flaschen weniger als die geforderten 2 Liter 
enthalten. (Entsprechend enthalten 5% der Flaschen mehr als 2.036 Liter.) 

In dieses Beispiel sind viele Annahmen eingegangen; die daraus gezogenen 
Schlüsse wären deshalb keine ausreichende Grundlage für eine Strafanzeige. 

Der Mittelwert unserer Population ist aus einer Stichprobe berechnet worden. 
Woher wissen wir, daß dies der tatsächliche Mittelwert der Population ist? Der 
Mittelwert hängt von der gezogenen Stichprobe ab. Wir können aber mit hülfe der 
Standardabweichung schätzen, wie weit unser Stichprobenmittelwert vom wahren 
Mittelwert entfernt ist. Wir können mit 95% Vertrauen sagen, daß der wahre 
Mittelwert zwischen folgenden Werten liegt; 

XM - 1.96*XD/SQR(N—1) bis XM + 1.96*XD/SQR(N-1), 

wobei XM der aus einer Stichprobe vom Umfang N berechnete Mittelwert ist. 
Dieses Intervall heißt 95%-Konfidenzintervall oder -Vertrauensintervall für den 
Mittelwert. Das 99%-Vertrauensintervall für den Mittelwert reicht von 

XM - 2.575*XD/SQR(N-1) bis XM + 2.575*XD/SQR(N-1). 

Genau genommen sind diese Berechnungen nur gültig, wenn der Stichprobenum¬ 
fang groß ist (sagen wir größer als 30). Für kleinere Stichproben sollten wir 
eigentlich die sogenannte Studentsche t-Verteilung statt der Normalverteilung 
verwenden. Aber das würde den Rahmen dieses Buches sprengen. 

Nehmen Sie die Konfidenzintervalle nicht zu ernst, und verwechseln Sie nicht die 
beiden Arten von Konfidenzintervallen. 
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Das nächste Programm berechnet 95%- und 99%-Vertrauensintervalle für die 
Population und den Mittelwert. 

700 REM VERTRAUENS INTERVALLE 

710 Cl=l.3B*XD:M1=C1/SQR<N):C£=2.575*XD: 

M£=C2/SQR<N) 

720 PR I NT : PP. I NT : PR I NT " VERTRAUENS INTERVAL 
35 X" 

730 PR INT:PR INT"POPULATION VON "XM- 

C 1 

735 PRINT" BIS "XM+Cl 

740 PR INT:PR INT"MITTELWERT ZWISCHEN "XM- 
M1 

745 PRINT" UND "XM+Ml 

750 PRINT:PR INT:PR INT"VERTRAUENS INTERVAL 
99 X" 

7S0 PR INT:PRINT"POPULATION VON "XM- 

C2 

7S5 PRINT" BIS "XM+C2 

770 PR INT:PR INT"MITTELWERT ZWISCHEN "XM- 
M£ 

775 PRINT" UND "XM+M2 

DER LETZTE SCHLIFF 

Verbindet man die Programme dieses Kapitels miteinander, so erhält man ein 
nützliches Programm zur Analyse von Daten. Das nachstehende Programm enthält 
die erforderlichen Verbindungsglieder. 

10 REM 47DATEN-ANALYSE 

20 PQKE53280,6:PRINTCHRf C 1 47 ) " ** 

* DATE N-A NALYSE ***" 

30 PR INT:PRINT"DIESES PROGRAMM FUEHRT EI 
NE ANALYSE VON DATEN DURCH, "i 
40 PR INT"DIE SIE DEM COMPUTER EIN- GEBE 
N. " 

50 PR INT:PRINT"DIE ANALYSE WIRD IN FOLGE 
NDEN PUNKTEN DURCHGEFUEHRT" 

60 PR INT:PR INT"1 > MITTELWERT DER DATEN 
70 PP. INT: PR INT "2 > MAX IMUM,MINIMUM UND SP 
ANNWEITE 
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30 PR INT! PR INT"3) STANDARTABWE ICHUNG 
90 PR I NT:PRINT" 4 ) VERTRAUENS INTERVALLE ( 
95 UND 99 IO 

100 PR INT!PR INT:PR INT"BITTE DIMENSIONIER 
EN SIE JETZT DEN SPEI-CHER. " ; 

110 PR INT"SIE KOENNEN DIE EINGABE JEDER2 
EIT MIT '-99999' BEENDEN." 

1£0 PR INT: INPUT"ZAHL ";M:M=M-1 
1 30 I FM< £ORM< > INT< M > THEN1 £0 
150 REM DATEN EINGABE £ 

160 D I MX < M) 

170 FORI=0TOM 

180 POKE53£S0,6:PRINTCHR$C 147>" #* 

* DATEN EINGABE * **"!PR INT 

£00 PR INT:PRINT"EINGABE VORZEITIG BEENDE 
N MIT -99399 " 

££0 PR INT:PR INT"SPEICHERKAPAZ I TAET TOTAL 
"M+ 1 

£30 PR INT:PR INT" NOCH VORHANDEN 

"M-I+1 

£70 PR INTCHR$< 144):PR INT:PR INT:PR INT:PR I 
NTI+1". WERT 

£80 INPUTXC I >:PRINTCHR*< 154): IFXC I)=-999 
99AND I>1THENN=I- 1: I=M:G0T0310 
£90 NENT 
300 N =I - 1 

3 10 PR INT:PR INT:PR INT"WEITER MIT LEER-TA 
STE " 

3£0 GETG$: IFG#< >" "THEN3E0 

400 REM MITTELEWERTE 

410 POKE53B80,6:FR INTCHR$( 147 ) " *** D 

ATEN-ANALYSE SEITE I * **":PR INT 
415 PR INT:PR INT"DATEN-ANALYSE VON"N+l"WE 
RTEN " 

4£0 X =0 : FOR I =0TON: X =X +NX I > : NEXT: XM=X/< N+ 
1 > 

430 PR INT:PR INT:PR INT"MITTELWERT 
= "NM 

500 REM MAX IMUM,MINIMUM UND SPANNWEITE 
5 10 MAN = -1E +37:MIN=1E +37 
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520 FQRI=0TON 

530 IF X< I>>MAX THEN MAX=X< I> 

540 IF XCIXMIN THEN MIN=X< I > 

550 NEMT 

560 PR INT:PR INT:PR INT"MINIMUM 
= "MIN 

570 PR INTSPRINT"MAXIMUM = "MA 

v 

580 PR INT SPRINT"SPANNWEITE = "MA 

X-MIN 

600 REM STANDARDABWEICHUN6 

610 X=0 s FORI=0TONs Y=X< I>-XMsX=X+Y*Y:NEXT 

620 XD=SQRC X,''N> 

630 PR INTs PR INT:PR INT"STANDARDABWEICHUNG 
= "XD 

640 PR INTCHR$ < 19>sFORI = 1TO20sPR INT:NEXT: 
PRINT 

645 PRINTCHRJC18>"SEITE 1 <E) SEITE 2 < 
2 > SCHLUSS < S ;■ “ 

650 GETGfs IFG$ < >"E"ANDG$< >"2"ANDGF< >"S"T 
HENS20 

655 IFG*="E"THEN400 
660 IFGf = "2"THEN690 
665 IFG$ = "S"THEN900 

690 PQKE53280,6sPR INTCHR$< 147 >" * * * D 

A TE N-ANALYSE SEITE II ***" 

700 REM VERTRAUENS INTERVALLE 

710 C1 = 1.96 *XD:Ml=C1/SQR<N>s C2=2.575«XDs 

M£=C2/SQR<N> 

720 PR INTsPR INTsPR INT"VERTRAUENS INTERVAL 
95 

730 PR INTsPR INT"POPULATION VON "XM- 

C1 

735 PRINT" BIS "XM+Cl 

740 PR INT:PR INT"MITTELWERT 2WISCHEN "XM- 

M1 

745 PRINT" UND "XM+Ml 

750 P RINT s PRINT s P RINT"VER TRAU ENS INTE R VA L 
99 X” 
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760 PR INT:PR INT"POPULftTION VON "XM- 

CS 

765 PRINT" BIS "XM+C2 

770 PRINTJPRINT"MITTELWERT ZWISCHEN "XM- 
M£ 

775 PRINT" UND "XM+M2 

800 PR INTCHR$< 19 > !FOR1 = 1TO20 s PR INT:NEXT '• 
PRINT 

310 PRINTCHR$< 18)"SEITE 1 (E) SEITE 2 ( 
Z > SCHLUSS <S>" 

320 GETGt s IFG$ < >”E"ANDGf< >"Z"ANDG*< >"S"T 
HEN820 

830 IFG$ = "E"THEN400 
840 IFG$ = "Z"THEN690 
850 IFGf = "S"THEN900 

900 PRINTCHR*< 147>:PR INT:PR INT:PR INT:PR I 
NT "NEUER STHRT < J/'N> ?" 

910 GETG*: IFG$< >"J"HNDG*< >"N"THEN910 
920 IFG$="J"THENRUN 

930 PRINTCHR$< 147 VHUF WIEDERSEHEN" : END 
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Mathematik mit dem C 64 


Diskette 


Inhalt: 

"0 1 RECHTU.DREIECK" 
”02DREIECKE" 
"03BRECHUNG" 
"04ERDTRIGON." 
"05POTENZEN 1" 

" 0 6 P 0 TENZ E N 2" 
"07KOMPLEXE ZAHL." 
"08QUAD.GLEICH" 
"09POLYNOME" 

"1©POLYNOME £" 

"11WURZELN" 

"12F0LGEN ARITHM." 
"13F0LGEN GEOMETR" 
■142INSESZINS" 

"15ZINSESZINS VAR" 
"1SFIBONACCI FOL.“ 
"17BASIS UMUANDL." 
" 1SE INH ER UMLIANDL " 
"19LI0CHENTAGE" 
"20KALENDER " 

"21 DATEN MANGE." 
"22GR0ESST.TEILER" 
" 23ERAST0STHENES1 " 
"24PRIMZAHLEN 1" 
"25PRIMZAHLEN 2" 
"2SERAST0STHENES2" 
"27PRIMZAHLEN 3" 
"2SPRIMFAKTOREN" 
"29PRIMZAHLEN 4" 

"30TRIFEL " 

"3 1 AUF 40 STELLEN" 
"32DEM0 MATRIZEN" 
”33 1NYER.MATR IZ EN" 
"34GLEICH.SYSTEM" 
"35EINFACHCODE" 

”3STAUSCHCODE" 

"37MATRIXCODE" 


Alle Programme dieses Ban¬ 
des sind auf der beigefügten 
Diskette gespeichert. Wir 
empfehlen, sofort ein Backup 
anzufertigen und dieses 
sicher zu verwahren. 

Die einzelnen Programme 
können bequem mit z. B. 

LOAD “01*”, 8 

geladen werden. 


"3SMATRIXCODE 3" 
"39K0PF O. ZAHL" 

" 40LJUERFEL " 

"4 1LJUERFEL MIT" 

”42 21-JUERFEL " 

"43 2UUERFEL MIT" 
"44 KARTEN" 
"45KARTEN MISCHER" 
"4SKN0EPFE" 

"47DATEN-ANALYSE" 
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NOTIZEN 




Finden Sie manchmal einzelne 
Teile Ihrer Programme schwierig 
zu schreiben? So daß Sie zwar 
wissen, was Sie wollen, aber sich 
nicht genau im klaren sind, wie Sie 
dabei vorgehen sollen? Kann man 
mit COS, ABS oder SGN etwas an¬ 
fangen? Was bedeuten sie, und 
wie lassen sie sich einsetzen? 

Das vorliegende Buch erklärt 
alle diese mathematischen Hilfs¬ 
mittel. Sämtliche mathematischen 
Funktionen des Commodore 64 
werden darin beschrieben. Ihre 
Anwendungen sind in kurzen Pro¬ 
grammen erläutert, die Sie wörtlich 
in Ihre eigenen Programme über¬ 
nehmen können. 

Das Buch enthält eine Fülle von 
Informationen zu so unterschied¬ 
lichen Gebieten wie Kryptogra¬ 
phie, Zufallszahlen, Logik, Folgen 
und Reihen, Trigonometrie, Prim¬ 
zahlen, Felder und Matrizen, Wahr¬ 
scheinlichkeitstheorie und Stati¬ 
stik. 

Czes Kosniowski lehrt an der 
Universität von Newcastle upon 
Tyne. Er hat mehrere Bücher und 
Artikel über Mathematik und Da¬ 
tenverarbeitung geschrieben und 
ist regelmäßiger Mitarbeiter der 
Zeitschrift Populär Computing 
Weekly. 
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